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ABSTRACT 


Traditionally, reserve estimates of plains coal have not 
included an error analysis, either of the input variables 
(average net thickness, area and density) or of the 
resulting tonnage. Relative sizes of the errors can be used 
to choose between different sources of data or between dif- 
ferent methods of calculation; the better data and method 
Should have a lower error. Density and seam thickness data 
require the use of French style geostatistics to determine 
their errors because they are regionalized variables. 

Data on the production for two years from the Estevan 
zone at the Boundary Dam Mine in southern Saskatchewan was 
studied to determine the sizes and sources of errors in the 
tonnage calculations. The results may be applicable to 
other plains coal mines. 

Net thickness is the main variable studied geostati- 
Stically. Unfortunately, there was not enough density data 
to estimate a variogram, so this variable, along with the 
area data, were analysed using traditional statistics. 

Two sources of thickness data were available: drill 
holes (at 150 m centres) and detailed (30 m spacing) meas- 
urements of the highwall. Two methods of estimating the 
average thickness were considered: ordinary kriging, which 
should be close to the best (=minimum variance) method, and 
the simplest, an arithmetic average. When comparing methods 
uSing pre-production data, the results were compared to the 
"true' tonnage determined from detailed in-pit measurements. 

The variogram models of the thickness were isotropic 
with two nested spherical structures and no nugget effect. 
The drill hole variogram had a higher sill and a longer 
range than the pit meaSurement variogram. 

Density and thickness contribute about equally to the 
error of the tonnage. If the error in the area is con- 
trolled by the amount of ground that can be opened in one 
year, then the error of the area is an order of magnitude 
Hanger than the other errors. Ini this case, 1t does not 
matter what data or what method is used. 

The drill holes produce a biased estimate, low by 6%, of 
the thickness when compared to the pit meaSurements. This 
is because the geological seam was logged in the holes not 
the seam as it would be mined. Better data to use 1s- 
highwall measurements from an opened pit near the projected 
pit; the resulting estimate is unbiased with a slightly 
lower estimation variance. 

The seam is so regular that the estimation errors are 
very small (1% to 5%). This means that the krige variance 
is only slightly smaller than the variance of the arithmetic 
average, which suggests that the extra effort required for 
Gecstatistical estimating is not justified.» In this mine 
geostatistics is only useful for determining the proper 
variance of the thickness. 
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Geostatistical papers are notorious for loose notati 
using a system which is unique to geostatistics. It can 
change several times within one paper, and sometimes the 
Same symbol has two referents; for instance, "a" is the 
weight in a linear combination as well as the range of a 
transitional variogram. Context is often essential to 
unravel the formulas. I have attempted to follow the 
notation used by David (1977) within the limitations of 
Scientific character set on the printer at the Universit 
Alberta. 


This study follows several conventions. 
1. Regular font is used for scalar variables, such as t 
Lagrange parameter u. 
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2. Lowercase bold is used for vector variables, such as the 
sample separation vector h; for instance |h| = h. 

3. Uppercase ITALICS are used for matrices, such as the 
krige matrix ok. 

4, Square brackets [ ] are used to indicate the elements of 
an array or the argument of a function or operation; for 
instance f[x] or E[R]. 

5. Any general variable is represented by x or R as above. 

6. All units are metric, though the mine still uses 
Imperial units. 

The page refers to where the variable is first used or 

described. 
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The size of the support of the data. The 
Support is the size (and shape) of the 


Sample over which the variable 1S averaged 
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calculate a variogram from irregularly 
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Psi is half the width of the angular class 


used to calculate a variogram from 
irregularly spaced data. See Figure 1.5, 


PUSOSy Re Hees hee: o Sere ates ates © Sele cisteteietecletels © 6.6 


Phi is the direction of the distance class 


used to calculate a variogram from 
irregularly spaced data. See Figure 1.5, 
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1. INTRODUCTION 


Matheron (1971, p.5) defines geostatistics as "+*+the 
application of the theory of the regionalized variables to 
the estimation of mineral deposits:*+*" Geostatistical 
estimates of coal reserves have been made for the Boundary 
Dam Mine at Estevan in southern Saskatchewan (Figure 1.1). 
In this mine lignite is stripped from the Estevan Zone of 


the Paleocene Ravenscrag Formation. 


1.1 OBJECTIVES OF THIS STUDY 
The objective of this study is to examine the methods 
of determining reserve estimates in plains coal mines. The 
basic formula is: 
Tonnage = average thickness x area x density 

The overall objective is divided into two major portions. 

1. Study the optimum drill hole spacing. It may be that 
hole separation could be increased for a less expensive 
Sampling program with little loss in accuracy. 

2. Study the methods used to estimate reserves, in particu- 
eta 
A. Choose the better of the two sources of data availa- 

ble on the average seam thickness over an area to be 
mined. The choice is between highwall measurements 
from a nearby pit and drill holes into the proposed 
pits 


B. Test geostatistically-determined estimates to see if 


selfdsitev suki 
: es J 
lanizettaseose atts 
= : a “— <a : 
abaued sf3 sl shear reac oved sev ise 3 ‘Saot 
Ps ee ac 2 
; ae ss 
t,t ssupi® nevada tgAage iiediuee ab a 


vy 
. fa 
728 
7h - 
w - 
+ 
eed 
= 
,=* 
E ) a 
] 2 =. - 
~ “ ~“ 
a ta s 
ne ives 
~ oo a # 
rN * ah J 
4 
~~ (‘~ ts 


font owt Otnd Bebivel at. avisdetde £1 ni 


© 


nittaxe of a2 cone 2 + 6 4 eee 


woos af agol siti djiw ots 


[iewipid asewjed ef sata od? rim 


ie hx .e 
~ A } oe 


3jea 9d9 moa? beqdisse ei mene: 


sae by 


Sl oe 
Astas sot paiseieve os 


: ’ i “7 . 
~ a ae ; é 


PLea? wie other 40 23 


is: 


7) 
p 
bh 
wT 


saenanina 2p S18VS ” 


« 


ti. pri cmge efor, Lita 6 suai igg sda 


10% Seaasvont ed © ives > nok sennaee 


2951 atentine of Beau sbonsae « 


<a 


¢ 
are 7 


2 ows afd Yo. aesded/ ex mee i> 
avg @eoniolda mose sos vers eit. om 


3) 


a. mae 
awe Bas 


“y 


if aetod Ifib bas’ sig | vmod 8 1092 


Saskatchewan 


Manitoba 


Saskatoon 
@ 


Regina 
& 


ESTEVAN 
@ 


0 100 


kilometres 


Figure 1.1 Location map 


Saag 


» 


KAVET2S 
.s 

or - oo 

et: 


ssi2zamoiia 


400 taae ne 


they are useful or better than the traditional 
methods. Geostatistical estimates are theoretically 
Supposed to be the best, but usually at a 
considerable cost in effort. There may be little 
accuracy lost in using simpler methods, simpler even 
than the methods presently used at the mine. 

C. Determine the size of the errors in the tonnage 
estimates and recovery estimates. This will indi- 
cate how much confidence to place in the numbers, 
how detailed to make any future calculations, and 
whether a given change in the recovery is real or 
FUSt a Statistical, fluctuation. 

D. Outline the major sources of the error in the 
tonnage figures. This should show where to concen- 
trate effort if the accuracy of the estimates is to 
be refined. If the error of the reserves is 
dominated by the error of one input variable, Say 
‘area’, then improving the accuracy of that one will 


increase the accuracy of the estimate the most. 


Determining the accuracy and errors above requires 
statistical analysis and because density and thickness are 
regionalized variablest, geostatistical analysis is 
required. All geostatistics must begin with variogram anal- 
ysis, so a preliminary objective is to determine the experi- 
mental variogram and to model it. 


+ See Appendix 1.1.1 for a definition of a regionalized var- 
lable. 
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1.2 OUTLINE 

A complete geostatistical study is done in several 
stages; beginning with a crude preliminary study, making 
many assumptions; finishing by refinements, correcting for 
each of the assumptions made in the first stage. This study 
is only the first stage, but the results were clear enough 
that the objectives were achieved without having to correct 


for all of the assumptions. 


Tonnage (T) is calculated from the averaget net 
thickness (Z)+, density of the coal (d) and area of the pit 


(A) using: 
T = Z-d-A caleae) 


Recovery (R) is calculated from the tonnage and the coal- 


loaded-out (CLO) using: 


Re= 7 100-CLO/T fio! 


The main type of data is thickness measurements of the 
seam, which have two sources: daily measurements in the pit, 
and drill hole measurements. Most calculations were done 
with both to determine the best data to use. The other 
input variables (density, area and coal-loaded-out) are from 
+ 'Average' in this study, is rarely used the traditional 
Sense of an arithmetic mean, but in the more general meaning 
of the 'best' single number with which to replace all the 


numbers in a subset of the data. ; 
+ Traditionally in geostatistics Z refers to a regionalized 


variable. 
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the mine records. 


Only thickness was analysed geostatistically. Although 
density is regionalized, there is not enough data to esti- 
mate the variogram. Instead ordinary (or classical) 
Statistics were used to analyse the density, and also the 
area and coal-loaded-out. Geostatistical analysis along 


with tonnage estimation form the bulk of this study. 


Any statistical analysis begins with exploration of the 
data, that 1S, estimating the basic statistical properties 
of the data. For the thickness, these include drift and 
variogram analysis because geostatistics is used. Variogram 
analysis begins with the calculation of the raw variogram, 
correcting it to form the experimental variogram, which is 
modelled. The variogram model should closely estimate the 
true variogram that is required to determine the estimation 


variances. 


Optimum drill hole spacing was considered in this study 
only as a trade-off of the expense of more closely spaced 
holes for a more accurate reserve estimate. In practice, 
other considerations also affect the choice of hole spacing. 
The trade-off was illustrated with a curve of the estimation 
variance against the spacing in the drilling grid. This 
curve shows at what spacing (70 to 220 m) there is a trade- 


off with the reserve accuracy. The particular balance 
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between the two within this range must be left to the 
judgement of the mine engineer. The curve was calculated 
for a representative pit, approximately two year's produc- 
tion, and using a square grid of drill holes that was 


expanded for each succesSive point on the graph. 


Coal tonnages, and in some cases recoveries, were 
estimated in several ways for several mined areas represent- 
ing production for one month to two years. Tonnage and the 
resulting recoveries are most strongly affected by the aver- 
age net thickness, on which most of the effort was spent. 
Tonnages for the area were calculated using the thickness 
determined in two ways: kriging, the best method from a geo- 
Statistical point of view, and an arithmetic average, the 
Simplest method. Whenever an estimate was made, its 
variance was determined in order to compare the accuracy of 
the different methods. Finally, the error of the input var- 
iables was used to determine the sources of the tonnage and 


recovery errors. 


Two types of tonnage were calculated: in-place coal and 
reserves. The in-place coal is the amount of coal that was 
in a bench prior to stripping and was estimated using 
thicknesses measured on the highwall after removal of the 
coal. It was used, along with the coal-loaded-out to 
determine the recovery. Reserves, on the other hand, were 


estimated for several locations after they were mined out, 
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but using only pre-production data, either drill hole or 
nearby highwall measurements. These two numbers were 
compared to the 'true' reserve figure, which was estimated 


from post-production measurements of the highwall. 


The error of the reserves for locations not yet mined 
depends heavily on the error of the area (VAR[A]) to be 
exposed. Two cases were considered: 

1. VAR[A] = 0, 

22" VARTAI~>"0" 

The first case is when the area iS given, such as a lease 
limit, and corresponds more closely to the strict meaning of 
geological reserves. In the second case, the area is deter- 
mined by the amount of ground that can be exposed in a given 


UNrteOn time, and so 1s haghly variable. 


1.3 PREVIOUS WORK 
Little previous geostatistical work has been done that 
applies directly to this study. “Most? or the previous work 


is geological. 


1.3.1 Previous Geostatistical Work 

Techniques of geostatistics have only in the last five 
to ten years slowly made their way into North America after 
twenty to thirty years development in France and South 
Africa. Journel and Huijbregts (1978) and David (1977) have 


the most extensive bibliographies of case studies. These 
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two, along with Clark (1979a), give the best practical 
descriptions of the methodology in English. Matheron (1971) 


describes the theory behind the applications. 


Over 200 deposits have been analysed by geostatistics, 
(David, 1977, p.70), but there are few previous studies sim- 
ilar to this one. Journel and Huijbregts (1978, p.595) list 
36 deposits that they used as examples of geostatistical 
analysis. Of these, half the studies are on tabular 
deposits of three types: bedded, residual soil and veins. 
Half of these tabular deposits are bedded, but only use 
thickness of the strata in conjunction with the accumula- 
tiont, never alone. None deals with the geostatistics of a 


coal mine. 


What work has been done on coal mines has for the most 
part been done on coal quality. The only person to have 
done a geostatistical study of plains coal was Flint (1978). 
He is referred to often in this study, so when no year is 
SuffLixed pwitushould.be.taken;to refer to his M.Sc. thesis, 
Flint (1978). He estimated the resources and their 
variances for three coal zones (including the Estevan Zone) 
in the Ravenscrag Formation, two of which are present over 


the entire Estevan Field. 


+ Accumulation is the average grade of an ore intersection 
times the length of the intersection. 
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Flint found the net coal thickness to be stationary and 
isotropic. The variogram model was spherical with ranges of 
200 m, 800 m to 1000 m and 4000 m to 6000 m; and a fairly 
large nugget effect. The coefficient of variation of the 
tonnnage was about 20%, with over two thirds of the error 
derived from error in the determination of the area 
underlain by coal. The error in the density contributed a 


minor amount. 


Sabourin (1975) examined the sulphur content in the 
Lingan Coal Mine of Cape Breton. Rendu and David (1979) 
studied the calorific value of the Hat Creek coal deposit of 
B.C. They introduced a new method of uSing geostatistics 
for studying thick sedimentary deposits by treating the 
deposit more like an irregular ore body than a regular 
Strata-bound coal deposit. Irvine et a]. (1978) have 
described a detailed coal resource inventory of the 
Ravenscrag Formation, using computerized traditional methods 
involving SURFACE IIt and arbitrary distance-proportional 


weights, rather than geostatistics. 


1.3.2 Previous Geological Work 

Previous geological work in the area is extensive, 
beginning with preliminary studies done for the Palliser 
expedition in 1857. G.M. Dawson did the first serious geo- 
logical investigations in 1875, but most of the geological 


+ SURFACE II is a grid manipulation and contouring package 
(Sampson, 1978). 
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work was done in the first third of this century on 
outcrops. More work has been done in the late sixties and 
early seventies useng drill hole data. ~Ah account of previ- 
CuGSeWOrKi 1s deravled in rlint’ (1978, pti). THe most 
extensive bibliographies of the previous work are in Irvine 


et a]. (1978) andflint (1978). 


1.4 REGIONAL, LOCAL AND MINE GEOLOGY 

The Boundary Dam Mine is located in the plains 
physiographic region of central North America. The 
topography 1s mostly flat with low rolling hills with deep 
and sometimes broad valleys. Geologically, the study area 


1S near the northern rim.of the Willystion Basin. 


1.4.1 Regional Geology 
The following discussion applies to the part of 
southern Saskatchewan underlain by the Ravenscrag Formation 


Wma ur Gut. 2s. 


1.4.1.1 Regional stratigraphy 

Paleozoic deposits, approximately 2000 m thick, are 
dominantly shelf carbonates, evaporites and basinal shales 
of Cambrian to Mississipian age (Figure 1.3), deposited over 


Churchillian basement. 


A major unconformity separates Mississipian strata from 


the Jurassic beds. Several transgressions and regressions 
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Figure 1.2 Regional geology 


i of -STep se 
Fan 
eS 23 oT rs es yee 


oor 


eon i = 
asitijemolia Ry 


devort aoitjuloa ifse -j{—- 


(A0@S! qaM ,OT2T) asipucd rosie Seidibe 
| | Tae, ae 
_ gpoloep LsacipeA 


a6 


12 


EMDPLESS Gp... os 
RAVENSCRAG FM. 


Frenchman £m, --— 
Fastend Fm. / 


PPPS PLE LILI AL ALA ALLL LIE is 


fo 
Cretaceous 
ce 


Bearpaw Fm. & 
Riding Mountain Fm. 


Mesozoic 


S 
i) 
&) 
sh & ss Jurassic 
Yn 
cc) 
ku 
4 
ie) 
i= . : . . 
Mississippian 
& 


dol & sh 
Devonian 
dol & sh , 
| siturian 
| ordovician 


sltst & sh 
Cambrian 


Figure 1.3 Regional stratigraphy 


Elk Poants Gpy 


Data from Douglas 
(O70 Chart. 1 Pi) 


Paleozoic . 


WSEPOSCTC 


w + Se ws = @ 


¢ 
‘ 
: 
4 
pone sciene i tliaietineettdainiae 


fi 


z ae ta we Natori 
lide eee eas 
5 
2 
“9 . 
; 7 


Pe LAR AIELLO GOES EM. 


| 
a ; 
| qe ioeved 
t 
! 
i 
| 


TH BOTs, ERE 


re eee 
; fe 2 Tom. Dede 
§ a ee eee i eh tals Mtoe wh, ° 


; : : 


esloued m a 7: ‘ 
(123. 3% over 


qiqsigiteite isnoteea 


As 


during the Mesozoic resulted the deposition of approximately 
1500 m of clastics, both marine and non-marine, shed largely 


from the Cordilleran orogen to the west. 


Cenozoic deposits are largely non-marine clastics 
containing several unconformities and one thin marine unit, 
the Cannonball Formation, which occurrs mainly in the 
adjacent United States. The Ravenscrag Formation of 
earliest Tertiary age (Danian and possibly a little later) 
was deposited over the central portion of the Williston 
Basin. The Ravenscrag coal measures are found in four 
coalfields (Figure 1.4): 

t.. Cypress, 

2. Wood Mountain, 

See wid VOW enune lh, 

4, Estevan. 

Occasionally patchy Eocene to Pliocene deposits are found 
above the Ravenscrag Formation. A thin veneer of Pleisto- 


cene glacial deposits caps the sequence. 


1.4.1.2 Regional structure 

The strata form a broad syncline, which plunges gently 
at approximately 3 m/km to the south into the Williston 
Basin. Locally this pattern is modified by solution of 
evaporites in the Prairie Formation (Middle Devonian) and 
collapse of the overlying strata to form gentle synclines 


(Peaceineml.2 peat). 
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Figure 1.4 Coalfields in the Ravenscrag Formation 
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1.4.2 Local Geology 
This description below applies to the Estevan Coal 


Field (Figure 1.5). 


1.4.2.1 Local stratigraphy 
Marine shales of the Riding Mountain and Bearpaw 

Formations (Figure 1.6) are overlain by the marginal marine 
Eastend Formation which in turn is overlain by the latest 
Cretaceous non-marine Frenchman Formation. The Cretaceous- 
Tertiary boundary is usually placed at the contact between 
the Frenchman and Ravenscrag Formations, though there is 
some disagreement. Total thickness of Cretaceous strata in 


the area is about 1000 m. 


The Ravenscrag Formation may be as much as 550 m thick. 
The lower contact of the Ravenscrag Formation is conformable 
and usually defined as the base of the first coaly material 
above the Bearpaw or Riding Mountain shale. The upper 


contact is everywhere erosional. 


Equivalents to the Ravenscrag Formation are parts of the 
Willow Creek and Porcupine Hills (Paskapoo) Formations in 
Alberta, and the Fort Union Group in the United States. 


W 


The lithology of the Ravenscrag Formation is "**-grey 
and buff-weathering shales, fine sandstones, siltstones and 


sandy shales--+-" (Flint, 1978, p.16), usually poorly 
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Figure 1.5 Local geology 
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consolidated, commonly containing lignite seams in up to 19 


zones. 


The depositional environment has been interpreted as an 
alluvial plain covered by forests, swamps and marshes ina 


humid temperate climate (Flint, 1978, p.17). 


Irvine et a]. (1978, p.49) have divided the coal zones 
in the Estevan Coal Field into two groups (Figure 1.7): 
fourteen zones of deep coal, too deep to mine at present; 
and five zones of shallow coal, which can be, or are being 
mined at present. The shallow coal zones (Figure 1.7, p.19 
and Figure 1.8), in stratigraphic order from the top down 
are: 
Wepesnort Creek, 
2. Roch Percée, 
a “Souris; 
4, Estevan, 
5.) Boundary. 
The Estevan and Boundary Zones, are the thickest, most 


widespread and most heavily mined. 


The Ravenscrag Formation in the study area is 
erratically overlain by unconsolidated sands and gravels of 
the Plio-Pleistocene Empress Group. The whole area is 
covered by Pleistocene till except where it has been removed 


by fluvial erosion. 
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Figure 1.8 Shallow coal subcrop in the Estevan area 
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1.4.2.2 Local structure 

The Ravenscrag Formation dips gently towards the 
southeast, but superimposed on it are a few very shallow 
southeast-plunging undulations (Figure 1.9). Solution of 
the Prairie Evaporite Formation has produced local collapse 


features. 


In places the poorly-consolidated Ravenscrag Formation 
has been structurally disturbed by thrusting from a Pleisto- 
cene ice sheet. This disturbance takes the form of either 
disorganized strata or imbricate thrust sheets, usually 
within 25 m of the bedrock surface, though it can go deeper 
(Flint, 1978, p.153). Most ice-push features occur beyond 


the study area to the northeast. 


1.4.3 Mine Geology 

The Boundary Dam Mine is located at 103° 00' West 
Longitude, 49° 05' North Latitude or Range 8 on the boundary 
between Townships 1 and 2, West of the Second Meridian 
(Figure 1.8, p.20). The mine is approximately 3 to 7 km 
south of the town of Estevan and 8 to 12 km north of the 
U.S. border. It is bordered on the north by the Souris 
River Valley and on fhe west by the valley of Long Creek, 
both of which cut through the Estevan Coal Zone. The 
southeast boundary is formed by thinning of the coal towards 
a major northeast-trending, sand-filled paleochannel cut 


into the seam (Figure 1.10). The study area (Figure 1.11) 
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Figure 1.9 Local cross section 
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encompasses the three pits active in the mine from 1977 to 


1978 along with the surrounding drill holes. 


Topography over the mine is flat with about 5 m relief, 
but incised by the Souris River in a broad glacial meltwater 
'channel, and Long Creek in a narrow Holocene ravine, both 


about 30 m deep. 


1.4.3.1 Stratigraphy at the mine 

Overburden, 15 to 20 m thick, 1S composed of glacial 
till, successively underlain by discontinuous deposits of 
sand and gravel of the Empress Group and poorly consolidated 


grey silty clay of the Ravenscrag Formation. 


Seams in the Estevan Zone (Figure 1.12) are named A to H 
from the top down (nomenclature is after Irvine et al., 
1978, p.60). Seams A and B are thin, uSually occurring as 
part of the fill in the major paleochannel. Seam C is a 
rider seam, too thin to mine in the Boundary Dam Mine. 

Seams D, E, F and G are mined, while seam H is a thin, 


discontinuous rider at the base of the zone. 


Maceral content of the coal is 65% to 85% huminite, 10% 
£o 2ovuinertinite and 5%) to» 15% exinite (irvine etal. 
1978, Figure 13, p.46). The average proximate analysis of 
coal from the seam is listed in Table 1.1. Average gross 


thickness is 3.9 m and average net thickness is 3.3 m 
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Figure 1.12 


Stratigraphic column of the Estevan Zone 


in the Boundary Dam Mine 
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Table 1.1 


Moisture 


Ash 


Volatile matter 


Fixed carbon 


Calorific value 


Rank 


Sulphur 


the Estevan Coal Field 


eS Thaey! 


24.6% 


28.8% 


15 100 kd/ke 


Lignite 


0.3% 


This data is from Adamcewicz (1963, p.101). 
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Typical proximate analysis of the coal from 
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Figure 1.13 
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(Figure 1.13). The descriptions below refer primarily to 
the portion of the mine studied in this project. Two main 
partingst of grey and brownish clay are recognized by mine 
Surveyors in the study area (Figure 1.14). The upper, prob- 
ably D-E, parting is present in the southern part of the 
study area, wedges out towards the north (Figure 1.15 and 
Figure 1.16) andiis 0.3m thick and 0.9 m fyom the top of 
the mined coal. The lower parting, probably EF, occurs in 
the north, wedging out towards the south and is 0.4 m thick 
and 2.1 m from the top of the mined coal. The F°G parting 
may have been measured ina few places in the south end of 
Dites ert so fit 1s) OJlym thicksand 0.6m Brom therbase of 


the seam. 


The pattern of splitst and partings seen by geologists 
in the drill holes from the study area is not quite the same 
as that seen by the mine surveyors. Different geologists 
seem to have their own criteria for identifying a parting 
and consequently, holes with one number of splits can be 
Surrounded by infill holes with a different number of 
splits. The floor consists of poorly-consolidated grey 


clay. 


A few linear areas (Figure 1.17) occur in the pit data 
where the parting and/or coal is thicker or thinner than it 
t A parting is a layer of rock within a coal seam. 


+ A split is a subdivision of a seam that is separated by a 
parting from another split. 
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Figure 1.14 Stratigraphic column for the study area 
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should be. These areas are probably minor paleochannels. 


1.4.3.2 Structure at the mine 

The coal zone dips gently to the southeast with minor 
undulations (Figure 1.18). No salt solution tectonics have 
been found, but there has been some collapse into a few of 
the old underground mines, particularly near the Souris 


Valley, outside and north of the study area. 


1.5 COAL MINING IN SOUTHERN SASKATCHEWAN 

Mining from small underground workings in the walls of 
the major river valleys began in the early 1870's. The 
first open pit mine started in 1927, and consolidation of 
the many small underground mines occurred about the time of 
the Second World War, with the conversion to Strip mining 
complete by 1955. Six mines (Figure 1.19) operate in the 
area at present. Seventy percent of the coal mined in 
Saskatchewan is used for power generation and most of the 
rest iS used in pulp and paper mills and for briquetting, 
(Adamcewicz, 1982, p.105). Coal from the Boundary Dam 
Minet, with a nominal annual capacity of 1.8x10° tonnes, is 
used exclusively by the Saskatchewan Power Corp. in the 
mine-mouth Boundary Dam Power Station. A detailed history 
of coal mining in Saskatchewan may be found in Irvine et al. 


(LOg6 men tas). 


+ The mine is owned by the Manitoba and Saskatchewan Coal 
Co., which is owned by Luscar Ltd. 
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Figure 1.18 Cross section through the mine 
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Figure 1.19 Operating coal mines in Saskatchewan 
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The Boundary Dam operation, like all those in 
Saskatchewan, uses Strip mining techniques. Walking 
draglines, working from a prepared surface, remove the 
overburden in one pass and cast the spoil onto the adjacent 
bench from which the coal has been removed. Scrapers and 
graders clean off the surface of the seam and remove any 
partings thicker than about 1/3 m. Partings less than this 
are not selectively mined. Finally, the coal is loaded into 
trucks by electric shovels for transportation directly to 


the power station. 
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2. DATA ACQUISITION AND PRELIMINARY ANALYSIS 


2.1 DATA SOURCE 
Data for this study was acquired from the Boundary Dam 

Mine of the Manitoba and Saskatchewan Coal Company. Most of 

the data consists of seams and parting thicknesses, but also 

included are some mine records of coal density and coal 
actually loaded out of a pit during a particular time 
interval. The thicknesses have two sources: 

1. Daily measurements of the mined seams in the pit 
highwalls by the mine surveyors (Figure 2.1). These is 
nominally 30 m apart along each bench, which are approx- 
imately 30 m wide. Data was collected from the three 
pits (1, 2A and 3) that were active in 1977 and 1978. 

2.) Poroducei om *drildang on 100, 150 and 200 m grids 
(Figure 2.2), covering three sections that surround the 
DIts, 

In both cases the original blue-line drawings of the mine 

plans with data location and thicknesses were xeroxed 

because the original drawings were too large to fit on the 

digitizing tablet, a Summagraphics Bit Pad. The X and Y 

coordinates of each data location were digitized and the 

thicknesses entered simultaneously on the keyboard of the 
attached HP2648a graphics terminal. All thicknesses 
pertaining to the coal were entered, both partings and 
splits. The pit measurements either had no partings or one 


parting, rarely aS many as two; whereas the drill holes had 


38 


.yosqme) ses Ae 
sgeogndoing paigvag ah. 


ttaneb isos 36 abieze2 


ape 8 ent mb tig.s “ = 


$5 wee ons. oven zsareendod 


7) 
bt 


bssselléo ecw S300 obiv n oe 


5 
e 
‘4 


sidgexesame2 6 , 3eidas Jue 


us 


sodjie 2inemewweasem tiq sat 


i aa 
+ 


9a3_ %0 Binen 
wy we “eq in ade we 


4 


2 ere fe 294q8 m-O8 | 


siaw jad? (& Soe aS ai 
OOF de oni tiie neds 
is QNLSEVOD), is. 
suid Lea terse “as asass Adee 
2¢42 Sas notzacel aish dgi i. 
4 @patwet® Pea tpi ae sf? —< 
igesol sisb doss to zotsnibt 
fauosnasiumie Saeagdne aoe i: 
cimist estdgsiup soeasen t 


sine s1ew [soo sd%9 of pniak 


cows as yvsm es qfeTst Vom 


She) 


e®veet Fe ese ete eeoe FF eevee eg vee ® 
eo eeeeeeeree® Feat Oe gga ce®e@ © 


¢©e®%®&® @e 8 @ @ 
ecoere® © © FG @ @ 
oe ® © @ &© ® 6 


metres 


Location of the pit data 


Figure 2.1 


[> 


eo3 tom. 


~ 


Figure 2. 


metres 


Location of the drill 


hole data 


40 


41 


up to five partings. A program that calculated map 
coordinates from the pad coordinates also calculated net and 
gross coal thickness. See Appendix 11 for an example of the 
resulting data set. Outlines for the pits were also 


digitized at points no further than 30 m apart. 


2.2 CLASSICAL STATISTICAL ANALYSIS 

Classical statistics is non-geostatistical statistics 
or what is usually referred to as simply ‘statistics’ and 
which does not deal with any regionalized variables.ft Any 
Statistical analysis must uSe data that is as free as possi- 
ble from all mistakes, so the data was cleaned, as outlined 


in Appendix 2.1. 


Before any geostatistical analysis can be done, the data 
must be explored with classical statistics. Not only is it 
good practice, but the geostatistics used in this study make 
two assumptions that can be checked only with classical 
Statistics. These assumptions are that the data is from one 
population that is normally distributed. In addition, the 
values of some classical statistical parameters affect some 
of those in geostatistics. The variogram sill is 
theoretically equal to the population variance, and the 
spatial distribution of the data determines how unbiased the 
arithmetic average is for a reserve estimate. 

+ Briefly, a regionalized variable is a variable whose 


statisticstare’affected by:its location... iSee Appendix 1.1.1 
for a more complete definition. 
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2.2.1 Statistics of the Net Thickness 

Tonnages are calculated from the product of the 
density, area and the average net thickness of the mined 
seam. The last variable can be calculated by taking the net 
thickness before or after averaging. The latter is the 
better way of determining the average net thickness, as 
outlined in Appendix 2.2. To do this the net thickness is 
required at each samplet point. Three classical statistical 
properties of the net thickness were examined, its distribu- 


tion, mean and variance. 


2.2.1.1 Distributions of the thicknesses 

If the thickness distributions are lognormal, then 
disjunctive kriging with de Wijsian variogram models should 
be used. If the population is normal, then more common 
variograms and kriging can be used. All the net thickness 
distributions are approximately normal except for some 
mixing of populations. Two methods were used to analyse the 


distributions: STPK and normal probability graphs. 


STPK is an interactive statistics package that can com- 
pare the empirical quantiles of the data to the best-fit, 
theoretical quantiles for any of the following cumulative 
distribution functions: normal, gamma, chi-squared, F, 
Student's t, beta, Cauchy, Weibull, half normal, uniform, 

+ 'Sample' is used in this study to refer to any measurement 


of the: seam thickness, either in a drill hole or on the 
highwall. 
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binomial, Poisson, geometric, and hyper-geometric. 


The fit of the data was compared visually with each of 
the above distributions. The only good fits were the normal 
and gamma distributions (Figure 2.3). Only one pit is 
illustrated due to data limitations in STPK, but the other 
pits show a similar pattern. The gamma function fits a 
little better because it is a three parameter function and 
SO can be skewed. However, the differénce is negligible, so 
the distribution was taken to be normal. The next best fit 
was the Cauchy function (Figure 2.4); all the rest were much 


worse. 


A program was written that plots the data ona 
Cumulative normal probability graph both to visually esti- 
mate the normality of the data, and to find any mixing of 
populations. Data from one normal population will plot ona 
straight line with its location determined by the mean and 
its slope by the variance. Multiple normal populations will 
plot on straight line segments. All the pit data plots on 
reasonably straight lines except for tails toward the thin 
values (Figure 2.5). This indicates a slight to moderate 
skewness to the left (negative). The points in the tail are 
from discrete areas in the seam that are distinctly thinner 
than usual (Figure 1.16, p.32), areally associated and often 
linear (Figure 1.17, p.33), suggesting the effect of small 


paleochannels. 
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The overall distribution of the data may be three param- 
eter lognormal. However, an adequate interpretation is that 
there 1S a mixture of two populations: 

1. most of the area, normally distributed, 

2. channelled areas, probably normally distributed. 

The pit data was considered to be one population since the 
departure from one population is small and the channelled 
areas could not be subdivided out and handled separately. 
These areaS are too small to be treated on their own, as 
they contain less than the 30 to 50 sample minimum for 


calculating a variogram. 


The drill data also plots as straight-line segments 
(Figure 2.6) but shows a much larger tail than the pit data 
(Figure 2.7). This longer tail is probably due to the 
larger sampled area having portions of the seam too severely 
channelled and thin to mine. The two populations in the 
Grill hole data were too intimately mixed over the study 


area (Figure 2.8) to separate into simple areas. 


2.2.1.2 Mean and variance of the thickness 

The mean of the net thickness of all the pit data is 
Seo Meandeot thegdril lsholerdatawis, 2.5m (Table ™2.1') 2 eThe 
Grill hole average is less because some holes were drilled 
where the coal was too thin to mine, generally in the west 
of the study area. In addition the drill geologists on some 


holes may separate out partings that will not be selectively 
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Figure 2.6 Cumulative normal probability curve for the 


drill hole data 


ee Repewrceemmeres ee ee MN Seem 


-< 


votive (igen bssim” 


a See 


olisivgeg bef fennan> _ 


svisz yiilidsdetq Ismion eviteigmey? 


ats stod ff9r3b 


slelee 


Se) 


(%) 
No) 
oO 


on 
jo) 


Cumulative Normal Probability, 
= 


Figure 2.7 


Wed 
a) 
DRGee HOLES 967 ey 
Wa = 


2.0 350 
Net Thickness, (m) 


Comparison of the pit and drill hole 


cumulative normal probability curves 


i 
es 


“ 
> 


eS ere 


° . 
ee” ee ed 


aoziasqmod 


evitsiums 


50 


1 “SPLIT e 
27 SPLITS 2 
3 SPLITS 90 


1) e 
sd e 
ee @ e 
eee @¢ @ 
eee @ @ 
e@©e¢ee @ 
eo e- 86 e 
eee¢ee 

Oo 


0 500 


metres 


Figure 2.8 Location of the partings in the drill hole 
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Table 2.1 

Area nr 
Pit 1 272 
Pree ca 201 
Prt 3 158 


All pits 631 


Drill holes 336 


= 


Summary statistics of the thickness for each 


pit and the drill holes 


Mean Variance Standard 
thickness Cm) deviation 
(m) (m) 
Siete) 0.084 O29 
See..6 Oe 0siG 0.19 
Shaler 0.063 Oig235 
Byler} OF073 Oeze 
3.29 ORL U2 55 


+ n = the number of samples. 
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In geostatistics, ‘population variance’ is the same as 
"variance' in classical statistics. The overall variance of 
the pit thickness data is 0.073 m* and of the drill hole 
data is 0.278 m*. The drill hole variance is larger for two 
reasons. Firstly, the geologists would be less likely to 
compare the meaSurements for consistency between holes. 
Whereas, the Surveyors may make the pit thicknesses more 
consistent and less variable since they are usually done 
sequentially, in both time and space, on a continuous expo- 
sure. Secondly, there are more splits in the drill data (up 
to four partings) than in the pit data (rarely more than one 
parting). The extra partings will remove a more variable 


amount from the the thickness of the drill hole data. 


2.2.2 Location of the Samples 

The spatial distribution of the data was determined by 
visual inspection of the data postingt, which shows that the 
pit samples are on a semi-regular grid at approximately 30 m 
Spacing with a few large areas of missing nodes (Figure 2.1, 
p.39). Postings of the drill data show that the holes are 
on three regular grids with several small areas of missing 
holes (Figure 2.2, p.40). The basic pattern is a 200 m 
square grid that has been infilled to a 100 m square grid in 


the northeast and a 150 m square grid in the southeast. 


+ Posting was done with SURFACE II (Sampson, 1978). 
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3. STRUCTURAL ANALYSIS OF THE DATA 


The term ‘structural’ is used in this study in the geo- 
Statistical sense: it has nothing to do with faults or 
folds. The statistical structure of a regionalized variable 
is its spatial variability, namely its drift and particu- 
larly its variogram. The theory of structural analysis is 


outlined in Appendix 1. 


3.1 VARIABLES USED IN THIS STUDY 
The tonnages in this study are calculated by the 


product of three variables: thickness, density and area. 


Thickness satisfies the conditions in Appendix 1.1.1 for 
a variable to be analysed by geostatistics and is the main 
geostatistical variable in this study. Thickness is a 
continuous variable distributed in two dimensional space and 


any linear combination of thicknesses is also a thickness. 


Density also satisfies the same conditions: it is 
continuous, distributed in three dimensional space and any 
linear combination of densities is also a density. However, 
only three densities were determined and at least 30 to 50 
measurements are required for geostatistical applications. 
Thus, even though density is regionalized and should be 
analysed by geostatistics, classical statistics had to be 


used. 
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Area is not a regionalized variable. It is not defined 
at a point and so does not meet the necessary conditions in 
Appendix 1.1.1. As a result it was analysed using classical 


statistics. 


3.2 DRIFT ANALYSIS 

Ordinary kriging, as used in this study, assumes that 
there 1S no drift. This assumption was tested by plotting 
on Figure 3.1 the one dimensionalt drift values estimated by 


the variogram program. 


The drift curves do not rise continually and have low 
maxima of 0.12 m (3% of the average thickness) for the pit 
data and 0.14 m (4%) for the drill hole data. Furthermore, 
the variograms (Figures 5.1 and 5.2 in Appendix 5.1) are 
bound and so show no drift effect. Both sets of data show 
little or no drift, so ordinary rather than universal 
kriging will suffice. Any small drift that is present 
should not affect the reserves eee Since David. (1974, 
p.114) states that the estimation is very robust with 


respect to a drift. 


3.3 VARIOGRAM ANALYSIS 

Analysing the variogram is by far the larger and more 
important part of structural analysis, if there is no drift. 
Variograms should only be estimated over a homogeneous 


+ A description of how a one dimensional function can be 
estimated is in Appendix 1.3.3.2. 
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region, which is an area of uniform geological, statistical 

and geostatistical properties. The pit data has two 

homogeneous regions; a north region consisting of pits 1 and 

3, ano @ South region formed from pit 2A, while: the drill 

holes had only one region. These regions are determined in 

Appendix 3. Usually variogram analysis is done in three 

Stages: 

1. calculation of the raw variogram, 

2. regularization or correction of the raw variogram to 
form the experimental variogram, 


3. modelling of the experimental variogram. 


3.3.1 Calculating the Raw and Experimental Variogram 

The raw variograms were calculated with the parameters 
discussed in Appendix 5.1 and drawn using the conventions 
listed in Appendix 1.3.1. Regularization is too small to be 
corrected for (Appendix 5.2), so the raw variograms are also 
the experimental variograms (Figure 3.2 and Figure 3.3), 
which were used for modelling. After modelling, the 
resulting variogram models were used as an input to the 


kriging program. 


3.3.2 Modeling the Fundamental Variogram 

Experimental point variograms are too ragged to use as 
input to the kriging program: a mathematical expression 
rather than the raw points is needed. Modelling is the pro- 


cess of fitting a mathematical curve to the calculated 
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points. It is largely an art because there are no 
algorithms to make a best fit of an ideal curve to the 
experimental variogram. Fortunately, "*-+*a visual fit is 
USval ly SUETIC lente.) Maccording (to David) (1977), to. 122) 
Parameters are chosen and a curve is calculated and compared 
to the experimental variogram. Then, the parameters are 
adjusted and the process repeated until a reasonable fit is 
achieved. Flint (1978) wrote a program to do this 
interactively; itS operation is briefly outlined in Appendix 
6. Some semi-automatic programs are on the market, but none 


is available on the University of Alberta computer. 


3.3.2.1 Variogram models 

Preliminary inspection of the experimental variograms, 
forsboth the pit <—(kigure 3.2; 0.57) and drill hole data 
(Figure 3.3, p.58), shows that the models are from the 
transition class. For the reasons listed in Appendix 1.3.6 
the spherical model was chosen for this study. The programs 
for calculating the variograms, modelling them and for 
kriging are based on the spherical model. Extra programming 
effort would have been needed to use any of the other two 
possible models mentioned below, with little prospect of 


increased accuracy. 


Two nested, isotropic spherical models (Figure 3.4 and 
Figure 3.5), with no nugget effect are used for the modelled 


variogram. Their parameters are listed on Table 3.1. The 
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Table 3.1 Variogram model parameters 

Parameter Souci p1it NOrtn pits All pits Drill holes 
Nugget 0.0 0.0 mae O20 

effect, Co 


Short range structure 
Stillenc; 03.02 amé 0.045 m? 02,035) m4 Oe, m- 


Range, a, S50 mz 50 mt 5 Ona 150 4m 


Longa nange, Structure 


Sil eCe 0. 01,5854 OnerO 3.5y sm" 0x.10 35. mez Ooi? 2s st 
Range, a2 180 m 450 m 500 m 1500 m 
Overall 0.0358 m? 0.080 m?* 0 ,0:7,0) an? Oyan3 Sagan ¢ 
Siu 7 Cy 
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is good within the limitations of geostatistical 


modelling, but it is possible that a combination of a hole 


effect model plus a short range spherical model would also 


Pact, 


not 


Saas 


A hole effect model (Figure 1.9 in Appendix 1.3.6) was 
used for the following reasons. 

The spherical model fits well over the distances used in 
the kriging calculations: a few hundreds of metres for 
the pits and about 1000 m for the drill hole data. 

The hole effect model would only improve the fit beyond 
these distances. 

Journel and Huijbregts (1978, p.171) state that "For 
estimation purposes, an experimental hole effect that is 
open to doubtful interpretation, or not very marked, can 
Simply be ignored::--" 

The modelling and kriging programs would have to be 


extensively modified to handle it. 


2.2 Nugget effect 


None of the models has a definite nugget effect. The 


intersection of the line through the first two points with 


the vertical axis was used to estimate C,. The north pits 


(Figure 3.4, p.60) and all the pits together (Figure 3.6) 


may have a small nugget effect of about 0.01 m*. The south 


pit and the drill hole variograms (Figure 3.2, p.57 and 


Figure 3.3, p.58) intersect the horizontal axis, which 
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suggests parabolic behaviour near the origin. 


A nugget effect can have two causes: human nugget effect 


and a structural nugget effect. 


David (1977, p.99) defines a nugget effect as "*++the 
uncertainty on the value of the sample itself." A human 
nugget effect, though small, certainly exists in this mine. 
All readings, pit and drill hole, were nominally taken to 
the nearest 0.1 ft (0.03 m), which has an error of +0.015 m. 


The mean square error is then; 
COn015)m) > =) 02000225 m- 


This iS an estimate of the human nugget effect (Gandin, 
19685 ep. 29) ou.SincestOS 12 ft sMthesbest that the -mine 


Surveyors can hope to do it is a lower limit. 


In practice, the Surveyors and mine geologists were not 
always able to meaSure to the nearest 0.1 ft, but sometimes 
had@t6. y6und to the néaresto0os 27 £tPyedots: fttor si *ftowhen*®the 
precise position of the contacts was uncertain. If all 
readings were accurate to +1 ft, then the human nugget 


effect would be 0.0225 m*, which is an upper limit. 


The true human nugget effect will be somewhere between 
0.0002 m* and 0.02 m*. The geometric mean of these two 


numbers, approximately 0.001 m?, is a reasonable order of 
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magnitude estimate of the human nugget effect. A more 
precise value of 0.0017 m* was calculated from estimates of 
the number of each type of approximation (+0.2 ft, +0.5 ft 
and +1 ft) that the surveyors made. This value is too small 
to model because the noise in the experimental variogram is 


larger by at least an order of magnitude. 


A structural nugget effect probably does not exist. The 
Grill hole data may be parabolic near the origin, which 
implies no nugget effect. The pit variograms also probably 
have no structural nugget effect for two reasons. Firstly, 
both pit variograms seem to have the same type of visible 
structure, both short and long range, even if the parameters 
vary. Therefore, it is likely that any invisible (nugget) - 
Structure would also be the same for both variograms. 
However, one pit variogram is parabolic while the other has 
an apparent nugget effect. This difference is probably 
random fluctuation with both true variograms having a simple 
linear behaviour near the origin. Secondly, any nugget 
effect would be due to a short range structure with a range 
less than 30 m? and a sill of approximately 0.02 m*. This 
means the thickness would vary up to +0.3 m (20+) within a 
distance of 30 m. This sort of variation cannot be seen in 


the continuous exposure along the highwall. 


+t o is used here and elsewhere in this study as an 
abreviation for the standard deviation. 
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3.3.3 Structural Interpretation of the Model 
Structural interpretation is the interpretation of the 


variogram model parameters in geological terms. 


3.3.3.1 Interpretation of the model parameters 

Transition models usually apply to phenomena that are 
fairly regular and 'well behaved', as sedimentary 
thicknesses often are. The spatial structure includes a 
zone of influence, inside which the data is better correla- 
ted the closer they are to each other and beyond which the 
Gata is independent. Such data has a finite population 
variance, which is usually the case with thickness vario- 


grams. 


The variograms show no anisotropies other than the 
Slight effect due to measurement techniques in the pits 
discussed in Appendix 5.3. Most geological anisotropies in 
a tectonically undisturbed area will align with the 
preferred paleo-environmental directions; in this case, the 
major paleochannel (Figure 1.10, p.23). Neither the pit 
(Figure 3.7) nor the drill hole variograms (Figure 3.8) 
change depending on the direction with respect to the major 
paleochannel. This is probably because there is no strongly 
preferred direction of deposition in an essentially flat 
Swamp. However, some anisotropy could be present nearer the 


paleochannel, or in other variables such as ash content. 
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As expected, the overall sill is well defined and gener- 
ally equal to the population variance. The long range 
Structure is probably due to the depositional history of the 
zone. When two drill holes are greater than 1500 m apart, 
the thicknesses are poorly correlated because they will tend 


to be in areas influenced by different partings. 


The overall sill is larger for the drill hole model than 
for the pit models because of the differences in the 
population variance listed on page 52. The structures of 
the pit data are probably different because of different 
depositional micro-environments associated with the two 
partings. The sill is higher on the north variogram because 
of all the small paleochannels, which will increase the 


population variance. 


3.3.3.2 Interpretation of the structures 

The short range models may be due to the method of 
measurement or deposition. In all three cases, the short 
range iS approximately equal to the sample spacing 
(Table 3.2). The differences between the range and the 
Sample spacing are because the ranges are only a rough esti- 
mate using one to two points from the experimental vario- 
gram, while a good eStimate requires at least three to four 
points. On the other hand, the short range sill (C,) is 
always approximately equal to the long range sill (C2) 


(Table 3.1, p.62), which suggests that the short range 
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Table 3.2 


Area 


North pits 


South pit 


Dribkl’ holes 


1s 


Comparison of the short ranges of the vario- 


gram models with the sample spacings 


Short ranget Nominal 
(m) sample spacing 
(m) 
50+ 30 
50+ 30 
150 150+50 


foerom Table 3.7), “p62. 
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Structures are also due to deposition rather than measure- 


ment effects. 


The short range structure of the drill hole model is 
approximately the same as the long range structures of the 
pit models, taking into account the very differerent sources 
of the data. The discrepancy may be due to an error in the 
drill hole short range, which was estimated from only one 
point. The sills are also moderately close: 0.17 m? and 


Op35%m-. 


Nugget effects (C,.) are not apparent on any of the 
variograms because the meaSurement errors are too small, 
about, 0.00 tim-> and) Phere 1S ino Vvisiblevstructure: with “a 


range less than the sample spacing in the pits. 


3.3.4 Comparison of the Models From This Study With Flint's 
models 
Flint (1978, pp.107 and 108) modelled the variogram for 
the net thickness of coal over the full extent of the 
Estevan Zone. His model was isotropic, spherical and nested 


with the parameters listed in Table 3.3. 


It should be possible to correlate the overlapping 
structures for two reasons. Firstly, the same variable was 
Studied: the net thickness of the Estevan Coal Zone as 


measured in drill holes. Secondly, the magnitude of sample 
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Table 3.3 Comparison of the variogram model parameters 
from Flint (1978) with the drill hole model 
from this study 

Parameter Flint's modelf This studyt 

Short range structure 

Nugget effect, Co O09 2m. 

Stel aC, 0.17) me 

Range, a, <2 50h. ty Loon 

Intermediate structure 
Sipe ll ed Gr Ones m7 O222am- 
Range, a2 1000 m 1S 007m 


Long range, structure 


uckcons (Flaite r Gd 7:8:,v ep -rlOvas. 

Pee rOlmlaO Le s.ul4 pi. OZ. 

tt This is the approximate minimum sample spacing and 
therefore, the maximum range for any short range structures. 
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separation overlaps: 250 m to 34 km in Flint's study and 

30 m to 3600 m in this study. The only structures that 
overlap are Flint's short range structure and the long range 
Structure from this study. Both the ranges (1000 m and 

1500 m) and sills (0.37 m? and 0.22 m?) are comparable. 
Flint's nugget effect corresponds to the short range struc- 
ture of this study. Discrepancies in the correlation could 
be caused by the ragged nature of experimental variograms, 
particularly Flint's, which means that there will be some 


scope for estimating the parameters. 


As a result of the comparison, an overall model for the 
Estevan Zone can be assembled from the models listed in 
Tablet3a3, pervs Pand®*Tablets¥ir*pl6z. “Mode structures 
chosen as closest to the real structures were always from 
the more closely spaced data when there is a choice, since 
they are better estimates of the true structure. The best 
structures are the following: 

1. short range model of the pit data, 

2. long range model of all pit data grouped together, 

3. long range model of the drill holes from this study, 

4, Flin tits long range model. ~ 

These structures have ranges of 50 m, 500 m, 1500 m and 
4000 m. The last three are in an approximate geometric 
progression, each about one third the next larger. The 
exception is the first one (50 m) which is 10 times smaller 


than 500 m. This discrepancy may be due to the smallest 
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Structure being a meaSurement rather than depositional 
effect, as outlined on page 70 above. It could also be due 
to the pit data having a different progression than the 
Grill hole data; in which case, the 500 m structure is also 


Spurious and there is no geometric progression. 


Serra (1968) studied the variability of iron concentra- 
tion in the Lorraine basin of France from the scale of 
petrographic slides to the whole basin. He also found the 
models to be spherical with ranges in a geometric 
progression. This effect may be common in data that can be 


modelled over a wide variety of scales. 


If all the ranges in the Estevan Zone are in a geometric 
progression, then the next larger one should be at about 
12 km. Flint's experimental variogram (p.108) shows a pos- 
Sible structure at this range, but this is not certain 
because the data presented on the graph end at 14 km, 
leaving only a short distance to determine if it truly is 
another structure. In addition, the variogram is becoming 
unstable with large fluctuation errors, since L = 34 km, so 


that 12 km = L/2.8 
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4. SAMPLE SPACING STRATEGY 


When estimating reserves, the uSual objective is to 
obtain the most precise estimate for the least expense. In 
this study the precision of the eStimate is measured in 
terms of its variance; the lower the estimation variance, 
the better the estimate.f The expense of an estimate is 
usually controlled by the number of samples, particularly 
when drilling, and the number of samples is determined by 
the sample spacing chosen. In general, the estimation 
variance will decrease in a nonlinear fashiont as the sample 
Spacing decreases with the precise relationship depending on 


the variogram alone. 


Since the variogram is known in this study, the effect 
of a Sampling program on the estimation variance can be 
calculated prior to further sampling. For instance, it may 
be possible to reduce the amount of planned drilling with 
only a slight increase in the variance of the reserves esti- 
mate, depending on the form of the function between the 
estimation variance and the sample spacing. This function 
was Getermined by kriging the estimation variance for sev- 
eral sets of data with different sample spacings, then 
plotting the results. The variance versus sample-spacing 
curve could have been determined exactly using auxiliary 
+ This asssumes that the reserve estimate is unbiased. 


¢ Non-linearity is due to the dimishing returns of more 
samples. 
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functions; instead, the numerical approximation by the 
computer program was used to demonstrate how this curve can 
be derived when the requirements of auxiliary functions are 


not met; namely irregular domains and sample locations. 


The estimation variance depends solely on the variogram 
and the sample locations, or sample spacing if the data is 
on a regular grid. The sample values themselves do not 
affect the variance (see Equations {7.12} and {7.14} in 
Appendix 7.2.4), which is why the variance can be determined 
prior ¢oO Griiling, once the varicgram is ‘known. For this 
study, aS in most, the variogram was determined in another 
part of the deposit and the intrinsic hypothesis assumed to 


be true. 


4.1 METHOD 
The variance versus sample-Spacing curve was calculated 

only for the drill holes, not the pit data. Measuring the 
highwall is cheap, so maximizing this sample spacing is not 
as criticar as fOr the drilling. In; addition, using the 
several hundred pit meaSurements required (see Appendix 
8.4.1) to calculate the curve by computer would be very 
expensive, and by auxiliary functions, awkward and very 


tedious. 


The data set used was artificially constructed with the 


sample locations on a perfectly regular square grid. This, 
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rather than an actual drilling pattern, was used so that the 

results would not be affected by the idiosyncracies of a 

particular grid. The data values were arbitrarily set to 

one, Since they had no effect. The theory of kriging and 
the operation of the kriging program are outlined in 

Appendix 7 and Appendix 8 respectively. Following are the 

input parameters for the kriging program. 

1. The domain was rectangular, 300 by 500 m; approximately 
the size and shape of pit 2A or a little over two years 
of production. 

2. The artificial drill hole grid was square, centred on, 
and aligned with the domain. 

3. The variogram model used was the drill hole model in 
Tabhe fal ,. p.62. 

The krige and extension variances for a given sample spacing 

were computed using the modified kriging program with the 

best calculation parameters as determined in Appendix 8.4. 

The calculations were repeated for other sample spacings 

Seenece 50 and 500 m by multiplying the sample coordinates 

in the original artificial data set by a scale factor. The 


artificial data and domain simplified this conversion. 


This curve has been calculated for the thickness 
estimation variance, so that it would not be dependent on a 
particular density or area variance. The overall shape of 
the curve will not change from the one for the tonnage, 


since density (for this study) and area are not regionalized 
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and their variances are not functions of the sample Spacing. 


4.2 RESULTS 

The results of the calculations are plotted on 
Figure 4.1 and Figure 4.2, which show, as expected, the 
variance increaSing with the spacing. The first figure is 
the general case for a seam with the same variogram; the 


second one iS a particular case, which may be easier to use. 


The krige variance is always less than the extension 
varliancet, except where both are essentially zero at sample 
Spacings less than 50 m. This means that a kriged estimate 
1s always more accurate than the arithmetic mean, down to 


50 m spacing. 


If the arithmetic mean is used, which closely 
approximates the calculation method used at the mineft, then 
there 1S a relatively rapid improvement in the reserves 
estimate as the hole spacing is decreased from 220 to 70 m. 
Within this distance the spacing can be chosen to trade off 


costs against estimate precision. 


The optimum spacing drill hole spacing must be chosen by 
the mine engineer using this curve and how much he is 
willing to pay for a reserve estimate of a given precision. 
+ The extension variance is the variance of the extension 
estimate, which is the arithmetic mean of the data. 


TALOPENSCCL  ODUD 1. 2et OlscmOeGchLDt LOReOGmarie resenves 
calculation method used at the mine. 
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At the drill hole spacing of 100 to 150 m used in the mine 
(Figure 2.2, p.40), the extension estimate has about ten 
times the variance of the krige estimate with the same hole 
Spacing. The same accuracy as the extension can be achieved 
by kriging samples approximately four times further apart, 


at 400 to 800 m spacing. 


Minimizing the estimation variance is not the only 
criterion for choosing the hole separation. Closely spaced 
Grilling 1S required at the Boundary Dam Mine to find and 
delineate coal lost to paleochannels and unrecorded 


underground mines along the Souris River Valley. 


The variance/sample spacing curve depends only on the 
variogram. The drill hole variogram is only defined for 
hy <"2000~m (Figuré” $35;"p.61) , so" thrs*variance/sample™ spac-— 
ing curve says nothing about sample spacings beyond 2000 m. 
In fact it must rise, otherwise one drill hole would give as 


good an estimation variance as a 2000 m grid of holes. 
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5. CALCULATIONS OF PRODUCTION AND RESERVES 


5.1 GENERAL METHOD 
In this study calculations were done on two sets of 
tonnage reserves: 

1. this year's in-place coal, or the amount mined this 
year, which is used to estimate the coal recovery, 

2. next year's in-place coal, or the amount to be mined 
next year, which is used to estimate reserves and future 
production. 

In both cases the krige estimate (Zf) and the extension 

estimate (Z%) were calculated and the results compared to 

determine if the extra effort involved in geostatistical 
calculations 1S worth the increase in accuracy. Both 
tonnages were estimated for the same domain to compare the 
actual production with the forecast production. Finally, 
the sources and the relative sizes of the errors were deter- 


mined. 


5.1.1 Calculation of the Tonnage Estimate and Its Variance 
The tonnage (T) is the product of the area, the density 


and the average thickness, or: 


T= sA<o 37" Ga 


Where: 
A = the area of the domain for which the tonnage is 


calculated. 
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the average density of the coal. 


z* = the average net thickness of the seam, either Z{ or 2%. 


The variance of the tonnage is needed for several 

reasons: 

1. comparing the accuracy of different methods for 
calculating the tonnage, 

2. determining the sources and sizes of the errors in the 
tonnage calculations, 

3. determining the range within which to expect next year's 
Production to fauna ; 


The tonnage variance is estimated by: 


VARUD). = G2g0°VARTAl + A~Z “VARLOY t A-d.VARUZ. | nen, 


Where: 


VARL ZS Jl = a, Or igs depending, on whether Z” = Z; or 2 


xx 


This equation is developed in Appendix 9.1. 


5.1.2 Calculation Methods Used at the Boundary Dam Mine 
Reserves at the mine are calculated by the polygon 

method using drill hole data. In this method the seam is 

broken down into 'polygons of influence" around each hole, 

with all the polygons meeting halfway between the holes. 

The seam is assumed to be the same thickness throughout the 

polygon as at the hole. Total coal volume is the sum of the 


products of the area and thickness for each drill hole 


polygon. 
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In-place coal is calculated using the highwall 
thicknesses and a modification of the polygon method. The 
drill hole polygons are used, but the polygon thickness is 
an average of all the pit meaSurements inside the polygon, 
rather than the drill hole thickness. From there, the 
method is the same as for the drill holes. Recovery is 
calculated as the percent ratio of the coal-loaded-out to 


the in-place coal. 


The calculation method used in this study is similar to 
the method used at the mine, except that the average 
thickness in this study is estimated twice (Z% and Z*) and 
the intermediate step involving the polygon of influence is 
bypassed. Since the pit data is fairly regularly distribu- 
ted (Figure 2.1, p.39) and since the miner's average is 
weighted on the areas of the polygons, the extension esti- 
mate (Z*) will give essentially the same answer as the 


method used by the mine engineers. 


5.1.3 Calculation of the Density Estimate and Its Variance 

The mine personnel have determined a ‘tonnage correc- 
tron factor’-of724°rt*/ton, “which 1s the reciprocal of the 
density (1.48 tonnes/m*), but they have not calculated the 


@rrop- or ‘this number. 


The error in this density value was determined by 


recalculating the density using the original data collected 
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by the mine engineers. They obtained these numbers by 
carefully measuring the volume and tonnage of the coal- 
loaded-out from three small (1x10* tonnes) sections of 
benches scattered over the mine. The variance of this data 
was eStimated using classical statistics and weighting by 
the tonnage of each sample. The density and its error is in 


Tabhe 1501. 


The coefficient of variation, sometimes called the 
relative standard deviation, is the Square-root of the 
relative variance, or the estimate divided by the estimation 
Standard deviation, and always expressed as a percent in 
this study. Since it is dimensionless, it is a straight 
forward way to compare the error in one quantity to the 


error in another quantity with a different dimension. 


Density 1S a regionalized variable, so proper analysis 
would estimate itsS mean and variance using geostatistics. 
Unfortunately, the amount of data is far short of the 
required minimum of 30 samplest and barely enough to to do 
classical statistics. As result, the variance estimated 
above is probably too high, since the method makes an 
implicit assumption of a pure nugget effect model for the 


variogram. 


+ For a brief discussion of this limit see Appendix 5.1. 
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Table 5.1 Data and statistics for the density estimate 


and error 


Density data 


Location Weight Volume Density 
(tonnes) (m?) (tonnes/m?) 

Ba teria wets ; 4723 3.133 leon 

bottom coal 

Phe df fcut e332; 14894 104982 1.46 

top coal 

Pires te Cut 47, 9570 6375 1ees0 

Density statistics 

Mean: 12285 tonnes 7m. 

Variance: 0.000435 (tonnes/m°) 7 

Goeclircient Of variatrvon: eran ice 

Relative variance: 0 40'20% 


+ A thin (0.1m) parting was removed from the volume and 


weight of pit 3, cut 4. 
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The density 1S assumed to be constant over the whole 
study area, so the same value is used for both this year's 
and-=next+year- s’tonnagesss Infact, "At* wiil* not be> constant, 
due to varying ash content, but there is not enough data to 


determine its spatial variability. 


5.2 IN-PLACE COAL AND PRODUCTION FOR THIS YEAR 
This year's in-place coal is needed to calculate the 
coal recovery, which is used by the mine staff to measure 


the efficiency of the operation. 


5.2.1 Domains Chosen 

The domains for calculating this year's in-place coal 
and recovery were selected on the availability of coal 
recovery data from the mine. Several sections of the mine 
had data on the coal-loaded-out as well as pit measurements, 
but the area of the domain could be determined accurately 
for only some of these. Dates on the mine plan, showing the 
location of some of the bench advances, often did not corre- 
spond to dates tabulated for the coal-loaded-out and 
therefore, for most benches, a direct comparison was not 
possible. Unfortunately, the rate of advance of the face is 
too erratic to make an accurate interpolation. For the 
remainder, a trial recovery was calculated, which should be 
in the range of about 70% to 100%. The recoveries fell into 
two groups: one within this range and the other grossly dif- 


ferent. The second group was either less than 50% or, more 
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commonly, far more than 100%, indicating that the wrong 
limits of the benches must have been used; consequently, 
these calculations were ignored. The following locations 
(Figure 5.1) had recoveries between 80% and 100% and 
therefore, were used for this study: 

1. pit 1, 1977 benches, 

2. pit 1, 1978 sbenches, 

3. pit 2A, 1977 and 1978 benches less the box cut, 


eee pits) DeNnchay. 


5.2.2 Calculating the Average Thickness 

The program modified from Flint (1978) calculated both 
a kriged average thickness and an arithmetic average 
thickness from the pit meaSurements; each is listed in 
Table 5.2. The average difference between the two estimates 
is about 1%, which is close to the size of the coefficient 
of variation of the estimates. Also the calculation error 
in the estimation variance (about 10%) is not much smaller 
than the difference between the coefficients of variation 
for the two estimates. As a result, the krige estimate does 
not make a large improvement over the extension estimate. 
The percent error of the krige estimate due to numerical 
approximations (that is the probable maximum difference 
between the calculated and the true value of the krige aver- 
age) is about 0.05% (see Appendix 8.5). Such small errors 
of approximation, almost two orders of magnitude less than 


the estimation error, can be ignored. 
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Table 5.2 Estimates and variances of the average seam 


thickness for this year's in-place coal 


Domain Estimate (m) 
Kriget lneb hoeye 
PLeiayeiorl os OU Se | 
Pate OS 5 Ou 3.61 
PIC 2A, 3.64 3.63 
1977 -& 1978 
Prcos, 3.00 Se Gal 
bench 9 


Avga,coerincient. of variation 


+ Krige estimate of thickness. 
$+ Extension estimate of thickness. 


Coefficient of 
variation 


Krige Ext. 
0.62% 0.62% 
0.44% ipso 
Zea Ser 
0.4% 223% 
0% Pe! 
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The coefficient of variation of a kriged estimate cannot 
be interpreted as it would be in classical statistics where, 
Given a normal distribution, the true value will fall within 
lo of the estimate two thirds of the time. The experimental 
distribution of errors in geostatistics usually has a non- 
normal distribution, with larger tails than the Gaussian 
distribution (Figure 5.2). Nevertheless, the usual 20 
limits have generally been found to include the usual 95% of 


errors, according to Journel and Huijbregts (1978, p.50). 


The error in the estimates for pit 2A are larger than 
the others because the model used had a nugget effect 
instead of a short range structure. The model was changed 
in order to experiment with the effect of a different 
Structural interpretation. The estimate is largely 
unchanged but the estimation variance 1S increased by a 


EACLOE Or aDOUtT. CWO. 


5.2.3 Area of the Domains 

The area of each domain is also calculated by the 
kriging program. It takes as input the X and Y coordinates 
of the vertices of a polygon defining the domain of interest 
and calculates the area of this domain exactly. The areas 


of the domains are listed in Table 5.3. 


Thesonly source Of errorain the area is, in the errors of 


digitizing the vertices, about 1 m tio, which is discussed 
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Experimental distribution of errors 


Normal distribution 


= 206. mean +1200. 


From Journel and Huijbregts (1978, Figure II.13, p.50) 


Figure 5.2 | Idealized distribution of the errors for 


most geostatistical estimates 
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Table 5.3 
Calculations 
Domain Area 
(m=) 
(xy 102) 
Pit 1, 1977 benches 2008 
Pit 1, 1978 benches les 2s 
Pita. total 3.836 
Bite 2A; A O23 
1977 & 1978 benches 
Pitas. bench & Chaz bso 
Paes total 1.555 


Areas of the domains used in the 
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time 
(eyata) 
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in Appendix 2.1.1. The relationship expressing the error of 
the area as a function of the errors of the vertices was not 
used, but instead a subjective guess at the error was made. 
The precise error was not necesSary since it is very much 


smaller than the errors in the density and thickness. 


It 1S possible to put an upper limit on the error of the 
area. If one assumes that the errors at all vertices are 
perfectly correlated, and act together to make the polygon 
larger or smaller at the same time, then the coefficient of 
variation of the area of a typical domain in this study is 
about 5%. But the errors at the vertices should be 
independent, not correlated. Most polygons in this study 
have about 100 vertices and this number of independent 
errors will tend to cancel, making the error in the area 
much smaller than the upper limit given above. In actual 
mine practice, the error of the area would be less than the 
estimates in’ this study since there would be no digitizing 
errors, only surveyor's errors, which should be much 


smaller. 


In order to determine the error of the area more 
closely, a subjective guess at the order of magnitude of the 
coefficient of variation was made A conservative estimate of 
0.1% was chosen, which probably overestimates the true 
value. This produces a relative variance of 1x10°°%, effec- 


tively zero when compared to the other relative variances, 
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which are uSually of the order of 1%. Therefore, 


equation {5.2}, p.84 reduces to: 


VARI Eo NAR Givd- eOVARL 2 177Z "2 {5.3} 


5.2.4 Calculation of the Tonnage 

The tonnage is calculated with equation {5.1}, p.83, 
using Z™ = Zk or Z%, depending on whether the krige or 
extension eStimate is used. The tonnage variance is 


estimated with equation {5.3}, p.96 slightly rearranged. 


The results of the calculations are listed in Table 5.4. 
As with the thickness, there is no real difference between 
the tonnages eStimated using the krige thickness and the 


extension estimate. 


Using the kriged average of thickness the coefficient of 
variation of the tonnage is 1.9%. Using the arithmetic 
average, the coefficient of variation 1S 2.3%. The 
difference between these two numbers (0.4%) is close to the 
approximation error so kriging is barely more accurate than 


the arithmetic average. 


The relative variance of the tonnage is the sum of the 
relative variances of the thickness and the density 
(equation {5.3}, p.96). The area contributes essentially 


nothing since its variance is so small. The density 
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Table 5.4 Tonnages of this year's in-place coal 
Domain Estimate Coeffiicirent of 
(tonnes x10°) variation 
Kriget Ext. Krige Ext. 
Pit 17,297 PETS Se 13.4 PRS Lpioy 4 
PRES TH; SOTLOTS OYE 51 6.69 mene EAA! 
Pites2A; Lo39 10.9 Pash § Sise 
1977 & 1978 
Pep 3, hes2 Vaso hays PA Ti 
bench 9 
Avg. coefficient of variation WD % Zu3% 


+ Calculated using the kriged estimate of thickness. 
+ Calculated using the extension estimate of thickness. 
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contributes 80% of the relative variance of the tonnage 
using the krige estimate and 30% to 60% of the relative 
variance of the extension tonnage. Thickness accounts for 
the rest. Pit 2A has a large eStimation error because of a 


change in the variogram model as explained on page 92. 


5.2.5 Calculation of the Recovery 
The recovery (R) is the percentage ratio of the coal- 
loaded-out to the coal measured by the pit data.t It is 


calculated by: 
R = 100°CLO/T (5/149) 


Where: 

CLO = the coal-loaded-out. 

The coal-loaded-out is the tonnage of coal coming out of the 
mine as weighed at the scales. Values for the benches used 

inthis section-are “listed iniTable .5.5., ;The formulasabove 


is discussed further in Appendix 9.2.1. 


The relative variance of the recovery estimate is 


estimated by: 
VAR[R]/R? = VAR[T]/T? {5.5} 


which is developed in Appendix 9.2.2. This equation means 
that the relative variance and the coefficient of variation 


+ The mine engineers also calculate recoveries using the 
Grill holes. But these were not considered, other than to 
note that they were generally worse (more variable and often 
over 100%) than the recoveries using the pit measurements. 
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Table 5.5 Tonnages of coal-loaded-out from the mine 


records 
Domain Coal-loaded-out 
(tonnes) 
Pit 1, 1977 benches be A084 567 
Pit 1, 1978 benches 587 774 
Pit 2A, 1977 & 1978 benches 1.092.569 


Pits oye bench,? 1288760 
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are the same for both the recovery and the tonnage of 
reserves. Thus, all the remarks above about the error of 


the tonnage also apply to the recovery. 


The recoveries from this study are listed in Table 5.6, 
which also has a comparison of the recoveries calculated by 
the mine personnel. The percent recovery averages about 90% 
according to all methods of calculation. The recovery 
calculated by the mine personnel using traditional methods 
is often within io of the geostatistical estimates and 
always within 20, except for pit 2A because of the nugget 
effect used in the variogram model. Again, there is little 
improvement uSing a krige estimate except for a slight 
reduction of about 0.4% in the coefficient of variation. 
Even the relatively simple calculation mehod used at the 
mine is more complicated than necessary since there is 
little difference between their results and a simple aver- 


age. 


Pie oliS oc. .1978) «coy 1008 (1977) elow) 100%) Sovalmost 
certainly the recovery is not perfect. The other two pits 
are within 1 or 20 of 100% so all the coal may have been 
recovered. The recovery from pit 2A is slightly above 100% 


due to statistical fluctuation. 
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Table 5.6 Comparison of the percent recovery 
calculated by the mine personnel with the 


percent recovery estimated in this study 


Domain Minet Estimate (%) Coefficient of 
(%) variation 

Krige¢ ExestT Krige Ext. 

PPE st 83720 S250) G27 125% WES HK 

O77 

Pie ad. 88.16 90:03 87.9 Bus oly 4 a7 % 

1978 

Bren2n, 90.91 HOO SZ LOO SZ 259% 333% 


oie) ee OO 


Pit 3; 96, 62 Shes: 94.7 120% Lik 
bench 9 
AVG se COGLINCTENt Ole Variatizon Vo9% Phe 


+ From the mine records as calcultaed by the mine engineers. 
$+ Calculated using the kriged estimate of thickness. 
+t Calculated using the extension estimate of thickness. 
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5.3 RESERVES AND PRODUCTION FOR NEXT YEAR 

Next year's tonnages and variance are also determined 
Byecouation 5,416 p.cs wand. equation {5.2}, p.84. The 
density is the same as that used for this year's in-place 
coal (estimated on page 85 above), but the thickness and 


area are arrived at differently. 


The reserves are estimated for a domain already mined 
out, but as if production had not yet started. The 
estimation used only data that would have been available 
before excavation, either pit data (from nearby benches) or 
Grill holes. These tonnage estimates can then be compared 
with the 'true' values, which are estimated by the best 
method possible (using kriged thicknesses of the highwall in 


the pit). 


5.3.1 Domains Chosen 
The domains used when estimating with the drill holes 
are all three pits (Figure 1.11, p.24): 
t. Dit lyre! Sige and 1978 \benches, 
2. Pato 2A, 
3... Wits. | 
The domains used when estimating with the pit measurements 
Bre west wouneives of pit 1 (Figure 5.1, p.90): 
f.. Oat bln-bOve -pbenches, only, 
2. Wa tether nl SS benches, only, 


All these domains were in production for approximately one 
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year (Table 5.3, p.94); as a result, they demonstrate the 
errors to be expected when estimating in-place coal or pro- 


duction for some future year. 


5.3.2 Calculation of the Thickness 

The average net thickness used to estimate next year's 
production can come from two sources: either drill holes 
into and around the projected pit, or pit measurements from 

a highwall beside or near the projected pit. Estimates 

using drill holes have the following @ priori advantages. 

1. Any drift will affect the results less since the 
estimates are interpolated rather than extrapolated as 
the pit meaSurements must be. 

2. The estimates can be made wherever there are drill 
holes, not just near a previously opened pit. 

3. There are fewer drill holes, so the computing will be 
cheaper. 

Estimates using pit measurements have the following @ priori 

advantages. 

1. They are less biased since they have approximately the 
Samesmeaneacwine § crue aivalue. (rable: 2. 1 sp.5d) 5.0 
Cont rastmto thes driv holes): 

2. There are more meaSurements, so the estimate may have a 
lower variance. 

3. The benches are generally advanced in a regular manner, 
so there is almost always an open bench near a proposed 
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5.3.2.1 Next year's thickness using drill holes 
The average net thickness for next year's in-place coal 
is determined from the drill hole measurements that are in 


Or near the proposed pit. 


Both the kriged and extension thickness were calculated 
and compared. As for this year's in-place coal, the 'true' 
thicknesses were determined by kriging all the pit measure- 
ments taken within the pit used as a projected domain. 
Kriging was used because it should give the best (in a geo- 
statistical sense) determination of the actual thickness 
against which to compare the projected estimates. The 
coefficient of variation is also tabulated to show how much 
latitude there is in the 'true' value; usually about one 


third to one half of the variability in the estimates. 


The results of the thickness calculations are listed in 
Table 5.7. Again, there is little to choose between krige 
and extension estimates except a slight increase in the 
estimation variance. However, both estimates consistently 
underestimate the ‘true’ thickness by 1.50 to 20 or about 6% 
on average. The 1o range on the -6% figure is +3%, which 
means that the bias is probably real. The bias 1s due to 
the drill hole thicknesses being thinner than the pit meas- 


urements on average (Table 2.1, p.51). 
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Table 5.7 Comparison of the 'true' thickness with the 


estimated thickness for next year, 


calculated using drill hole data 


Domain Estimate, (m) 


Krige Ext. 


Pith i. Covad 
'"True'T S52 


Drill holes Ep esis is Spee We, 


Pitech Iino), & 1978 
BT RUe 2a] ed A, 


Drill holes Smo Silon!| 
Pites, ;totow 
hed ts og BY a Slee 


Drill, noes Shey ZA Sez 


Avg. coefficient of variation 


Avg. error ral 8 ce fare Wy Saltalea Sysi sere ip 65 


Coefficient of 
Variation 


Krige Bx 
0.99% 
ox 2.4% 
Oey 
2e9% 3.4% 
0.48% 
Sr oy 4.8% 
Zeik 3.5% 


+ 'True' thickness is kriged from pit measurements within 


the pit. 


ac «wt 


3x2 


iss 
- L} 
i. 


to svateyraak® 


tb 


noiteiaey « 
egizt 


noissitev io Sag 


#8. meee 


=| 


- A, ~ i 
~CIG.0> 


metvese@ Jig morl bepizA-si eseasore: 


106 


5.3.2.2 Next year's thickness uSing pit measurements 
The same domains were eStimated as above but this time 
using the thickness estimated from pit data of the benches 


near the domain. 


The thicknesses of domains in pit 1 were estimated using 
two types of benches: 
1. all the benches excavated in the previous year (the pre- 
vious year's pit), 
2. just one bench, which was considered for two cases: 
A. the bench immediately beside the proposed pit, bench 
38. Since this bench straddles the pits for the two 
SuccesSive years, only the half that was outside the 
proposed pit was used. 
B. a bench more removed from the proposed pit, benches 
36 and 39% 
These benches and pits are shown on Figure 5.3. The average 
thickness along with the 'true' thickness determined from 


the measurements made in the pit are shown in Table 5.8. 


In contrast to the drill holes, the estimates using 
nearby pit data are unbiased. The actual examples presented 
in Table 5.8, p.108 tend to overestimate the true thickness 
by about 4% but the lo scatter of the estimates (about +5%) 
includes the 'true' value, thus the apparent bias is proba- 
bly just a, candom fluctuation. Using data from a full year 


reduces the coefficient of variation by only 0.5% compared 
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BENCH 38 
BENCH 3S 


BENCH 36 
BENCH 38 


Sen Lamhe PS Wie: 


0 500 


metres 


Figure 5.3 Benches used to eStimate next year's in- 


place coal and production 
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Table 5.8 Comparison of the ‘'true' thickness with the 
estimated thickness for next year, 


calculated using nearby pit data 


Domain Estimate (m) Coefficient) of 
variation 
Krige Ext. Krige Ext. 


Patan, 199777 
'True't 3.50 0.62% 


Estimated by: 


1OVeR pit 3.66 B08 1Wa6% 1 0% 
Bench 38 Sic 3.94 2100;% 2heaux 
Bench 39 S87 3335 PLD Say § Orne 


Pit iene O78 
‘True! + 3.51 0.44% 


Estimated by: 


WOVE Dat 36. 6F Cay 3S) 1.0% 1.4% 
Bench 38 Sav 3.94 2260% Daan 
Bench 36 Say S58 1.7% W. SX 
Avg. coefficient of variation VDE Woe 


Avg. error 3.8+4.4% 4.7+6.6% 
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to the data from one bench. Also, it makes no difference to 
the coefficient of variation if the bench is beside (bench 
38) or removed (benches 36 and 39) from the proposed pit. 
This insensitivity to the location of the data is because of 
the lack of a large drift and because the distance 
separating data from the pit to be estimated is generally 
between the long and short range structures of the vario- 
gram. This portion of the model has a low slope 

(Figure 3.4, p.60), so the 7 terms and the resulting esti- 
mate will not change much with a change in sample separa- 
tion. Finally, it makes little difference if the krige or 
extension estimate is used since the improvement using geo- 


Statistical methods is less than 1%. 


The to of the distribution of actual estimation errors 
is 5%, whereas the average coefficient of variation is 2%, 
for both krige and extension estimates. These two numbers 
Should be the same; the difference suggests that the 
calculated estimation variances underestimate the 'true' 
estimation variance. However, the discrepancy may be a 
Statistical fluctuation since it is close to the errors that 
are introduced by the calculation approximations. If the 


discrepancy is real, then there could be a local drift. 


5.3.3 Area of the Domains 
When considering the area and its variance for next 


year's tonnage, there are two possibilities. Firstly, the 
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area can be assumed to be known perfectly. That is, given a 
domain, not necessarily for one year's production, what are 
the reserves underlying it? The domain could be an 
arbitrary block of land, a lease limit or some other 
predefined area. This domain would be used to estimate the 
in-place coal for next year. Alternatively, the area is 
what will actually be opened up next year, which is subject 
to considerable error. This domain would be used to esti- 


mate the projected production for next year. 


For the first possibility, the area is perfectly known 
Wien NO.ertob, in particular VARL[A] = Os8and equation 15.33, 
p.96 applies. For the second possibility, the variance is 
controlled by the efficiency of the excavation operation, 
hence very large. How to determine the variance in this 


second case 1S discussed in Appendix 10. 


5.3.4 Calculation of the Tonnage 
The tonnages are calculated in the same manner as the 
tonnages for this year's in-place coal. The tonnages 
calculated from the drill holes are listed in Table 5.9 and 
the tonnages from the pit measurements are in Table 5.10. 
In each table, the error of the tonnages is estimated for 
the two cases for the variance of the area: 
1. VAR[A] = 0, when the area is predefined, for example, 
next year's in-place coal, 


2. VAR[A] > 0, when the area is what will actually be 
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Table 5.9 Comparison of the 'true' tonnage with the 
estimated tonnage for next year, calculated 


using drill hole data 


Domain Estimate Coefficient of variation 
(tonnes x10°) VAR[A] = 0 VAR[A] > 0 
Krige Ext. Krige Ext. Krige & Ext. 


Pitt 1, total 


“irae 5.42 Pas 4 34% 
Drie: 5s, 16 Sia6 XRG * 2 sary 3% 
hole 


DitercAy tod de Lo LO 


ar were 5.80 uae Se 
Drill 5.54 5.49 3.2% 37% 31% 
hole 


Patt Gc evotel 


'True't 52:56 l,. DS 31% 
Droid 5.02 4.95 3..9% 5. OF 31%, 
hole 

Avg. 

coefficient of- variation 3.0% eugtaled By 
Avg. 

error Pints ti 9 Oot oi Oe 


+ Calculated from Table 5.97, 1. 105.. 
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Table 5.10 Comparison of the 'true' tonnage with the 
estimated tonnage for next year, calculated 


using nearby pit data 


Domain Estimate Coefficvent. of variation 
(tonnes x10”) VAR[A] = 0 VARTA) > .0 
Krige Ext. Krige Ext. Krige & Ext. 


Pre 1977 
“SEY eLay 5. 39 LES % S15 


Estimated by: 


ROU SS Bite 5.63 556 Po eile 3% 31% 
Benen 38 _5.86 6.06 2.4% 5% 31% 
Bence our owt Sy he) 2.6% 6% Sa ae 
Prete i619 7.8 

“True T 5.40 anew 3% 
Estimated by: 

CAE oy ie mr asied els Seog hy | 0% Sle 
Beneh-38°).5.77 6.06 2.4% 7% 31% 
Benchs6_ 5.249 sees 22% 3% 31% 
Avg. 

coefficient of variation pleas 4% 31% 
Avg. 

error 3.84+4.4% 4.7+6.6% 


tT \Catculated from Table 5.87) p. 108. 
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Opened up next year, or next year's production. 

The 'true' tonnage is calculated from the 'true' thickness, 
which is kriged from the pit meaSurements inside the domain 
estimated. The area used is always the estimated area 
opened up in one year, 1.04x10° m? (see Appendix 10.3). 
Even though no pit has exactly this area, a constant area 
was used so that the tonnages could be compared easily, 


Giffering only in the thickness used. 


The error of the projected production (VAR[A] > 0) is 
dominated by the error in the area that will be opened up 
during the next year. The proportion of the relative 
variance of the tonnage due to both the density and the 
thickness is so small (about 0.5%) that virtually any data 


and any calculation method 1S good enough for estimating. 


The error of projected reserves (VAR[A] = 0) is similar 
to that for the in-place coal and recovery for this year. 
The density contributes about 20% to 40% of the relative 
variance, but slightly less if the extension estimate or the 
Grill hole data is used, since both have a higher variance. 
Also, the krige coefficient of variation is less than 1% 
better than the extension coefficient of variation. Once 
again there is little improvement in the estimates using 


geostatistical methods. 
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6. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


6.1 SUMMARY OF RESULTS 
6.1.1 Classical Statistics 

In the Boundary Dam Mine area, seam net thickness meas- 
urements are fairly regularly distributed on square grids of 
30 m for the pit data and and'100, 150 and 200 m for the 
Grill holes. 


Average net thickness of the seam in the pits is 3.6m 
but thinner (3.3 m) when measured in the drill holes because 
more splits were logged than were selectively mined. The 
Grill holes have about twice the coefficient of variation, 
probably due to more variable logging styles on the part of 
the drill geologists and more consistent reporting by the 
mine Surveyors in the pits. Both kinds of data are approxi- 
mately normally distributed with tails on the thin side, 
possibly due to paleochannels cutting into the coal. The 
tail on the drill hole distribution is large enough that a 
three-parameter lognormal distribution may fit better. The 
pit data could be separated into two populations (controlled 
by the partings), primarily on statistical evidence and 
secondarily on geology. On the other hand, the drill hole 


data was one homogeneous population. 
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6.1.2 Geostatistics 
No drift waS apparent in either set of data, except 
possibly at distances greater than the largest distance 


(about 2000 m) used in the kriging. 


Transition variogram models with no drift were fit to 
the data. The best model of net thickness for both pit and 
Grill hole data is the spherical model, usually as two 
nested structures. A hole effect model may be a better fit 
to data that includes the thickness of any rock, in particu- 
lar, the partings and the gross thickness. There is also 
the possibility that a Gaussian model may fit the net 
thickness. Fortunately, kriging is insensitive to the model 


chosen so the spherical model gave good results. 


No nugget effect was used, though it is possible to 
model the variograms with a small nugget effect. This 
change in the model has little effect on the estimate but a 


moderately large effect on the estimation variance. 


The sill is approximately 0.07 m?* for the pit measure- 
ments and about four times as high for the drill hole meas- 
urements. This higher variability is probably due to the 
random inclusion of thin rider seams and other differences 
in style of interpretation between geologists logging the 


holes. 
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Ranges from this study are of the order of several 
hundred metres. These, along with the ranges found by Flint 
(1978, p.107) form an approximate geometric progression: 
5O0im: W50empon500sm,242000em Cand> 12/000! m? jee (The ‘structures 
are probably due to depositional factors, though the 50 m 


range could be due to measurement methods. 


Pit variograms show a small amount of zonal anisotropy 
perpendicular to the benches. The anisotropy ratio is about 
1:2 for both the sills and the ranges, with the larger 
values lying in the direction perpendicular to the bench. 
This is probably due to the mine Surveyors remembering the 
last measurements they took along the bench when confronted 
with a gradational contact. The drill hole variograms are 
isotropic, as is the variogram for the aggregate of all the 


pit data. 


6.1.3 Drill Hole Spacing Strategy 

Sampling strategy is determined only from the point of 
view of minimizing the estimation variance, without 
considering the necessity for outlining mining hazards. 
Sample spacing is a trade-off between cost and estimation 
error. Surveying the highwall is comparatively cheap, so 


the pit sample spacing was not analysed. 


The error of any estimate will increase as the sample 


Spacing increases. Figure 4.2, p.81 illustrates this 
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relationship for the drill holes at the Boundary Dam Mine. 
If the extension estimate is used, essentially what is done 
at the mine, then the extension variance is only affected by 
drill hole spacings between 70 and 220 m. Beyond these 
limits the variance curve is flat with respect to spacing, 
which can therefore be maximized with no effect on the error 
of the estimate. On the other hand, the krige estimate has 
a lower variance and a steady trade off between spacing (or 
cost) and the error over the full range of practical sample 


Spacing above 70 m. 


At the drill hole spacing used at the mine (100, 150 and 
200 m), the extension estimate has about three times the 
error of a kriged estimate using the same data spacing 
Conversely, for the same estimation variance, kriging 
requires samples four times further apart than the simple 
average of the data, the extension estimate. However, as 
will be seen in the next section, the percent errors are so 
small (1% to 5%) that there is not much practical difference 


between the two estimation methods. 


6.1.4 Tonnage Calculations 

Tonnages in this study were calculated two ways: by 
kriging and also by a simple arithmetic average of the 
thickness, similar to the methodt used at the mine. As a 
result, the conclusions from this study should be applicable 


+ The mine engineers use a modification of the 'polygon of 
influence'. 
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to the regular method of calculation used by the mine 
engineers. Tonnages were calculated both ways for two 
cases: tonnage of coal that was in the pit, but is now 


recovered; and the tonnage of coal in a proposed pit. 


6.1.4.1 In-place coal and recovery for this year 

Conclusions for the reserves calculations are the same 
as those for the percent recovery, because the coal-loaded- 
out has negligible error. The recovery averages about 93% 
and ranges from 83% to 100%. Since the coefficient of vari- 
ation for the recovery is about 2%, most of these are a 


Statistically significant reduction from complete recovery. 


Only a small improvement in the error can be made with 
the more sophisticated methods of geostatistics. The 
extension estimates (the arithmetic means) have only a 
Slightly larger error (2.4%) than the krige estimates 
(1.9%). In all four areas, recoveries estimated both ways 
fFaliewithine 1-or 2¢"0f each other... Considering the extra 
effort of geostatistics, there is no advantage in kriging 


the estimate when a simple arithmetic average will suffice. 


Density variation contributes 80% of the error in the 
tonnage or recovery when a kriged estimate is used and about 
half of the error when an extension estimate is used. The 
error of the average net thickness contributes the rest 


since the area is essentially error-free. 
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6.1.4.2 Reserves and production for next year 

No large advantage iS gained by using sophisticated 
calculations of average thickness and tonnage. The simple 
arithmetic mean of the thickness is adequate for calculating 
reserves. Kriging only reduces the estimation variance to 
3% or less, from about 4% for the extension estimate. The 
reserves calculated using meaSurements from a nearby pit are 
an unbiased estimate of the 'true' tonnage to be expected in 
the projected pit. On the other hand, the drill holes 
consistently underestimate the the 'true' tonnnage by about 
6% due to the drili hole net thickness being thinner on 
average than the pit meaSurements. The 'true' tonnage was 
determined by kriging highwall thicknesses from the pit. 
The largest error by far 1S introduced by the variable 
efficiency of the men and equipment, which can expose at one 


face somewhere between 7x10*% and 13x10 m? 


of ground per 
year. All other sources of error in the projected tonnage 
are negligible by comparison. On the other hand, if the 
area to be mined is known without error, then the density 
contributes about a third of the relative variance with the 


rest coming from the thickness. Calculation errors, such as 


the numerical integration, are minute compared to these two. 


6.2 CONCLUSIONS OF THIS STUDY 
The following conclusions strictly apply only to the 
studied portions of the Boundary Dam Mine (pits active in 


1977 and 1978). They can be applied to other portions of 
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the mine if the properties of the data (particularly the net 


thickness variogram) do not change too much. Possibly the 


conclusions can be used at other plains coal mines with sim- 


ilar data properties, though this must be done with care and 


the results should be tested. 


he 


The error of the tonnage estimates increases when drill 
hole spacing increases from 70 to 220 m. Below 70 m 
Spacing the error is negligible, while beyond 220 m and 
out to at least 2000 m the error is constant at about 
10% and unaffected by the drill spacing. Between these 
limits hole spacing can be chosen from Figure 4.2, p.81 
for the desired error. Other considerations, such as 
searching for potential mining problems, may determine 
the spacing actually used in practice. 

Methods of determining reserve estimates could be 

altered slightly to reduce bias and also simplified 

without seriously increasing the error. 

A. Drill holes produce a reserve estimate that is 
Slightly biased on the low side (by about 6%) and 
with a larger variance. The best data to use are 
highwall measurements from a pit or benches that are 
no more than a few hundred metres from the proposed 
pits 

B. Geostatistical methods of estimating reserves of 
coal at the Boundary Dam Mine are only marginally 
better than the simplest method, which is arithmetic 


averaging of the data, essentially the method used 
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at the mine. The main advantage of the more sophis- 
ticated methods is a slightly smaller estimation 
variance, so it 1s doubtful if the large effort of 
geostatistical estimation is worth the minor 
improvement in the accuracy of the estimate. Geo- 
Statistics are only useful if the estimation 
variances are required, as they were in this study. 
Errors in the tonnage and recovery estimates are 
equal and very low, about 4% to 6%. These are 
approximate 95% confidence limits, so that any 
change in the recovery of less than 4% has only a 1 
in 20 chance of being real, rather than a random 
Staristacalsrluctuations 

If the area of the proposed pit is given, and 
therefore error-free, density and thickness 
contribute about equally to the error of the 
tonnage, but the proportion of each can range from 
one quarter to three quarters. Calculation 
approximations are miniscule by comparison. 

If the area of the proposed pit is unknown and must 
be estimated from the area that can be uncovered by 
the stripping operation in one year, then the error 
Of this area.is,sovlarge, (about, 30%) that wit 
overwhelms all other sources of error. In this 
case, any data and any calculation method are 


adequate. 
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6.3 RECOMMENDATIONS 

This section contains two sets of recommendations; 
beginning with further research that could be done in the 
area of this study, followed by suggested changes to the 
method of calculating reserves and recoveries at the 


Boundary Dam Mine. 


6.3.1 Further Research 

Research that could be done in the future falls into 
two areas. Firstly, a more thorough geostatistical analysis 
of this data or of the Boundary Dam Mine in general could be 
done, primarily to check some of the assumptions that had to 
be made. Secondly, more general geostatistical work could 
be done. This would explore some of the possibilities that 


were not pursued in this study. 


6.3.1.1 Detailed research to extend this study 

This study was detailed enough to determine the best 
method for calculating reserves and expected production; 
nevertheless, it was only a preliminary analysis of what 
would be a more complete geostatistical analysis. A much 
more involved study could be done to determine the validity 
of the assumptions and the Sensitivity of the results, not 
for practical purposes but for the intellectual exercise. 
Suggestions for potential studies are listed below. 
1. Regress the drill hole data against estimated 


thicknesses at the drill hole locations. The estimated 
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value should use the pit data in order to reduce the 

bias. This regression could depend on: 

A. the location within the mine, 

B. the number of splits, 

C. the drilling program, either preliminary or infill. 
Contreabling formtherdniliing programemaycontrol.<for 
the different drill geologists and their variable 
logging styles. 

The drill holes, corrected by the regression, may 

produce a less biased estimate though the estimation 

variance may not be improved. 

Study the thickness variogram in more detail. 

A. Take detailed measurements in the pit to eliminate 
the possibility of a nugget effect and to determine 
if the parabolic behaviour at the origin is real or 
not. 

B. Consider anisotropic models for each pit separately. 

C. Use other models, such as a Gaussian or a hole 
effect model. 

D. Check the intrinsic hypothesis by gathering and 
analysing data from more widely scattered portions 
of the mine. Does the variogram depend on location 
or the sedimentology of the floor, seam or roof? 

Examine more closely the kriging and its effect on the 

results by considering the other methods to estimate the 

net thickness in Appendix 2.2. 


A. Use smaller domains with more local variogram 
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models. For instance, krige the area underlain by a 
Single split and combine several such areas. Next, 
compare this result to an overall estimate. 

Use universal kriging, particularly split by split 
rather than on the total thickness as was used in 
this study. The individual splits tended to show 
more drift and a more individual variogram than the 


aggregate thickness. 


Determine the density variogram and evaluate more 
precisely the effect of the density variance on the 
tonnage estimates. 

Determine the variability of the coal-loaded-out and its 
effect on the variability of the recovery. 

Determine explicitly the error in the area of a polygon 


and check that it really is small. 


6.3.1.2 General geostatistical research 


Below are listed some of the more general topics that 


could be done using this data, with or without geostatistics 


Bnaswith or without the’ data from Flint’ (1978)% 


Flint's kriging program could be expanded to handle 


other types of kriging: 


Ts 


universal kriging, in which case the drift could be 
compared to a trend surface analysis, 
lognormal kriging, 


punctual kriging, 
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4, random kriging, 

Bead licezs kriging, 

6. changing data and/or model parameters interactively, 

7.  non-spherical models, in which case the modelling pro- 
gram would also have to be altered. 

Universal and/or lognormal kriging could be used with the 

drill hole data, which exhibits a drift and may be 

lognormally distributed. If the kriging capability is 

expanded, then the variogram modelling program should be 


also. 


The ranges of the variogram model from this study and 
Flint (1978, p.107) are in an approximate geometric 
progression, as are the ranges for the Lorraine iron 
deposits in France (Journel and Huijbregts, 1978, p.167). 
The estimates of the ranges could be improved with a merged 
Gata set. Then, if this progression is not just chance, it 
Should be possible to predict the next larger and smaller 
ranges. There is already a Suggestion of the next larger 
range in Flint (1978, p.108) at 12 000 m. These two 
predicted ranges could be checked, though gathering the data 
would be difficult in both cases. Estimating the shorter 
range would require detailed meaSurements in the pits, prob- 
ably at several locations because the indeterminacy of the 
upper and/or lower contact may Gepend on the immediately 
adjacent facies. Estimating the larger range may call for 


unavailable or nonexistent data, necessitating the expansion 
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of the variable under study to include other laterally 


equivalent coal zones. 


Only one variable, the net thickness, was analysed geo- 
Statistically in this mine. Variograms could also be 
estimated for any of the variables in either the proximate, 
ultimate or maceral analysis. These could be compared to 
see if the model parameters, particularly the range, were 
Similar and if any regularities had a geological control. 
More detailed data, Such as proximate anlyses, could be 
cokrigedt with more scattered data, such as ultimate analy- 
ses, to improve the estimate of the latter. In particular, 
the relationship between the density and the thickness could 


be determined. 


Any of the above analyses could be done for several 
mines and any Similarities compared across zones, formations 
Or regions. In particular the variogram model could be 
compared, and if it turns out to be fairly constant, then it 
could be used for future exploration in mines where there 


was insufficient development to define a stable variogram. 


The data set used for this study is large and quite 
détailteds rangingespatiablymirom 30,.to 2500 m. If 1t canbe 
combined with Flint's data set, then the span is from 30 to 
14 000 m, which appears to include several transition models 


t+ Cokriging is the simultaneous kriging of two or more 
regionalized variables. 
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with ranges in a geometric progression. Such conditions can 
be modelled by one logarithmic scheme (Serra, 1968) and this 
could be explored. The data could possibly be fitted 

equally well with a generalized linear model (Figure 6.1) of 


the form: 
wih) s= ah” “<0. 0° <n < 2 {6.1} 


The generalized linear model with n < 1 is the theoretical 
variogram of an object with a fractal dimension, where the 
power (n) is the non-Euclidian dimension correction factor. 
A large span in order of magnitude is rare in most data 
sets, but is necessary for any fractional dimension analy- 
Sis, which has not apparently been calculated for a coal 
seam or any stratum in general. The variogram and the non-_ 
euclidian dimension correction factor are both measuring 
approximately the same thing: the variability of the data. 


It would be interesting to compare the two approaches. 


6.3.2 Suggestions For The Mine Engineers 
This section contains suggestions to the operations 
personnel at the Boundary Dam Mine for possible changes in 
the way reserves are calculated. These suggestions may also 
be applicable to other plains coal mines, particularly if 
they have the same variogram. 
1. Re-interpret the drill hole thicknesses to the thickness 
that will actually be mined, or determine a correction 


factor by regressing it with the pit data or coal- 
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From Clark (1979a, Figure 1.6, p.9) 


Figure 6.1 Generalized linear model 
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loaded-out. 

When estimating reserves, use thicknesses measured from 
highwalls up to 500 m from the proposed pit. Those 
numbers are less biased than the drill holes 

Include all the data inside the area to be estimated 
plus any data less than a fixed distance from the 
margins of the area. The distance is determined by the 
range of the variogram, about 500 m for the Boundary Dam 
Mine. 

Average the net thickness over the whole seam. It is 
not necessary to break it into splits for separate 
calculations. 

Use a simple arithmetic mean when calculating the aver- 
age thickness of the seam over a given area. Even the 
best method, kriging, is only a marginal improvement 
over the simplest method, an arithmetic average, if the 
Sample sites are fairly evenly distributed. 

Do not bother doing calculations to a precision better 
than about 1% (other than guard digits), since this is 
the lowest error in the input data (net thickness, area 
and density) and some errors go up to 5%, particularly 
the area. 

Refine the estimate of the density and its error, since 
it contributes about half of the error in most tonnage 
calculations and it should be fairly easy to lower the 
uncertainty. 


For projected production, refine the estimate of average 
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man/machine productivity, or how much ground will be 
opened by the operation in a given amount of time. This 
contributes so much variability that if the error in 
this area cannot be greatly reduced, then any data and 
any calculation method is adequate. 

9. Recovery is only accurate to about +4% for 95% of the 
time; therefore, changes smaller than this are probably 
not real but random statistical fluctuations. 

In general, the tonnage calculations can be simplified 


because the seam is so regular. 
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APPENDIX 1 
OUTLINE OF THE THEORY OF STRUCTURAL ANALYSIS 


The following is a brief introduction to elementary 
variogram theory, with the emphasis on those aspects 
required for this study. David (1977) gives a good intro- 
duction to the subject, while Journel and Huijbregts (1978) 
give a more complete description. For a definitive treat- 
ment one must read the papers of Matheron and his colleagues 
at Fontainebleau; for example Matheron (1972). 


1.1 INTRODUCTION 


The term 'geostatistics' is used in two senses: North 
American, or more general sense; and a French, or more 
specific sense. The North American usage of the word is 
analogous to the definitions of geophysics or geochemistry. 
It means the application of any method of statistics toa 
problem in geology. In contrast to this, the French usage 
of geostatistics refers to a particular method applied to a 
particular problem. The method is any application of 
Matheron's theory of regionalized variables. The problem is 
to determine the best estimate of the average of a 
regionalized variable over some area. 'Geostatistics' is 
used in this study in the French sense. 


1.1.1 Regionalized Variables 


A regionalized variable is one that is distributed over 
a region in space or time, which means it exists at every 
point in the region. In practical terms this means that all 
Samples have a location and this location has an effect on 
the statistical properties of the variable. Usually the 
location is in 2 or 3 dimensional space, but if the location 
is in time, then signal theory 1S more commonly used. To be 
usable in geostatistics the variable must also be linearly 
additive, which means any linear combination of the variable 
must have the same meaning as the variable. Also most of 
the theory has only been worked out for continuous measures. 
Thickness and density are regionalized, area 1S not. 


For the purposes of this study, classical statistics is 
defined as any non-geostatistical technique. Almost any 
Statistical study must make assumptions; unfortunately, 
classical statistics makes one that is almost never true in 
geology. It assumes that the data used is independent; in 
other words, any one data value is not correlatedt with any 
other data value. However, in geology the data at one loca- 
tion is usually correlated (sometimes very highly) with the 


twiComnelated’ sisvgenerally used in its statistical, not 
geological, meaning for the remainder of this study. 
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data at nearby locations. Classical statistics can give an 
unbiased average, if the sample distribution in space is 
very regular, but it will not be the best, except by chance. 


The best average is defined in terms of the BLUE, which 
is an acronym for Best Linear Unbiased Estimator. 


Best means that the variance of the estimation is a min- 
imum. This means that over the long term this estimate (2*) 
will come closer to the true value (Z) than any other esti- 
mate. The '*' indicates that the quantity is an estimate, 
usually of the thickness, while the symbol for the true 
value has no modifier. 


Linear means that the estimator is a linear combination 
of the data. Mathematically the linear estimation is 
expressed as: 


Zn = La 2 | ir ies 
Where: 
Zero che summation, for @li if= 1 to n,; n is' the number “of 
; samples. 
ea, = the ¢" tn'werqhnt4 
z, = the i'th sample. 


It 1S not necessary that the estimator be a linear 
combination; it can be any function of the data, such as: 


Deen pe eT a er Neate ene ne feo 


Where: 

Riza o cee ec, | = eSomerarbpitrary=erunctilon Of ald 
sample values. 

An arbitrary function will, in theory, give a better esti- 

mate with a lower variance. More often though, a simpler 

nonlinear formula 1s used: 


/ A PE il Waa Cleo 
Where: 
f,{z;] = some function of one sample value, the same func- 


tion is used for each sample but with different 
parameters. 
However, this method, nonlinear or disjunctive kriging, is 
G@ifficult and the theory is not yet fully worked out. 
Linearity is a mathematical convenience for ease of 
computation and which also requires fewer assumptions about 
the statistical properties of the data. 


An unbiased estimator 1s one that over the long term 
neither over- nor underestimates the true value. It is an 


+ Geostatistical convention is to use the abreviated form of 
the summation sign. 
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estimator because the objective is to statistically estimate 
an unknown average thickness. 


1.1.2 Fundamental Measurements in Geostatistics 


The two most important parameters of the distribution 
of a random variable are the first two moments, in classical 
statistics these are: 
iemeans pF 
2h iGvuariances of: 

The corresponding parameters of a regionalized variable are: 
fevedei ttl om(h)s 

20trvarioogram>.7¢h) < 

In geostatistics these are not defined on the variable 
denecuiy out On gts inerement, Z2(xth).—9Z(x). The adriftryrs 
defined by: 


mony etal 2 (xthy noes (x0) {4244 

Where: 

E[R] = the expected value of R. 

R = a general variable. 

x = the location expressed as a vector.f 
im panticuvarsex = [urvily a is sthe: easting and vy 
the northing. 

Z(%) = the value of the regionalized variable at location 
esc 

Z(xth) = the value of the regionalized variable at location 


xt+h,which 1s separated from location x by the 
vector h. 
The variogram, y(h), is defined by: 


Veh) TSCIVAR(Z(xth) =t20x4 io 


Where: 

VAR[R] = the variance of &. 

The drift and variogram are functions of the vector between 
each sample pair. Usually they are presented as a graph of 
the statistie, plotted ‘on the vertical, axis, against the 
Sample separation, along the horizontal axis, rather than 
the scalar representation used for the moments in classical 
statistics. 


1.1.3 The Intrinsic Hypothesis 


In order to model reality, particularly with mathemat- 
ics, it is always necessary to make assumptions about the 
mathematical behaviour of the real world. The main assump- 
tion made in the geostatistics used in this study is the 


+ The convention used in this study is to represent any var- 
tableethat cant bévalvectorgin boldface: 


~ 
Sisticiige 
a bale rye 
f B42 & 
r 
r 
: a t 
’ 
va 
= 
= aa 
J fed — | Caw Be 
¢ 
Sry. 2zuOoGcse 
_— —— — ~~ 
op Aad ead pape S 
7 ; 
an 4 r 
a) « oi “ 
= ‘ 
~Vev a6 3429 
a2 ¥u au~ 
= ; } 


r 


19. szepeinsyeg hie ind, 32 


(5, sa 
. gtaeeoM ana 28223. edd. 9 5 aided Bb: 
: see vu ane th ae 
becitaneipes a te 2 tes9mai6 
. y 
) ) ings 
/ Sad. no. Beni te Jon, st6 -@3en- 
(=) = tie) =. -tesamsatoai 2 


ea 
- “jer he otisv Bedoeg 
: eee Sees 


‘a if 
ie & 26) Geresia 
"S21. hee “sani 
> - 
«fits 
iv & ae | az kia. 
t{ooi €¢a3 30 
a } a G32 rae 
SigsesS-a: <\ 
i 5 Lam 
- é tee 
& 
sm 6190). SS) 
> «au-«2 od. =) Le ae T6& MEIC 
‘ F J — t oe 7 
fasaerd ets yeds ufieue 
« © . 
tea tagiitay eft ao Becta 
¥a 7 TOSLZION acid pacis: 
< in . > + > 
* Ke { one Noistss 


iasfusidxeq .villes? fete” 03 3a 
GWwWLS6 Geet oF YIseReren, BAS ae: 
| «biter ises ef2 ie “yoivaded. Gan 
side ai Bear scé: siaetaoeg ads ak es v 
% at ™ a 


7 O07 at yOete 2 
a 


7 2a »>¢ 
aa ai pe nok? ns" 4 103 
_ -8Oet Bied. 103 . 


fev sed a 65 


137 


intrinsic hypothesis. It assumes that the drift and the 
variogram remain constant over the study area. 


If this hypothesis does not hold, then other more 
complicated methods in geostatistics must be used. In fact, 
the data in this study follow the slightly stricter assump- 
tion of weak, or second order stationarity, because the 
variogram has a sill and because VAR[Z(x)] is defined. 
However, only the quasi-stationary assumption was tested 
Since the data is only known to be stationary within a 
neighbourhood defined by h < H. H is the maximum distance 
used in the kriging equations, about 1000 m in this study. 
However, for simpicity the basic assumption used will be the 
intrinsic hypothesis. 


WZ, DRIFT 


The drift is shown as a graph of the average difference 
between all pairs of data separated by a vector plotted 
against the size of the vector. As suggested by David 
(1977, p.267) the drift is estimated using: 


ah ee aia NGM eke ie ey) balers ds: 
Where: 
N (h) = the number of sample pairs separated by h. 
Zxe)) =a Sample at J6cation. x, . 
Z(x,+th) = a sample at location x,+h, which is removed from 


. Sample Z(x;,) by the vector h. 
Note that, like all basic formulas in geostatistics, the 
formulas for the drift come in at least two varieties: a 
definition and a estimation formula. Equation {1.4}, p.136 
defines the drift, but it must be translated into 
equation {1.6}, p.137 to actually calculate an estimate of 
the quantity. For the intrinsic hypothesis to hold the 
drift must be zero, or: 


Mh) 0 forall n AL 


According to Clark (1979a, p.23), the condition usually used 
im practice is. 


Mineo sO te fOrera lnc {Ajit 


If the intrinsic hypothesis does not hold, then ordinary 
kriging will give a biased estimate and universal kriging or 
the method of intrinsic random functions must be used to 
obtain an unbiased estimate. Examples that have a large 
Ggitc are shown in eidgures tela The curve of thesdritt moves 
continously away from zero axis in either the positive or 
negative direction and with increasing slope. The vario- 
grams are unbound, which means they have no sill or steady 
horizontal portion for large h. The unbound portion is 
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Sulphur in sa, coal*-seam Leadwinmawe LeadszLoc mine 


Modified after David (1977, Figures 190-193, p.268) 


Figure 1.1 Examples of the drift and variogram from 


data with a large drift 
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parabolic, with the slope increasing as h increases. 


1.3 VARIOGRAM 


"Variogram' is used in this study to refer to the semi- 
variogram, which is defined by equation {1.5}, p.136, 
whereas, the variogram is defined by the same formula, less 
the factor 4. Most geostatistical authors use 'variogram', 
and a few even use both terms, when they mean 'semi-vario- 
gram'. Ambiguity is rarely a problem if care is taken. 
Assuming the variogram exists, a formula for estimating it 
1S: 


qth) ae SPN hy ele tey th) - 2(x,))2 Leg! 


Where the terms are as defined for equation {1.6}, p.137 and 
y™ estimates y. By convention in geostatistics the '*' is 
usually dropped from y*, so y refers both to the unavailable 


real variogram (y) and to the estimate of it (7%) 


Variograms usually slope up to the right and a common 
behaviour is to level off after a certain distance 
(Figure 1.2) called the range (a), while the height at which 
He tlevelseort vs the sili (Cc): 


1.3.1 Conventions Used to Illustrate Variograms 


Graphs of variograms in this study generally follow sev- 
eral conventions; one set for variograms of the coal 
thickness, another for theoretical variograms. 


When the variogram is of data from this study, the units 
are always the same. The horizontal axis (the lag or h) is 
in metres, while the vertical axis (the variance of the 
nett thickness) is in metres*. Standardized scales for the 
axes would have facilitated easy visual comparison of vario- 
grams between figures, unfortunately, the pit variogram 
would have been too small, obscuring any detail. All 
calculated points are plotted, except those using fewer than 
50 sample pairs. The plotted points are connected or the 
model is drawn out to L/2, beyond which there is no rela- 
tionship between the ideal and the estimated variogram.¢ L 
is the maximum distance between sample pairs. 


When the variogram is theoretical, such as an 
illustration of a model type, no data is used so the curve 
is drawn to the limits of the graph. Direct comparison of 


+ Two variograms also contain a gross thickness variogram. 
+ See Appendix 1.3.3 for a discussion of the variogram 
estimation errors. 
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Modified after Clark (1979a, Figure 1.3, p.7) 


Figure 1.2 A transition variogram, spherical model 
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different variograms is done by scaling the axes, the hori- 
zontal with the range (a) the vertical with the sill (C). 


1.3.2 Anisotropic Variograms 


Anisotropy in a variogram is a substantial difference 
in the shape of the variogram depending on the direction of 
the vector h. In sedimentary variograms it is almost always 
present between the vertical and horizontal directions. 

This is because there is more variability vertically through 
the section than along it. It is also often present between 
Girections parallel and perpendicular to the depositional 
Strike, because of greater variability down the depositional 
dip. In variograms there can be two types of anisotropy 
(Figure 1.3): geometric anisotropy, sometimes called affine 
anisotropy; and zonal anisotropy, sometimes called 
Stratified anisotropy. Zonal anisotropy is the general 
case, while geometric anisotropy is a particular case of it. 


Geometric anisotropy is present when the variograms for 
different directions have the same sill but different 
ranges. A simple linear transformation or scale change of 
the coordinates of the data can correct for this, following 
which kriging can be done using an isotropic variogram. The 
linear transformation of the coordinates has the effect of 
Stretching and/or shrinking the different ranges to the 
average range. Variograms parallel to the depositional 
Strike often show a larger range than the perpendicular 
variograms. 


Zonal anisotropy is present when the variograms for the 
different directions have a different sill for bound models 
and usually also a different range or possibly a different 
model altogether. This cannot be corrected by a simple 
linear transformation of the data coordinates. Instead it 
is corrected by subdividing the data into groups that are 
each isotropic or geometrically anisotropic. Variograms 
parallel and perpendicular to the bedding often show zonal 
anisotropy. 


1.3.3 Estimation Errors for the Variogram 


The variogram as estimated by the data is rarely the 
variogram as required by the theory. Estimating a variogram 
produces a smoothed version, that must uSually be corrected 
to remove the effects of the smoothing. This smoothing has 
sources: firstly, each sample almost always has a finite 
volume, resulting in a regularized variogram. Secondly, the 
data is grouped into distance and angle classes, analogous 
to, and with the same effect as the construction of a 
histogram. 
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y(h)/C 
y(h2) = y(a,h,/az) 


; - - ee _ 7 — = 
GEOMETRIC (AFFINE) ANISOTROPY 


y(h) 


ZONAL (STRATIFIED) ANISOTROPY 


ModifiedhaftereOlea.(1975, Figures 2.2.2%and 2.223) p.22) 


Figure 1.3 Types of variogram anisotropy 
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1.3.3.1 Regularization of the variogram 


A regularized variogram results when the data has a 
finite support, rather than infinitesimal or point support. 
This means that each data value is the average of the varia- 
ble over some small but finite volume or area, for instance 
the average grade of a piece of core. Support also includes 
the shape and orientation of the sample but this does not 
apply to thickness measurements, which are made at very 
close to an ideal point. This smoothing will tend to lower 
the sill and increase the range of the raw variogram with 
respect to the ideal point variogram (Figure 1.4). The 
effect of the regularization is measured by the ratio l/a, l 
is the size of the support, and ais the range. The larger 
l/a is, the more the smoothing, but unless l/a > 0.1 the 
smoothing will not uSually be apparent. 


1.3.3.2 Grouping of the data 


Data is often irregularly distributed on a plane, so it 
has to be grouped into classes in order to calculate a 
variogram. There are two ways that data must be grouped 
(Figure 1.5): lag (Ah) determines the distance classes, and 
psi (Ww) determines the angular classes. Optionally, data 
can also be grouped into directional classes, determined by 
p. 


The lag (Ah) is the linear regularization. It is the 
average separation between the points on the graph of a 
variogram. It is also the width of the distance over which 
h is averaged for one point on the experimental variogram. 
A larger value of the lag will smooth the variogram and 
detail will be lost, particularly near the origin. A 
smaller value will produce fewer data pairs for each point 
on the variogram, making it unstable and unreliable. 


Psi (W) is the angular regularization. It is the half 
angular width of the directions of h within which data pairs 
are grouped to calculate a point on the variogram, and is 
the angular equivalent to the lag. If it is small, then few 
pairs can be found for the first points on the variogram, 
and these points on the variogram will not be stable. If it 
is large, then an excessive number of pairs will be found 
for the more distant points on the variogram, which will 
increase the computing costs. Large values will also tend 
to smooth any geometric anisotropies. 


If y = 90°, then all pairs are used irrespective of the 
direction of h, and the average variogram will be estimated. 
This variogram is called the one dimensional variogram, even 
ff the dava is distributed in 2 or 3 dimensional space. 
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YC Cal 


POINT VARIOGRAM 


REGULARI ZED 
VARIOGRAMS 


0 1 h/a 


Modified after David (1977, Figure 100, p.131) 


Figure 1.4 Effect of regularization on a spherical 


variogram 


: F 
' ie 
~ 
4 ‘ 7 
2 
7 1 
“4 
“Bs 
Pe ha * 
MA, tx * 4 P - 
soe re 7 ae a a 
(i dworteton pow wo! ow ay be te - oS SES SS Oy 
- ‘or a a or ¢ 
¥ } la ba ~ 
'- o™M - : 
rote we : ecm as 
GST THAIS 
rAd Fa Bet” Bets 9 fd <> 
oMASTROT HAM 
, en ae ae 


‘ 


ee ee ee eae a ee eT 


re) { 


a 
ie 


i) Bivad wets Beliibow ; 


Pn, 
i 
~ 
c 
c 

4 

1 

c 

7) 

1 aed 
ou 
- 
al 
it? 


rns i 


ja 6 80 Soliesizsivges to 259393 


145 


To calculate the variogram two points are needed. 
The’ fyrst- point Sela Ae. 


The second point is any that falls inside box B. 


Meod@ined after Flint (1078 Figure 20; p.59) 


Fv guresl.(5 Grouping of the data for calculating a 


variogram 
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The direction of the angular classes (¢) is used to 
search for any anisotropy. At least two directions should 
be chosen, perpendicular and parallel to the suspected 
anisotropy in the data. Good practice also suggests examin- 
ing two more directions rotated 45° from the first two, in 
case the original ¢'s happen to bisect the anisotropic axes. 


1.3.4 Behaviour of the Variogram Near the Origin 


The variogram can have four different behaviours at the 
origin (Figure 1.6). “In decreasing@order of probability of 
occurance, these are: 

1. linear behaviour, 

2. (‘nugget effect, 

3. parabolic behaviour, 
4, pure nugget effect. 


Linear behaviour is the commonest property at the 
Origin. The variogram is continuous, but not differentiable 
at the origin. This means the variogram intersects the 
Origin with a non-zero slope: 


y (h) = 0 “pas hi. .0 1104 
dy(h)/dh > 0: #4 as th 3.0 via cra pal 
Where: 


dy(h)/dh = the first derivative of y(h) with respect toh. 
In this case it is the slope of the variogram in 
CHE™ATreEeCTIOn- OL it. 


A nugget effectf is present when the variogram does not 
pass through the origin as on Figure 1.2, p.140. 
Theoretically, the variogram always passes through the 
Origin, since the covariance between a sample and itself is 
zero. In practice, the experimental variogram often does 
not. This paradox is resolved by stating that the value of 
the variogram at zero is not the value of the variogram as 
it approaches zero: 


y(h) # 7(0) =: as h > 0 a Beard 


The value of the variogram as h » 0 is called the nugget 
effect, (Co): wittis defined by: 


ER OGRE Oe ENS Bohan ie) elraahs a; 
The nugget effect is the random variance, uncorrelated from 
+ The term 'nugget effect' originated along with mining geo- 


Statistics in the goldfields of South Africa where nuggets 
of gold produce this effect on the variogram. 
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Modified after Matheron (1971, p.58) 
Figure 1.6 Types of variogram behaviour near 
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sample to sample, while the rest of the variance is correla- 
ted. When the nugget effect is present the slope of the 
variogram near the variance or Y axis is non-zero, as in the 
linear case in equation {1.11}, p.146. 


The nugget effect adds a constant to the model: 
With) = "Co. & mode. flare | 


According to Rendu (1978, p.20), it can have two causes. 

The first is due to sampling or meaSurement errors, some- 
times called the human nugget effect. The second is 
irregularities in the variable with a range shorter than the 
lag or the sample spacing. The relative nugget effect (C,) 
is the ratio of the nugget effect to the sill: 


Cy, = Co/C a at 
Where: 
Gla ate Overall Sli. main bale uCUlal +. € =" CotCqt Coe bs ee FC. 


Gaecanuvary £rom.0,).Or no nugget effect sto li for. apure 
nugget effect, which is discussed below. 


Parabolic behaviour near the origin is rare in geostati- 
Stics since it is the property of a very regular 
regionalized variable. Sedimentary thicknesses and some 
meteorological data, according to Gandin (1965, pp.49 and 
54), can often be this regular. The variogram is continuous 
at thesorigin (equation tb 10), bp. 146) ssSimitar to the 
linear behaviour, but unlike the linear behaviour it has a 
zero slope at the origin: 


dy(h) /dh y= "0 7 as th 4°60 Cea: 


A pUre wmugget ectece pisatne limiting ymodeletwlth. Cu =" 1) for 
a transition variogram when the sample spacing is larger 
than the longest range. Of the several hundred deposits 
that have been analysed with geostatistics, only one showed 
a pure nugget effect. With this model classical statistics 
is aS good as geostatistics since it gives an unbiased esti- 
mate with the same minimum estimation variance. 


1.3.5 Confidence Intervals on the Fundamental Variogram 


The experimental point variogram is not the true vario- 
gram that is required in the kriging calculations. The true 
variogram is defined uSing every point in an area of 
infinite extent. In contrast, the experimental variogram 
(y*) estimates the true variogram (vy) by using (1) a sample 
of all the potential data from (2) an area of limited 
extent. These two finite approximations produce the 
estimation error and fluctuation error respectively, which 
have been explicitly determined only for the case where the 
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model is linear and the samples are normally distributed. 

In this study the samples are close to normally distributed, 
but the model is spherical rather than linear. The results, 
however, should give a rough indication for the present 
data. 


The estimation error 1S due to a finite rather than an 
infinite number of samples being used to calculate the 
variogram. The error is usually small, particularly when 
the number of samples is large as in this case (n = 200 to 
400); consequently, it can be safely neglected. 


The fluctuation error is due to the variogram being 
calculated over a finite area rather than an infinite area. 
This error can be very large. In the linear/normal case 
outlined above, the error is small for small distances. 
However, when h becomes larger than about half the total 
distance over which the variogram has been calculated (L), 
then there is almost no relationship between the true and 
the experimental variogram (David, 1977, p.116) 


Nothing can be done to improve the fluctuation error 
except to ignore all points with h > L/2 when modelling the 
variogram. More conservative geostatisticians use only the 
datawitheh < L/45 for instance: Clark .(1979b,. 0.92). 


1.3.6 Variogram Models 


The kriging system of equations needs the value of the 
variogram calculated between many pairs of points. A mathe- 
matical expression of the variogram is needed to do this, 
rather than the usually ragged set of points from the 
experimentally determined variogram. Modelling is the art 
of fitting a mathematically defined curve to the data, much 
as a histogram is fitted by a probability density function. 
Unfortunately, there are no automatic fitting methods. 

Davis and Borgman (1979) give a test for the goodness of the 
fit of a one dimensional variogram, but it is worked out in 
detail for a test of only one point at a time. More useful 
cases were not described. Another method is cross 
validation, which is choosing the model that produces the 
lowest eStimation variance. Iterative visual fitting, the 
commonest method, was used in this study. 


Not all functions may be used for a variogram model: 
they must be conditionally positive definite, which has two 
consequences. Firstly, any estimation variance is greater 
than or equal to zero, and secondly, the variogram must have 
less than a parabolic increase for large h, or: 


Vim ythy 7h) =0. + aas hese Taller hy et: 


All the models discussed in this study are conditionally 
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positive definite. 


Variogram models fall into two broad classes: transition 
models, and non-transition models. Transition models are 
bound variograms and therefore have a sill, while non- 
transition models are unbound and lack a sill. Preliminary 
inspection of variograms from the Boundary Dam Mine show 
that they are of the transition type so non-transition 
models will not be discussed further. 


Transition models rise to a value, then stay steady at 
this value (Figure 1.2, p.140). The distance from the 
origin to the steady value is the range (a). The height at 
which it remains steady is the sill (C), which is 
theoretically equal to the population variance (07) and 
usually taken to be equal to the sample variance (s?). In 
other words: 


y(o) = s? eer ae 


If the sill and variance are not equal this can be a result 
of two causes: either the presence of drift in the data, or 
a hole effect in the variogram. 


If there is a drift, then the variogram will 
continuously rise above the sill, approximately 
berabolically CRigure tl, p.138). 


A hole effect can be produced when rich/thick zones 
alternate with poor/thin zones. If the alternation is 
approximately regular, then this periodicity will show as 
oscillations in the experimental variogram above and below 
the population variance (Figure 1.7). 


There are three types of commonly-used transition model 
(Figure 1.8): 
1. spherical model, 
2. exponential model, 
3. gaussian model. 


The spherical modelt is by far the most commonly used 
model in geostatistics. Its formula is: 


yinere Ci shyca - h-/2aewee torn h sa 
eG ehOneh 2 1a eet 
Where: 
Eve athe. si 1); 
h = the sample separation. 


+ The spherical model is occasionally called the Matheron 
model in honour of Georges Matheron, the father of 
theoretical geostatistics. 
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The hole effect is due to alternating 
rich and poor layers in a uranium mine. 
Modified after David (1977, Figure 89, p.111) 
Figure 1.7 An example of a variogram with a hole 
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Modified after Rendu (1978, Figure 4.1, p.18) 


Figure 1.8 Common types of transition model: 


Spherical, exponential, gaussian 
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a = the range. 

Its ubiquity is a result of several reasons. Traditionally, 
everybody has used it, so people continue doing so. 
Considerable theoretical work has been done on it, and as a 
result many of the problems have been solved. It usually 
seems to fit the data better than the other two. The 
exponential model rises too slowly to the sill and the 
gQausSian model is usually too flat or parabolic near the 
Origin. David (1977, p.106) explains the origin of the term 
"spherical'. "A process having a spherical variogram can be 
generated as follows: first take a random (Poisson) distri- 
bution of masses in space. Then assign to each point of 
Space a grade given by the number of masses within a radius 
a of the point (in other words take the number of masses 
within a 'sSphere')." 


The exponential model is more commonly used in 
biometrics. In this science, data tend to be more regular, 
which causes the variogram to rise more slowly to the sill. 
The formula is: 


y¥(h) = CH1 = exp{-h/al} 20% 


The gaussian model is rarely used in geostatistics 
except where the parabolic behaviour at the origin is impor- 
tant. Whe formula is: 


y¥(he=) Ci | oauexpl ahs Aas (elpee2aiah 


Occasionally it has been fit to the variogram of sedimentary 
thicknesses (Clark 1979b, p.94), since they are usually 
highly continuous and often parabolic near the origin. It 
has also been used in meteorology for sea level pressure 
(Gandin, 1965, p.54). 


Hole effect models are another class of model, either 
transition or non-transition. They are usually 
trigonometric functions with a decaying amplitude. A common 
and simple transition model is the sine model (Figure 1.9) 
which is an inverted sinc function: 


y(h) = C{t - sin[lzah]/zah} alee 22e 


Other hole effect functions are described by David (1977, 
p.111) and by Journel and Huijbregts (1978, pp.252 and 169). 


A nested model is a linear combination of several simple 
models with different parameters (Figure 1.10): usually two 
models, very rarely three, almost never four. It has the 
formula: 


y¥(h) = y,(h) + yoCh) + yeChjyt «++ + >, (h) ae 3h 
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Figure 1.9 The sine model, an example of a transition 


hole effect model 
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Modified after Rendu (1978, Figure 4.6, p.20) 


Figure 1.10 Two nested models 
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Where: 

y,(h) = any valid variogram model. 

Any linear combination of conditionally positive definite 
functions is also a conditionally positive definite func- 
tion, so it also is a valid variogram model. The sills of 
each model are the combination weights. 


1.3.7 Robustness of the Model 


Geostatistical results are uSually very robust with 
respect to the variogram model, both type and parameters 
(Journel and Huijbregts, 1978, p.167). Most models have 
Similar graphs even though they may have different func- 
tions. For instance Journel and Huijbregts (1978, p.246) 
show some data that is fitted equally well witha 
logarithmic model, which has no sill, and two nested 
Spherical models. They also show on page 234 how close an 
exponential and two nested spherical models can be 
(Figure 1.11), much closer than the scatter in most data. 


de de ae Two nested spherical models 


One exponential model 


(hy /C 


ha, 


From Journel and Huijbregts (1978, Figure III.30, p.234) 


Figure 1.11 Comparison of an exponential and two 


nested spherical models 
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APPENDIX 2 
PRELIMINARY DATA ANALYSIS 


2.1 CLEANING THE DATA 


Once the data was in machine readable form mistakes 
were checked for and corrected. Each data point, pit and 
Grill hole, has three main variables: ID, location and 
thickness, all of which were checked. ID's not in 
sequential order after sorting indicated missing data which 
was redigitized. 


2.1.1 Location of the Samples 


Incorrect locations are due to either mistakes or 
errors. Mistakes, points digitized at a wrong location and 
missing or duplicate points, were corrected in two ways. 
Firstly, the posted data was overlaid on the original blue- 
line mine plan then missing or incorrectly located data 
redigitized. Secondly, duplicate points were found by 
visually scanning a sorted list of coordinates then 
eliminated. 


Errors in location have several sources: 
1. Surveying errors, 
2. slight scale changes introduced in the blue-line and 
xerox copying, or paper stretch, 
a) Glarei zing errors, Which are Of three kinds: 
Ae INCOLLeCt TequStrarvren, or lOCatlam OL, the map son, the 
table, 
B. not placing the mouse directly over the point, 
C. the tablet not recording the precise location of the 
mouse. 


Nothing can be done to check the surveying errors; 
however, they are likely to be at least an order of 
magnitude less than the +1 m error due to digitizing and so 
can be ignored. Scale change errors are automatically 
compensated for when registering the maps before digitizing. 


The digitizing errors could not be totally ignored or 
eliminated; however, they could be measured. Several 
points, whose true locations were known, were digitized 100 
times, to give an average +1o error of approximately 1 m. 
This is small enough to be ignored for the variogram and 
kriging but not for the error of the area, which is examined 
in more detail on page 92. 


2.1.2 Thicknesses of the Seams and Partings 
Thicknesses, which may have been recorded incorrectly 
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by the mine personnel, were checked by searching for extreme 
values. Only one fell beyond 30 from the mean, but it was 
retained since it was in an area of channelling and rapid 
thickness variations. 


Thicknesses that were incorrectly keyed into the 
computer were checked for by hand. Six readings in pit 2A 
were found to have a peculiar stratigraphy, where a 
thickness was recorded for the bottom coal and the parting 
but not for the top coal. These were discarded since it was 
not possible to determine the net coal thickness. 


2.2 METHOD TO ESTIMATE NET THICKNESS 


To calculate volume and tonnage it is necessary to 
estimate the average net thickness of the seam over a given 
area or domain. 


For the pit data this can be done in three ways: 

1. estimate average net thickness directly, 

Subtract the average parting thickness over the whole 

area from the average gross thickness, 

3. add the average thicknesses of the bottom and top coal. 
Then calculate a weighted average of this and the aver- 
age thickness where there is no parting. 

Method 1 was chosen for the following reasons. 

1. Calculating directly 1s simpler because less work is 
involved and it is less costly as the kriging program is 
expensive to run. 

2. The net thickness variograms are simpler and better 
behaved than the gross thickness variograms, which have 
a strong adrift (Figure 2.1 andsriguue 2.27) 4 As a 
result the net variograms can be modelled more easily 
and more accurately, which will tend to produce a better 
estimate. 

A. Method 2 requires use of the parting variogram 
(Figure 2.3), which should be modelled by at least 
one and possibly two nested hole models. 

B. Method 3 requires the variograms for the bottom and 
top coal and for the seam with no parting 
(Figure 2.4), all of which are more ragged than the 
net thickness variogram (Figure 2.1, p.160 and 
Figure 2.2, p.161) and would probably require hole 
effect models. 


For the drill hole data, method 1 or 2 above could be 
used, while method 3 would be impossible because the study 
area cannot be subdivided into simple contiguous regions of 
consistent stratigraphy (Figure 2.8, p.50). Method 2 would 
be very poor practice since the parting variogram would 


+t See page 139 for a description conventions used to 
illustrate variograms in this study. 
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Ferlect Severe parti OemanOcm d..Si1g le phenomenon. The 
variogram would also have a larger sill, which increases the 
variance of the net thickness estimate. For these reasons 
plus the time and expense involved the drill hole net 
thickness was also estimated by method 1. 
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APPENDIX 3 
SELECTING HOMOGENOUS REGIONS FOR ANALYSIS 


Geostatistical analysis must be done on regions that are 
homogeneous both geologically and statistically. Although 
the geological evidence should have primacy over the 
Statistical, in this study there are contradictions which 
must be resolved. 


3.1 SELECTION CRITERIA 


Homogeneous regions were selected using three criteria: 
the geology must be the same, the data must be from one 
population and the intrinsic hypothesist must hold over the 
region. 


The geology must be uniform for the results to be mean- 
ingful. For instance, it is not very useful to estimate the 
average thickness of a zone that includes both coal and a 
Sand channel. 


The region must also be from one statistical population 
if any statistical inferences are to be made. For instance 
it is difficult to make any inferences about a region that 
both includes and excludes a rider seam. 


Finally, the geostatistical analysis assumes that over 
the region of interest the intrinsic hypothesis is true. 
This hypothesis assumes that there is no drift in the data 
and that the variogram is constant over the region under 
Study. 


The pit and drill hole data are not homogeneous, so were 
considered separately for the following reasons. 

1. Both sets of data have different sources: the highwall 
and core or electric logs. 

2. Both have different geology. Even though the measure- 
ments were made on the same zone, they were measuring 
two different things: the mined zone for the pit data, 
and the geological zone for the drill holes. 

3. Both sets of data are from significantlyt different 
populations (Table 3.1). 

4, The variograms are different (Figure 3.1). 


For each set of data, pit and drill hole, the regions of 
interest were selected uSing these steps: 
1. choose a region of uniform geology, 


+t See page 136 for a description of the intrinsic 
hypothesis. 

+ Generally whenever 'significant' is used in this study a 
Significance test was done, in which case the confidence 
level will be reported ina footnote thus: a = 0.05 
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Table 3.1 Comparison of the summary statistics of the 


thickness for the pits and the drill holes 


Data n Mean Variance ene F 
thickness (m?) 
(m) 


All pits 634 3.64 O07 
Loa Sie Onl 
Hira) 336 B29 D276 
holes 
Gravicalevalues. (ae=— 0.705) 1.96 (Aci aes 


+ The t test for equal means assumes that the variances are 
egual, which 1S not the case here. However, the variances 
are not grossly different and the test statistic is so large 
(13.9) that the conclusion is probably cotrect. In eithe€® 
case the variances are significantly different, which is a 
sufficient condition for the populations to be different. 


+ The '+' on this and subsequent tables indicates that the 
value has been very crudely interpolated from the 
statistical tables. 
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2. check that the data is from one population using 
classical statistics, 

3. check that the intrinsic hypothesis is valid over the 
region. 


An important assumption of both the t test and the F 
test is violated in this study for the thickness, almost 
always the variable to be tested. Each assumes that the 
data samples are random and independent. The data has not 
been selected at random, but on regular grids (Figure 2.1, 
p.39 and Figure 2.2, p.40). Neither is the thickness 
independent, it 1S regionalized with a variogram that cannot 
be modelled with a pure nugget effect. The test results are 
not strictly correct, but are calculated as a rough 
LNGi Cat OTs 


3.2 SELECTING HOMOGENEOUS PIT REGIONS 
3.2.1 Geology of the Pit Regions 


The only geological data available from the pits is the 
thickness of the coal and the partings. The thicknesses do 
not “changé drastically (Figure 1013, pe28”and Figure. 116, 
p.32). The study area was subdivided on the basis of pres- 
ence/absence of the D°E and E*F partings (Figure 1.15, 
p.31). This resulted in three areas: 

Ve “a Trorth area@split by tChe-E-F parting; 

2. a middle area of solid coal with no partings, 

3. a south area split by the D-°E parting. 

Not enough points (3 to 7) are present on the F-G parting to 
separate it out from the north area. 


The mine surveyors only meaSure a parting when it is 
selectively removed in mining. It is possible that the 
partings extend into the middle area but were too thin to 
separate out. Figure 2.8, p.50 shows many holes in the 
middle area that have a parting logged in them. The 
Ppartings appear to terminate more by wedging out, as 
Prgure 1 ?t6, ">.32 suggests, rather than by shaling out. 


If this is the case, then the middle area will contain a 
mixture of two or three populations: 
to" north area population, 
2. a possible overlap or gap between the two wedge-edges of 
the partings from the north and south areas, 
3. south area population. 
The statistics of the middle area should be a mixture of the 
Statistics of the north and south areas. 


3.2.2 Classical Statistics of the Pit Regions 


The cumulative normal probability curve of the middle 
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area has several properties of mixed population. It plots 
midway between the curves for the north and south areas 
(Figure 3.2), with its ends almost parallel to the curve 
that dominates the extremity. In this case the thin end is 
Sub-pata lel to themnonensarea curve. Finally, itis 
concave towards the higher-variance curve (north area) 


Mixing 1s confirmed by calculating what the summary 
Statistics should be, assuming mixing with no overlap or 
gap. The population mean and variance of the middle area 
should be functions of the population mean and variance of 
the north and south areas. This was checked by using the 
Sample statistics. 


The sample mean of the middle area (Z) should be a 
Simple weighted average of the means for the north and south 
areas, Or: 


ie oy A GaN i ad a A 0s am ba 


Where: 
Zn, = Sample mean of the northern area. 
s = Sample mean of the southern area. 

n = Np/(n, + n,), the proportion of all middle area samples 
that belong to the north area population. 


fo=) mh, /in, + M,.), the=proportion of ‘All middle arearsamples 
that belong to the south area population. 

Nn, = Number of samples from the northern population in the > 
middle area. 

n, = number of samples from the southern population in the 


middle area. 


The variance of the middle area (S?*) should be a 
weighted average of the north and south variances corrected 
for the spread between the means, or: 


SParknss 4 kiss + kk Cayo) 25) Usiee 


Where: 

kn, Ks, 2n, 2; have the same meaning as in the formula for 
the mean. 

s? = sample variance of the north area. 

s? = sample variance of the south area. 

The above two formulas are derived in Appendix 4. 


The results of the calculations are shown in Table 3.2 
on the line 'North+South mixed'. The calculations assume 
there are two populations in the middle area that meet 
midway between the recorded wedge-edges of the partings, 
which is between pit 2A and pits 1 + 3. Statistics measured 
for the middle area are the sameft as the statistics 
estimated assuming a mixed population. 
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the north, middle and south areas of pit 


data 
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Comparison of the 
thickness for the 


areas of pit data 
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summary statistics of the 
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Variance 
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O72 0756 
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Standard 
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Test statistics from comparing the middle area and the 


North+South mixed sample 
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To further test the assumption of a mixed population the 
middle area was divided into two subareas: a north-mid 
Subarea and a south-mid subarea (Figure 3.3). The sample 
statistics for these subareas were compared to each other 
and to those for the north and south areas (Table 3.3). The 
two subareas of the middle area are from significantlyf dif- 
ferent populations; but on the other hand, there is no sig- 
nificantt difference between the north area and the north- 
mid subarea or between the south area and the south-mid 
Subarea. 


3.2.3 Geostatistics of the Pit Regions 


The variograms for the middle subareas were also 
compared to each other and to the variograms for the north 
and south areas. The variogram for the north-mid subarea is 
more Similar to the variogram from the north area than to 
the one from the south area (Figure 3.4). The sills are 
almost the identical, though there is some difference in the 
range. The variogram from the north-mid subarea has a sill 
three times that from the variogram from the south area. 

The variogram from the south-mid subarea is most Similar to 
the variogram from the south area (Figure 3.5). The ranges 
are Similar and the sills are close. The variogram for the 
south-mid subarea has a different range and sill from the 
variogram for the north area. 


3.2.4 Homogeneous Pit Regions 


Since the middle area is composed of two mixed 
populations, it was subdivided between the pits and each 
portion combined with the north and south areas to give two 
homogeneous regions of pit data: a southern region, composed 
of pit 2A and a northern region, composed of pits 1 and 3. 
(Figure 3.6): 


Within these regions the intrinsic hypothesis holds. 
There is no drift (Figure 3.7), except possibly in the 
northern region. This drift does not begin to become appar- 
ent until about 1000 m, which is at, or beyond the limit of 
any estimating done in this study. The two subregions (with 
and without a parting) within each homogeneous pit region 
have approximately the same variogram (Figure 3.4, p.175 and 
Figure 3.5, p.176); therefore, the variograms are probably 
Stationary. 


Even though the homogeneous regions have mixed geology 
(the partings wedge-out part way across), the statistics 
Suggest that in fact the geology is uniform and that the 
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Figure 3.3 Location of the areas and subareas of the 


pit data 
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Table 3.3 Comparison of the summary statistics of the 


thickness for the areas and subareas of pit 


Gata 

Area n Mean Variance t F Fecr Re 
Or thickrmess (m7} (0.05) 
Subarea (m) 
North TZ Beeols Or07 36 

1285 Tes Giese 
North-mid 258 Sie (Syd O07 ir4 

Ss bot Zio 1 Sipe 
South-mid We & Aube OR O27 8 

5 a | Weed 16.52 
South 79 Sera. 0.0461 
Perr. 00.05)) 1.96 


+t See footnote f under Table 3.1, p.166. 
+ The variances are equal at a = 0.02. 
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partings probably extends across the regions, but are too 
thin to selectively mine and measure in the centre of the 
study area. 


These two homogeneous pit regions cannot be combined for 
these reasons. 

1. The geology is different Since there are different 
partings in each region (Figure 1.14, p.30; Figure 1.15, 
Destianad Figure 12065 p.3s2). 

2. The populations are different since the means and 
variances are significantlyt different (Table 3.4). 

3. The intrinsic hypothesis does not hold since the vario- 
grams are different (Figure 3.7, p.178). 


3.3 SELECTING HOMOGENEOUS DRILL HOLE REGIONS 


The drill hole data cannot be subdivided into separate 
homogeneous regions. The geology varies slightly, but the 
classical statistics and the geostatistics are the same over 
the whole area. 


There 1S a general tendency toward more splits in the 
north and south portions of the study area with fewer in the 
centre (Figure 2.8, p.50), but holes with different numbers 
of splits tend to alternate, so that it 1S not possible to 
separate areas underlain by a particular parting. 


The drill holes were subdivided into subsets based on 
the number of splits, irrespective of their location and 
Summary Statistics were compared (Table 3.5). In general, 
there is very little difference; the means tend to decrease 
as more partings are separated out but there is no trend in 
the variances. 


The statistics of each subset were compared against the 
population statistics estimated by the sample statistics of 
all the rest of the drill hole data without the subset in 
question (Table 3.6). This was repeated for each number of 
splits from one (with no partings) to five (with four 
partings). In general the comparison between the various 
Subsets was either good (both mean and variance fit) or fair 
(either mean or variance did not fit); there are no gross 
discrepancies. The subsets of holes with one and three 
splits of coal did not fit very well, but the difference is 
not large. It is reasonable therefore from the standpoint 
of classical statistics to combine the drill hole data even 
though they include different numbers of splits. 


If there are two populations of drill hole data similar 
to the pit data, then they should have about the same loca- 
tion as the pit populations. The drill holes were divided 
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Table 3.4 Summary statistics of the thickness for the 


homogeneous regions of pit data 


Region n Mean Variance e F 
thiekness (m? ) 
(m) 
North 430 S| Aye, O207'59 
Sieben Zits 
South 201 SUIS 0.0356 
Grrtacal values (Cac= 0A05) 1.96 1.3+ 


t+ See footnote T under Table 3.1., p.166. 


Table 3.5 


Number 
of splits 
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Summary statistics of the thickness for the 


drill holes, subdivided by the number of 


splits 

n Mean Variance 
thickness (m? ) 
(m) 

80 BS 0.298 

143 BG Os 07 

94 Ay, 0.160 

v7 SOS, (he O33 


2 300 OV 00577 


ak monet BY 
Pict ("m3 
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Comparison of the summary statistics of the 


thickness for the drill holes, subdivided by 


the number of splits 


Table 3.6 

Number n Mean 
of (m) 
splits 

| 80 3.48 
Alia. ~ 256 ees 
i 143 Sao 
Al L=2 193 34:29 
3 94 Chea 
All-3 242 SESS 
4 1 33 
Aelia 319 3°30 
5 Zz SOO 
Ae 334 3229 


erie (0.05) 


Topol al 


Solty 


Variance t 
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into two areas (Figure 3.8), at approximately the same 
division used for the pit data (Figure 3.6, p.177). The two 
areas are statistically homogeneous (Table 3.7) since the 
means and variances are the same.t The variograms are 
essentially the same with no drift (Figure 3.9); 
consequently, the intrinsic hypothesis is considered valid. 


It may seem odd that the same phenomenon (the Estevan 
Coal Zone net thickness) measured one way (by drill holes) 
has one thickness population, but measured another way (in 
the pits) has two populations. This difference has two pos- 
Sible reasons. Firstly, the two phenomena measured are not 
really identical. The drill holes measure the whole zone as 
logged by a geologist, whereas the pit measurements describe 
only the mineable portion of the zone as seen by the 
excavation equipment operators. Secondly, the variance of 
the drill hole data is large enough to include two 
populations. At the 0.05 level, a difference in the means 
of 0.13 m between a separate north and south drill hole 
population would be barely significant. This is very close 
to the difference (0.18 m) between the north and south 
regions of the pit data (Table 3.4, p.180). 
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Figure 3.8 Location of the north and south areas of 


the drill hole data 
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Table 3.7 Comparison of the summary statistics of the 


thickness for the drill hole areas 


Area n Mean Variance t F 
thickness (m2) 
(m) 
North 169 Sige 2 0e3Ts 
0.92 Tks 
South 167 Sake 0.369 


Critica values (a = 0205) 1.96 1.32 
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APPENDIX 4 
DERIVATION OF THE FORMULAS FOR A MIXED MEAN AND VARIANCE 


It 1S assumed that the pit data from the middle area 
with no parting is made up of two populations: one from the 
northern and the other from the southern portion of the 
middle area. 


4.1 MIXED MEAN 


The overall sample mean of the data from the middle 
area is defined as: 


Z = (1/n) Lz 


For simplicity Z is unsubscripted, but it should be taken to 
mean the summation over all values of the subscripted varia- 
bles. Separating the data from each population, this 
becomes: 


Deen) oe Mle oe fan 43 


The sample mean for the data from the north portion of the 
area only is defined as: 


ri ets MEVAG A IB rae 
This can be rearranged to: 
Deas = Rna2e {4.2} 


Substituting this and the analogous equation for the south 
portion of the area into equation {4.1} above, it becomes: 


Cai (el Ce et e/a) pes 


(ayy Aas Ate abel Cg RTA gy tae aot 
Define: 
ky =onn/ (eet on, ) 


and since; 


therefore: 
Ene cA dae 2 {4,4} 


Substitute this and the southern analog into equation {4.3} 
above: 
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lea Kato Fe. 
Z, is the mean of the sample from the northern portion of 
the middle area. But since it was assumed that this sample 
is drawn from the same population as the north area and 
Since the sample sizes are large (n > 100), Z, is also the 
mean of the sample from the north area, Similarly for z,. 


4.2 MIXED VARIANCE 


The overall sample variance of the data from the middle 
area 1S defined as: 


Cease (1/7) ee eae ai) (az ey eat 


Since the sample is large (n > 100) the biased form of the 
equation is as good as the unbiased form; therefore: 


t/n = i/(n—1) 


Separating out the data from each population as was done for 
the mean, this becomes: 


Se WE NG Age soBpam AS Pla cee tal ae ipa rare gp hyar, \be i: {4.5} 


The sample variance for the data just from the north portion 
of the middle area is defined as: 


ears al A ale ata tye ead «ley nbs 1 2 AN ae 
This can be rearranged to: 
re acain=e Gay Apes) (8am ae Baa MeV 


Substituting this and the south analog into equation {4.5} 
above along with some expansion, equation {4.5} becomes: 


S2icaP ish nai Sct e DSi A Me cez aetes (elf) (2. 23) ) 
SA Gli dati he 2batay hares uP Dea Sp An snot (Op ane) male, 


Substituting equation {4.2} and equation {4.4} above into 
this with some further rearranging, it becomes: 


SiasekSis ke Sc ate Gmiee/ mim eine eee Ol) G/N ea) ae 
EAE Grey Aon) bree Geta OPA NS ee Aho ras vate Meet GO pare Mate 
Collecting similar terms this becomes: 
Stee es) ia 9k) Sc ot a Neen ene Alien) eo 
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Substituting equation {4.4} above into this, it becomes: 


eke ea Soe eee Reg ers roe {4.6} 


Recall that: 


eee er rai nat ia) 


Keke GIN Ce. I hv w=. Ki, 


Substituting this and the south analog into equation {4.6} 
above, it becomes: 


Sa noe ek cen ee PO oe o26 ee eK ee 
This can be rearranged to: 

Soyo SHS LS hl See aimee ae 
s? is the variance of the sample from the northern portion 
of the middle area. But Since this sample is drawn from the 
Same population as the north area and since the sample sizes 


are large (n > 100), s?# is also the variance of the sample 
from the north area, Similarly for s?. 
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APPENDIX 5 
CALCULATING THE EXPERIMENTAL VARIOGRAM 


The raw variogram is calculated directly from the data, 
then the experimental variogram is calculated from the raw 
variogram by correcting for the smoothing. 


5.1 CALCULATING THE RAW VARIOGRAM 


The raw variogram is calculated uSing a program written 
by Flint (1978), which was modified from one written by 
David (1977, p.149). The program takes as input three 
parameters: lag (Ah), phi (¢) and psi (Ww), all illustrated 
on page 141. The program breaks all the pairs of points 
into classes depending on the values of the lag, ¢ and y, 
then calculates an average value of the variogram for each 
class (using equation {1.9}, p.139), which becomes one point 
on the raw variogram. An example of some output from the 
variogram calculating program is listed in Appendix 12. 


Not all points on the raw variogram are usable. If 
there are less than 30 to 50 pairs of samples used to 
calculate a point, then the value of the variogram is 
unstable (Journel and Huijbregts, 1978, p.194). Points with 
less than 50 pairs in the calculation were not plotted, for- 
tunately the number of pairs was typically in the 1000's, 
except near the ends. Since the lag was chosen to be just 
greater than the nominal sample spacing, no points were lost 
at the beginning of the variogram. The only points 
discarded were at the far end when h > L/2 and therefore not 
useful anyway. 


Various values of the parameters Ah, W and ¢ were used 
to calculate the raw variograms in order to find the best 
numbers to use. 


The lag (Ah) of the pit data was varied from 320 to 20 m 
in geometric increments. As the lag became shorter, the 
variogram became more ragged and unstable (Figure 5.1 and 
Figure 5.2), particularly when it approached the minimum 
Sample spacing of 30 m in the pits and 100 m in the drill 
holes. The reason is that fewer pairs were used to estimate 
each point on the variogram. 


The most important part of a variogram is near the 
Origin, where a nugget effect may be present. This part of 
the variogram can only be estimated by using a small lag. 
The lag must also be small enough to provide several points 
to estimate the shortest range; Journel and Huijbregts 
(1978, p.211) suggest Ah < 1/3 or 1/4 of the short range; 
the ranges are 50 m for the pits and 150 m for the drill 
holes. This is a conflict between the requirements of a 
Stable variogram (large Ah) and an accurate determination of 
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Co (small Ah). The rest of the variograms in this study 
were eStimated using a compromise value of 40 m for the pit 
data and 100 m for the drillholes. This did not produce a 
variogram which was too ragged and it showed the behaviour 
near the origin adequately. 


The usual choice for the size of the angular classes (Ww) 
f50c0 4 but, ts value vs not toc critical’, “Againg data from 
pit 3 and the drill holes is used to illustrate the effect 
Ciay 1 Figure Sosrancerroure, 5.4. With v= 10° and 80°, 
there is little change and therefore, the standard wy = 30° 
was used for all variograms in this study. 


The last calculation parameter, ¢, is the direction of 
the angular grouping and is used to any anisotropic struc” 
ture to the data. Each anisotropic variogram was calculated 
along four directions, 0 -, 45°, 90°). -and?1S5> tombe 
Suspected anisotropic axes. The suspected axes were 
parallel to the benches (Figure 5.5) and parallel to the 
major paleochannel (Figure 3.7, p.68 and Figure 3.8, p.69). 


5.2 CORRECTION OF THE RAW VARIOGRAM TO FORM THE EXPERIMENTAL 
VARIOGRAM 


Smoothing of the variogram is the change in its shape 
due to the samples not being at mathematical points and to 
the grouping of the samples when calculating it. This 
change must be reversed to caclulate the experimental vario- 
gram from the raw variogram. 


All thickness meaSurements in this study were taken 
almost at a Single point. Each pit thickness was probably 
an average over the area (1 to 10 cm wide)f across which the 
mine Surveyor visually estimates the limit of the seam; in 
this: case l/at=1s at most 1051 m/ 100 mer 05001. the drili 
hole thicknesses were averaged over the diameter of the core 
(0.1m) or the depth of investigation into the rock by the 
geophysical tool (about 1m), which gives l1/a = 0.001 and 
0.01 respectively. In each case 1/a is so small that the 
raw variogram and the variogram corrected to a point support 
are indistinguishable; thus, the variograms were not 
corrected for finite support. The effect of different l/a 
ratios is illustrated on page 143. 


Grouping smooths the variogram due to non-zero Ah and y. 
Experience has shown (David, 1977, p.146) that the 
corrections are minor if these two quantities are small. A 
small lag equals the average distance between samples and a 
small Ww is less than 45°. These constraintsS™»are 'satisfied 


+ All the calculations in this paragraph are order of 
magnitude only. 
+ A discussion of 1l1/a can be found on page 143. 
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by the values chosen for the parameters in this analysis. 
The sample separation is 30 m in the pits and 100 to 200 m 
for the drill holes, while the lags used were 40 m and 100 m 
respectively. The angular regularization (W) was always 
30°. The corrections can be computationally complex and the 
final results probably do not justify the extra effort. 
Since no corrections have to be made, the raw variograms are 
also the experimental point variograms (Figure 3.2, p.57 and 
Btqure -3..3,, DvoSun 


5.3 ANISOTROPY OF THE EXPERIMENTAL VARIOGRAM 


The experimental drill hole variogram is isotropic, 
while the pit variograms show a slight zonal anisotropy with 
the axes perpendicular and parallel to the benches 
(Figure 5.5, p.196). When h is parallel to the benches, the 
Sill may be lower and the range shorter than for the vario- 
gram perpendicular to the benches. However, the 
perpendicular variogram is so short and close to the L/2 
value that the range may not be real. If it is real, then 
it 1S probably due to the surveyors’ tendency to remember 
the last measurement along the bench and consciously or 
Subconsciously adjust the current one, particularly when the 
contact 1S questionable. On the other hand, the closest 
measurement across the bench would have been taken several 
months earlier and would not have any mental effect on the 
present reading. This effect will increase the range and 
decrease the sill measured parallel to the benches, but it 
is not large, approximately half of the sill and maybe up to 
half of the range. 


The variogram was considered to be isotropic for several 
reasons. Firstly, the anisotropy is not large and kriging 
is robust with respect to the variogram. Secondly, the 
range, with the largest anisotropy, does not affect the 
model as much as the sill with the lesser anisotropy. 
Finally, the pits are more or less perpendicular to each 
other (Figure 1.11, p.24), which means that the aggregate 
variogram of all pits together will tend to be isotropic. 


The drill hole variogram is also isotropic for several 
reasons. Firstly, the drill holes were not logged sequen- 
tially in one direction, so there is no 'memory effect’. 
Secondly, any anisotropy due to location, for example 
perpendicular and parallel to the paleochannel or to the 
wedge-edges of the partings, is too small to be measured. 
This suggests that the underlying variogram in the pits may 
also be isotropic. 
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APPENDIX 6 
VARIOGRAM MODELLING METHOD 


Flint (1978) wrote a program to model variograms, which 
was modified to remove some of the inherent guessing in 
modelling. A range of points is input to the program, which 
then estimates the structures. Finally, the parameters are 
adjusted by trial and error to improve the fit. 


For a Simple variogram with one structure, the process 

is as follows. 

1. The value of the population variance is used for the 
Siig 

2ePeThe first few’ points are input” into the’ program,. which 
fits a least Squares line through them. This line 
intersects the sill at 2a/3 (Figure 1.2, p.140) and the 
vertical axis at the nugget effect (Co). 

3. The parameters are adjusted until a reasonable fit is 
achieved. 


When there are several nested structures, as in this 
study, the best method has these steps. 

1. The process above is begun at the largest structure with 
the sill equal to the population variance. This time 
though, the best fit line intersects the vertical axis 
at the sill of the next smaller structure (Figure 1.10, 
pets) ¢ . 

2. The process above is repeated using the sill determined 
in the previous step. The last step determines Co. 

3. The initial guess at the parameters 1S adjusted to 
improve the fit. 


The variogram parameters Cy, C and a, are not equally 
easy to eStimate. The range is difficult because it 1S at a 
tangent pornt + *furthermore® it is particularly difficult 
when the structure is defined by only one or two points, as 
are the short range stuctures in this study. A good esti- 
mate requires at least three or four points, but fortu- 
nately, the range does not greatly affect the model. On the 
other hand, the sill of the long range model is very easy to 
estimate, as are the other sills, though they all usually 
need some adjustment from the first estimates. The nugget 
effect has a large effect on the model and is the hardest to 
estimate in this study, as it is the result of an 
extrapolation of a few points to intersect the vertical axis 
at a small angle. It must uSually be guessed by visual 
estimate rather than trusting the results of an algorithm. 
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APPENDIX 7 
OUTLINE OF THE THEORY OF KRIGING 


Kriging iS a method for estimating the average value of 
a regionalized variable from data distributed over a region. 
Each data value can be either at a point or an average over 
a region. Similarly, the kriged average can be at a point 
Or over a domain. 


7.1 VARIETIES OF KRIGING 


Many varieties of kriging are available, the choice 
depending on the assumptions made and the properties of the 
data. The three mentioned in this study are listed below. 
1. Ordinary kriging: the expected value of the data is 

unknown but stationary. It is the same unknown value 
over the whole domain, which in mathematical terms is: 


E[Z(x,)] =m iolerac} 


This variety of kriging is used in the present study. 
2. Universal kriging: the expected value of the data is 
both unknown and not constant but varies depending on 
the location, in other words there is a adrift. 
3r, esbognormal kriging: sthis. is. sordinary..kriging .on 
lognormally distributed data. 


7.2 OUTLINE OF THE THEORY OF ORDINARY KRIGING 


The operation of kriging has two parts: eStimation or 
the calculation of some weighted average value of a variable 
in an area, called the domain, followed by the estimation of 
the variance of this estimate. 


The average value calculated by kriging is the BLUE. 
The calculation is done by finding those weights for averag- 
ing the samples that will minimize the variance of the esti- 
mate. This is accomplished by solving a system of linear 
eguations called the kriging system. A detailed description 
of the system and its development is in Journel and 
Hui 1 bGLegts 1019/8 wD u5 0.3) panda David 21S hie. e3.1 be 


The mining industry uses a variety of other estimators: 
1. the polygon method, used at the Boundary Dam Mine, 
2. simple average, called the extension estimator in geo- 
Statistics, 
3. distance weighted averages. 
These all use more or less arbitrary weights; therefore, 
they cannot guarantee a minimum variance estimator. In 
addition, they take into account neither the statistical 
Structure of the data, in particular the variogram, nor the 
Spatial arrangement of the data. 
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7.2.1 General Problem: to Find the BLUE 


Samples Z(x,), i = 1,. 2, 3, °** "nm have been” measured 
and Z is to be estimated: Z is the true average value over 
some domain. The estimate (Z*) should be the best possible, 
which in this case is defined as the BLUE.f The estimator 
is determined by the krige system of equations which will be 
developed below by considering the BLUE in this order: esti- 
mate, linear, unbiased and best. 


The objective is to estimate the true average value (Z) 
by Z*, which will be some combination of the sample values 
Zixee) AS in equation (lect D. 1s 5. 


Z* could be any combination of Z(x;); for example, 
Seat One oe De Ss onOrmeduation leo}, ‘Delos. “But. in this 
Study a linear combination is used: 


Zee na (2x4) iret 
Where: 
a; = the weights for a linear combination. 


A linear combination is used for the reasons listed on 

page 135 above and also because nonlinear estimates require 
extra assumptions beyond the intrinsic hypothesis which are 
ditticuley tot verily... Thesoroblem now reduces to finding 
those weights (a,) that will give the best estimate. If 

a, = 1/n, then Z* is the simple arithmetic mean of the 
Sample data or the extension estimate. 


The estimator should be unbiased, which means that in 
the long run Z* will systematically neither over- nor 
underestimate Z. The unbiased condition is satisfied 
(Davide 977," De coo) 1 tne aucrinsic, nypothesi's 15 Satis— 
fied and if: 


T= yeah Ee ht: 


Finally, the estimate should be the best one, which is 
called the krige estimator (Z{%). The best estimator in this 
case is defined as the estimator with the minimum variance 
of estimation. This means that over the long run the esti- 
mate (ZX) will come closer to the real value (Z) than any 
other eStimate’ (Z™). 


The variance for any estimation (02) is calculated by 
(David, 19775. Dp. coo" 


eee agy( he a) ananty (out). yu VV) {7.4} 
i eS 


+ The Best Linear Unbiased Estimator is defined on page 135. 
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Where: 

V = the domain over which Z is averaged, in general 
LieiSsa volume ;nbut».inw this study sites an.area. 

UV + aviow the average variogram between regions V, and V2. 

This term is estimated by: 


IR 


yb Viren Vi ) CUA hig) Eee Kapex gh) tego: 
Ue 

Where: 

x jj iSs an Vi«s 

aS ain. Ve. 

The summations are over all pairs of points, one from each 

region. The average variogram 7(V,,V2) is calculated by the 

following steps. 

1t. Consider all pairs of ,points},.the firstex, .in,V,:and. the 
Otherex, in V2. 

2. Calculate the vector separating them, x;-x;. 

3. Calculate the variogram at this separation, y(x;-x;). 

4 Average all these variogram values between each pair of 
Points, Cl/nin, )zhy (x oxy) 

is 


The two regions V, and V2 can have these forms: 

1. V;, = V2, in other words, the average variogram between 
all pairs of points in one area, the domain. This is 
called the average variogram of the domain and abbrevi- 
ated to 7(V,V). 

2. V;, = X;, in other words, one of the regions has been 
collapsed to a point x;, usually a sample location. 
This 1S called the average variogram between the sample 
ix.) Vana eene omMaineana 1S written, ase y(x7,)). <The 
average variogram is properly defined in terms of all 
points inside the sample being one of the pair of 
points. However, if the sample is small enough with 
respect to the domain, then one point in the centre of 
the sample is uSually accurate enough. The smoothing 
introduced by this approximation is small. It is the 
Same as the smoothing due to non-point support 
(Figure 124). 1p.144)7 

3. If both regions have been collapsed to points or samples 
at x; and x,;, then this is called the average variogram 
between the samples Z(x,) and Z(x;) and is written as: 
Ae ae cee 


7.2.2 Auxiliary Functions 


Calculating an average variogram that includes a non- 
point region is very difficult.+ It has been explicitly 
solved for only a few cases, which must have all the follow- 
ing characteristics: 

1. particular model, either the spherical, exponential, 


+ In the general case it involves a sextuple integral of the 
average variogram, one for each dimension in each region. 
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linear or log (de Wijsian), 

2.) “rectangular “domain, 

3 Simple spatial relationship between the sample and the 
domain. For instance the sample must be at the corner 
or along the side of the rectangular domain. 

The resulting average variogram called auxiliary functions. 

Unfortunately, the samples and domains in this study do not 

fall into the simple patterns assumed by these functions, in 

which case it 1S necessary to evaluate the ¥ terms using 
numerical integration. 


7.2.3 Kriging System of Equations 


For the estimate to be the best, o2% must be a minimum 
with respect to the weights aj;, or: 


OO Ge, Nee OS Orem = 2,8 348 see on ia. Oo 
Where: 
d02/da, = the partial derivative of o2 with respect to aj. 


These equations, along with the unbiased condition of 
equation {7.3}, p.200, are n+1 equations in the n unknown 
a,'S. This over-determined system of equations 1S minimized 
uSing the standard Lagrange technique, which introduces an 
extra variable called the Lagrange parameter (u) and 
minimizes a new system of n+1 equations in n+1 variables: 
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These reduce to: 


Za, 7(x;,x;) OPE Axe, AVS) , for in = li, Ve oF a 8 8 
J 
a ; = | {77,163} 


These. can also betwritten’ in Matrise Lorn: 
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Or in simpler notation: 
KeA = D 7294 


Where: 

K = the krige matrix: ~See equation {7.5}, p.201 for a 'defi- 
nition “ofthe 7 (x7 7x ,)*terms: 

A= [a,, Sree een ult 
the solution, a column vector of the Lagrange parameter 
and the required weights for the krige estimate. 

D= Lyx Vv V7 (xs, V2), oo WAX 5 ee 
a column vector of the average variograms between each 
Sample x; and the domain V. See equation {7.5}, p.201 
for a definition of the 7 (x,,V)oterms: 

This is called the krige system of equations. The required 

weights (A) are determined by solving this equation: 


A = K''eD 710i 


Where: 

K-) €-the°Pinverse of K: 

The components of A are uw and the weights (a;) needed to 
calculate the reserve estimate that will minimize its 
variance. 


7.2.4 The Estimators and Their Variances 


The krige estimator (2%) is calculated by: 


Zk = La, 2(x;) Cre 1t 
Where: 
a; = the elements of A, the solution vector, without u. 


The variance of the krige estimate (07) is calculated by: 
Of Ma LaF VCR EV TT Coe ya) P12 3 


Equation {7.4}, p.200 could also be used but that equation 
has more summations and is slower to compute. 


The extension estimator (Z%) or the arithmetic average 
of the data is calculated by: 


ZX = (1/n)2Z(x; ) PUES HESS 


The extension variance (0%) is calculated by substituting 
aTPe="17n-into-equation@t724ty pr200°totgive: 
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APPENDIX 8 
DESCRIPTION OF THE KRIGING PROGRAM 


This appendix contains an outline of the operation of 
the kriging program with the emphasis on the modifications 
madescto, original version from Flintw(1978), brief lists of 
the input and output of the program, the best values to use 
for most of the input parameters, and the size of 
inaccuracies in the results due to the calculation method. 


8.1 OPERATION OF THE KRIGING PROGRAM 


The basic kriging program as written by Flint (1978) 

has the following steps: 

1. a input of the parameters and data as listed in Appendix 
8.3.1 and Appendix 11 respectively. 

2. calculation of the krige:matrix/A,;" and vector D, using 
these steps: 
A. construction of the numerical integration grid, 
B. computation of the 7 terms by numerical integration, 
Caigeonstructiony ofy hands De from, thes 7A x j4'Xy and 

¥(x,,V) terms respectively, 

3. inversion of the matrix and calculation of the weights 
vector 4, using IMSL subroutine LEQTIF, 

4, calculation of the krige estimate Z{ and its variance 
Ox, 

5. output of the results as listed in Appendix 8.3.2 with 
an example in Appendix 13. 

In addition, the program calculates the area of the domain 

and a default cell size. 


8.2 MODIFICATIONS TO FLINT'S KRIGING PROGRAM 


Flint's kriging program was modified to perform the 
numerical integrations twice instead of once, and also to 
calculate the extension estimate and its variance. 


Auxiliary functions as described by Clark (1976) cannot 
be used because the domains to be estimated are irregular, 
not rectangular. Instead, to calculate 7(x;,V) and ¥(V,V), 
numerical integration must be done which introduces 
approximation errors. 


This numerical integration is not the normal variety 
where a function is evaluated at each point in a region; 
instead it is the evaluation of a function (the variogram) 
for the distance between each pair of points in a domain. 
Integration is done by setting up a grid over the domain and 
evaluating the variogram either between the sample and each 
node for y(x,,V) or between each pair of nodes for ¥(V,V). 
The smaller the grid cell size the closer the calculated 
value is to the true value. Normally the grid is completely 
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inside the domain, and consequently, the calculated value 
will always underestimate the true value. This is because 
the grid will tend to use smaller h's than if all of the 
points inside the domain were used. Shorter h's mean that 
the variogram values will tend to be less and the average 
variogram lower. 


The program was modified to calculate the value twice: 
once as above with all nodes just inside the domain and then 
again with the same grid extended by one node to just 
outside the domain. In the second calculation the average 
variogram should be overestimated, assuming no approximation 
errors. Thus, the two calculated values ought to bracket 
the true value if the cell size is small enough (about 20 m 
in this study). Finally, the program averages the two 
values for the best guess and reports the percent difference 
between the two estimates for judging the accuracy of the 
result. 


Ordinarily, calculating the average variograms twice 
would not be necessary but in this study the seam thickness 
is so regular that the estimation variances are very small. 
The krige variance is largelyft the result of the subtraction 
of two nearly equal large numbers: ¥(x;,V) and ¥(V,V) in 
equation -7ei2y, "Oe20or. Sltigne, approximatron errors in 
these two numbers, particularly the latter, can produce 
large errors in the krige variance, to the point of a nega- 
tive variance. 


8.3 INPUT/OUTPUT OF THE KRIGING PROGRAM 


Listed below is the input to, and the output from the 
modified kriging program. 


8.3.1 Input Parameters and Data for the Kriging Program 


Input may be entered interactively from terminal or 
from a file in batch. The input consists of: 

1. Sample meaSurements in a file which consist of: 

A. net thickness of the pit measurements or drill 
holes, 
B. location of the samples as X and Y coordinates, 

2. the domain over which the average is estimated. This is 
in the form of a file of the X and Y coordinates of the 
digitized vertices of a polygon outlining the domain. 
The program also needs the FORTRAN format of the above 
two files. 

3. a description of the variogram model, consisting of 
these parameters: 

A. the number of nested models, up to four, though only 


+ The Lagrange parameter (u) is usually small in this study. 
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one or two were required in this study. 
Bees Cone the DUGGeL Grrlecr, 
C. CC; -*°* Cy, the sills of the nested models, 
D. R, *** Ry, the ranges of the models, 
. the cell size to be used in the numerical integration, 
whether the elements of the inverted krige matrix are to 
be tested for accuracy and, if so, for how many digits. 


O1 > 
° 


8.3.2 Output From the Kriging Program 


The program has two types of output: the kriged results 
and a list of the sample covariances, y(x;,,V), with sample 
weights, a,;. These are listed below, with an example in 
Appendix 13. 


The kriging results are presented (both on the screen 

and written to a file) in two portions: 

1. an echo of the input parameters, including those listed 
above along with the following results of some 
intermediate calculations: 

A. the number of samples used, 

B. the area of the rectangle that encloses the domain. 
This is the size of the numerical integration grid; 
it helps in choosing an appropriate cell size that 
will balance computing speed and accuracy. 

C. a default cell size based on the size of the 
enclosing rectangle, 

D. the actual numbers of rows, columns and cells of the 
grid inside the domain, 

E. whether the accuracy test passed or failed for the 
specified number of digits. 

Zum the, results, of, thevkrigang calculations: 

A. a calculation of the area of the domain over which 
the average was estimated, 

B. the kriged estimate, 

C. the extension estimate, 

Both of the estimates above include the following: 

a. the estimate itself: Zi or Z&, 

b... the. estimation variance;. of or o£Z, 

c.4 Lne,covariance. of, the domain: 7(V,V), 

dad. the average covariance of the samples to the 
domain; Za;,7(x,,V), } 


e. the Lagrange parameter; u. 
The last three are used to calculate the estimation variance 
and their values are listed in order to judge the accuracy 
of the variance. If the values are large and close 
together, then the variance will be less accurate. Each of 
the above five values also has a percent difference 
reported. This is used to estimate the approximation error 
in the numerical integration. 
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The percent difference is also reported as a variance 
and standard deviation, calculated as if the over- and 
underestimation were a sample of two. This error in the 
estimate will add to the estimation variance for the total 
error in the estimate. However, the percent difference 
expressed aS a variance waS uSually several orders of 
magnitude less than the esStimation variance and for this 
reason, was neglected. 


The sample covariances and weights are written to a file 
in two portions: 
1. a brief echo of the important input parameters in order 
to identify the results, 
2. the following information on each sample: 
“<n cud. Numbe te, 
ues Gl OC OTHION *ntXi ia Sp eenands Yi rCOOrdaunates, 
sample value; Z(x;), 
. the covariance of the sample to the domain; 7(x,,V), 
the krige estimate weight; aj. 
This output is used to verify the input and to determine the 
contribution of each sample to the estimate and estimation 
variance. Some samples may be screened or otherwise 
contributing little, and so could be dropped from future 
estimations. 


OO > 


8.4 SENSITIVITY OF THE KRIGING RESULTS TO THE INPUT 
PARAMETERS 


Kriging and the program to carry it out are complex, as 
result, it 1S not possible to accurately predict the effects 
of input parameter variability on the output. Therefore, 
the program was tested using either artificial data with 
known results or data from this study for a more realistic 
teste. 


The following input parameters were empirically tested: 
1.  mumber of samples to be kriged, 
2. accuracy of the variogram model parameters, 
Bee pOTcanGe lisse. 
The other inputs, the data and the domain, were thoroughly 
checked and assumed to be correct. 


8.4.1 Number of Samples to be Kriged 


The number of samples to include in the kriging is a 
compromise between accuracy and computing cost. More 
samples will produce a more accurate estimate, but the 
number of samples is one less than the number of 
Simultaneous equations to be solved, so the costs rise 
rapidly with the number of samples included. All samples 
inside the domain should be used, plus those samples very 
near the domain. But data beyond the largest range is 
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uncorrelated with the domain, and having little effect on 
the reserve estimate, can be ignored. Between these two 
limits, as samples get further from the domain, they are 
screened out, that is, their effect is reduced by closer 
Samples (David, 1977, p.258). In this study the screen 
effect 1s strong because the relative nugget effect is low. 
The minimum number of samples that produces a stable esti- 
mate and variance was determined experimentally. 


Computer hardware and software impose some upper limits 
on the number of samples. The program can handle a maximum 
of 200t samples on the regular CPU and 180¢ on the much 
cheaper array processor. Unfortunately, most domains 
studied had more than 180 samples. As the number of samples 
increases, the more ill-conditionedtt the krige matrix 
becomes, until eventually it is effectively singular and the 
Subroutine cannot invert it. This was only a problem when 
kriging pit 1 because it is so large. 


Each of the pits was kriged using all the drill holes 
inside the pit plus all the drill holes inside successively 
larger aureoles around the pit (Figure 8.1). The values of 
the krige and extension estimates and their standard 
Geviations were plotted against both the distance of the 
aureole from the pit boundary and the number of samples 
kriged. 


Pit 2A (Figure 8.2) is representative of all the pits. 
The krige estimate and variance settled down after a certain 
aureole size and number of samples was reached. Between the 
different pits, this point was more constant in terms of the 
aureole size than in terms of the number of samples 
(Table 8.1). The kriged estimate and variance settled down 
at about 40 samples. The extension estimate and variance do 
not settle down, possibly because the numbers become more 
biased with more samples, though 40 samples is also a 
reasonable number to calculate them as well. 


Samples out to 200 m away from the domains were used as 
a conservative limit. When estimating with drill holes, 
this aureole around any pit or bench included the minimum of 
40 samples required for stability. However, this criterion 
produced more than the software limit of 200 samples for 
estimates using pit data, so only the nearest 200 were used. 
If all available pit samples were used, the results should 
not change very much for the following reasons. 
1. The seam is very regular and there is no drift. 
2. The nugget effect is low, so distant samples will be 


+ This limit is easily changed. 

+ This limit is impossible to change. 

tft 'Ill-conditioned' means the matrix is nearly singular, 
making A very sensitive to small changes in the elements of 
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Figure 8.1 Aureoles of data used to determine the best 


number of samples for kriging 
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(m 


EXTENSION 


(m) Estimated Thickness, 


Estimation Standard Deviation, 


OF 
Sample size 
t2 18 25 34 42 53 63 78 93 
Inside 100 200 300 400 
pit 2A Aureole size, (m) 


Figure 8.2 Thickness estimate and standard deviation 


versus the aureole size for pit 2A 


Table 8.1 


Pit 


Pit 


Pit 


Pat 


2A 


2 


Sample size and aureole distance required to 


produce a stable estimate 


Kriging 


Sample 
size 


42 


38 


42 


Aureole 
distance 


(m) 


50 


200 


200 


Extension 
Sample Aureole 
Sizef distancef 
(m) 
50 70 
26 150 
38 1:70 


+ The sample size and aureole distance are averages of where 
the estimate and itsS variance stabilizes. 
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screened. 
3. The estimate is stable for this number of samples 
(Pigurée 8.2; 3p. 10). 


8.4.2 Accuracy of the Variogram Model Parameters 


Kriging iS very robust with respect to the model chosen 
for the variogram (Journel and Huijbregts, 1978, p.233 and 
David, 1977, p.113) and is moderately robust with respect to 
the parameters chosen for the model. It was not possible to 
test the effect of different models on the results because 
both the modelling and the kriging program used can handle 
only a spherical model. 


The effect of the parameters was tested experimentally. 
Different sets of parameters were input to the kriging pro- 
Gram with all other inputs kept constant and the results 
compared for the pit data and again for the drill hole data 
(Table 8.2). Fairly large changes in the parameters could 
be made without large changes in the results. The estimates 
are within 0.jio0 to 0.50 of each other, or expressed in 
thicknesses they are within 1.5 cm over a total of 3.5 m; in 
other words identical. The estimation variances, expressed 
as a standard deviation, were much more variable, in some 
cases up to 50%. However, the values are so small that the 
absolute differences are also small. | 


See i3171€ell ‘Size 


Cell size is the distance between the nodes of the grid 
used for the numerical integration. The choice of the cell 
size 1S a compromise between increased accuracy and rapidly 
increasing computing costs as the cell size is reduced. The 
effect of the cell size was checked experimentally using 
data for which the results were determined by exact 
calculation using auxiliary functions. Brooker (1979), 
Clark (1979a, p.106) and some hand calculations provided the 
worked examples. The program actually uses the number of 
cells rather than the cell size. The number of cells was 
varied and the results compared to the exact values using 
auxiliary functions. Below about 50 cells the results were 
very erratic. Above about 200 cells the results were 
Stable; for this reason, 200 to 300 cells were used in this 
study, which produces a cell size of about 20 m. 


Flint (1978, p.67) also needed about 200 cells. Clark 
(1979a, p.121) states that by general agreement the number 
of cells should be in the range 64 to 100. David (1977, 
p.280) writes that 100 cells are necessary to calculate j¥ 
terms with a few percent precision but that fewer than 10 
are needed for a stable estimate. This study required more 
cells because the thickness population variance is so low, 
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variogram model parameters 
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study 
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Kriged drill hole data 
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Robustness of kriging with respect to the 
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which makes the calculated results more sensitive to the 
inputs. 


8.5 INACCURACY OF THE RESULTS DUE TO NUMERICAL 
APPROXIMATIONS 


The accuracy of the computed results is affected by 
numerical approximations at two places in the program: 
numerical integration of the ¥ terms and inverting the krige 
matrix. 


8.5.1 Inaccuracies Due to Numerical Integration 


The accuracy of the numerical integration was tested 
using the same data as in Appendix 8.4.3. so that the 
computed value could be compared to the exact value. Krige 
estimates and variances were computed for several cell sizes 
(Figure 8.3 and Figure 8.4)t. The percent absolute 
difference between these results and the result calculated 
by Clark (1979a, p.109) using auxiliary functions was 
plotted as the 'error'. The 'difference' is the percent 
absolute difference between the two approximations computed 
by the machine; the over- and under-approximation as 
described on page 205. To judge the accuracy of the 
approximation for real data the ‘error’ must be known, but. 
the machine can only compute the 'difference'. Fortunately, 
if a large enough number of cells is used the two machine 
answers should bracket the exact value. AS a result the 
‘difference' should be greater than the ‘error', so that the 
‘difference’ can be used as a conservative estimate of the 
accuracy of the results. 


The variance behaves as it should but the percent error 
of the krige estimate iS Somewhat erratic and not always 
less than the percent difference. This may be due to the 
particular location of the grid with respect to the domain. 
However, the ‘error’ and the 'difference' for the estimate 
are always within an order of magnitude of each other. The 
absolute values on the log-log plot are so small that the 
program can be considered to be accurate and the percent 
difference reported by the program a good indicator of its 
accuracy. 


This study generally used about 200 cells for which the 
kriged average was always within 1% of the true value and 
usually closer. The krige variance was always within 10% of 
the true value and often within 2% of it. 


The extension estimate was always exact of course. The 
extension variance is estimated in a Similiar, though not 


t Note that both of these figures are on a log-log scale. 
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Figure 8.3 
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identical, manner to the krige variance (equation {7.4}, 
p.200 and equation {7.12}, p.203 respectively). Therefore, 
it will have about the same size error as the krige 
variance. 


8.5.2 Inaccuracies Due to Matrix Inversion 


The accuracy of the krige matrix inversion is tested by 
the inversion subroutine and reported as the number of sig- 
nificant digits. As the number of samples increases the 
number of significant digits decreases; it is usually two 
digits (about 1%) and occaSionally one or three digits 
(about, 10% or 0.1%) Thiswis the accuracy of the individual 
weights (a;) calculated, of which there are 50 to 100 for a 
kriged average in this study. 


The estimate iS a linear combination of these a;'s and 
David (1977, p.281) states that "+*+-linear estimates are 
deceptively robust." This is because the calculation errors 
in the weights should tend to cancel out, with the result 
that the calculation error of the kriged estimate and 
variance will be less than the significant digit error and 
much less than the inaccuracies due to numerical 
integration. For this reason and because the results are 
close to the exact values for the test data, any 
inaccuracies due to inversion of the krige matrix by the 
Subroutine were ignored. 
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APPENDIX 9 
FORMULAS FOR THE PROPAGATION OF ERRORS 


Estimates and variances of the reserves (T) and the 
recovery (R) are calculated from parameter estimates of the 
input variables; A, d & Z and CLO & T respectively, using 
the formulas below. The general expression is considered 
with any dependent variable (T or R) related to the input 
Var VaDLeS. (X1,— Xo moe, Kaus Dyutne relation: 


mmeils Xr) ck cya Xe ce eX el LO} 
Ors 

MeN tie X21 kor eX) eyes 
Where: 
T, ®, x; = general variables. 


E[-] and VAR[-] are estimated for two cases, depending on 
whether the input variables (x;,) are correlated (for ®) or 
not (for T). For the tonnage they are uncorrelated, but for 
the recovery they are highly correlated. These formulas, 
given below and developed by Hahn and Shapiro (1967, p.252), 
all assume that fourth and higher order terms effectively 
disappear. The high order terms in the formulas below are 
either zero or negligibly small, so higher terms can be 
ignored. 


9.1 ESTIMATE AND VARIANCE WHEN THE INPUT VARIABLES ARE 
UNCORRELATED 


The thickness (Z*) has no drift; consequently, the 
expected value is everywhere the same, in which case it will 
not vary with the area of the domain. The density was 
treated aS a non-regionalized variable, which implies that 
it too has no drift and therefore will not vary with the 
area of the domain. Finally, Z* and d must be assumed to be 
uncorrelated since there 1s not enough density data to 
determine .COVI Za)... Lm tact,athere will be a oslioght corre- 
lation because portions of the seam with more clay bands 
will probably be both thicker and denser, since clay 
compacts less and is heavier than coal. However, the coal 
at the Boundary Dam Mine has a low ash content (Table 1.1, 
p.27) and therefore, the covariance term should be small. 
The input variables; A, d and Z”, are uncorrelated and so 
any covariance terms vanish. 


When the input variables are uncorrelated, that is: 
COM LXeaXar ese eta tOlemel Elie) 


E[T] is estimated by: 
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Pi EX ay Soe LX nl) 
+ 22(0°T/0x7)VARLX; J These U3 
In this case, where T = T; 
Pes Ix 
therefore: 
a-T/ox* = 0 
Which means that the estimate of E[T] reduces to: 
BOreh= SAE Yr) peex 2 tee eB ard {9.4} 
The variance is estimated by: 
VAR(T}] = Z(0T/dx;) *VARLX; ] 
+ L(0T/0x;)(8*T/ox7 ELXi-ELxiJ]* + ++> ise 
Where: 
07T/dx? = second partial derivative of T with respect to x, 


As for equation {9.3}, p.219 second and higher order terms 
vanish to give: 


VARDT Pe= Sz CeNPE-T A) WARIO] AP POPOL SET MeZT hes 1976) 
Where: 

= T 
Ny =A 
X2 =a 
a) 42 +7 el bher SZ PuOriezZy 


This can be expanded to: 
VAR(T) = d2725-VARLA) 0A 2" VARIG ot AZd-VARIZ 19244 


Where: 

WVAR(Z" } = of “or 0% depending on whether 2" = Ze or Zt. 

An alternative version of the above is expressed in terms of 
the relative variances: 


VAR(T AT #='"VARL A] JA* 86 VAR[ d h/d= “P°vAR[ 2 O)Y 7" {9.8} 


Where: 
VAR[x]/x? = the relative variance of x. 


Relative variance, always reported as a percent in this 
study, is a convenient way to express the error in a varia- 
ble. This error can be split into simple additive 
contributions from several sources, unlike the raw variance. 
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9.2 ESTIMATE AND VARIANCE WHEN THE INPUT VARIABLES ARE 
CORRELATED 


The mine engineers expend considerable effort to raise 
the correlation between the in-place coal and the coal- 
loaded-out as close to one as possible. Consequently, the 
recovery and itsS variance will be biased if calculated using 
formulas from the previous section. Unbiased values require 
the formulas from this section. 


9.2.1 Mean of the Recovery 


When some of the input variables are correlated, that 
is: 


GOV bx i), x; Part 7] fon isome. ij #63 
then E[R] is estimated by: 


E[R] = FIBLX‘a] -Elx ak, a ABAX ads 


i $22 (37R/8x 8x ;.) COVI Xi nx; ] tera {9.9} 
ba 
When ox? c= GLO, ithen 'd*RAGCLO* «= 40 pvanddahl, termstinvoiving 
it vanish. When x; = T, the terms are negligibly small. 


This can be demonstrated with an order-of-magnitude 
calculation using the following data derived from 

Bates i hO se pabia: 

1 eCLO= Tre 25410 tonnes 

22° VRE) = bc10P/etonness: 

SB meCOVIT, CLO] = txit0mestonres +. 

COV[T,CLO] is estimated by remembering that the correlation 
coefficient between reserves (T) and the coal-loaded-out 
(CLO) will be approximately one, therefore: 


1 = COV[T, CLO] //( VARI EIVARICEO?) 
Ores 
COVET CLOIS y (VAR TIVARLCLO]) 


The coal-loaded-out figures are reported to six significant 
figures (Table 5.5, p.99), which implies that: 


VAR LCLO. =e Oat 


There are two terms involving T in equation {9.9}, p.220. 
When i = j then: 


CO2R/ oT) VARUT = 02007 BaVARET |) = iixi0s) 


and when i # j then: 
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CO2R7 OTOCLO) COViT sy chOdm=e— (400772) COV LT. CLO] =F 1 eae 


I 


Both terms are negligibly small when compared to a typical 
recovery of about 100% with an estimation variance of about 
4%. As a result, equation {9.9}, p.220 reduces to: 


E[R] = FLELXS Bb ex2e ke, BCXR) J £9410} 


9.2.2 Variance of the Recovery 


The variance is estimated by: 


VAR[R] = LE(8R/3x;)(8T/3x;)COVEx:,x;] + Bt ee* {9.11} 
ij 
Where: 
COVERY XS =-VAR[X,J]°9ife2°= J 
9) = higher order terms in (dR/dx;)(8?R/0dx jax, ) 


This formula makes the same assumptions as were made above, 
and are met for the same reasons. The first term can be 
decomposed into a sum of variances (when i = j) anda sum of 
true covariances (when i # 3). Using the same data as was 
used for the mean, the true covariance terms are approxi- 
mately 1x10°’. The partial derivatives in 2 evaluate to 
approximately 1x10°'*, which the attendant expectations will 
not substantially increase. All these and higher order 
terms are negligible compared to a typical recovery 
estimation variance of 4%. Equation {9.11}, p.221 with the 
appropriate variable substitutions becomes: 


VAR[R] = (CLO/T?)*VAR[T] “9.124 
or in terms of the relative variance; 


VAR[R]/R? = VAR[T]/T? 19433 
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APPENDIX 10 
VARIANCE OF THE AREA OF PROJECTED PRODUCTION 


The reserves estimates require the variance of all the 
input variables, in particular the area. If the area is not 
predefined, then it has an error associated with it and 
therefore, VAR[A] > 0 and equation {9.7}, p.219 or 
equation {9.8}, p.219 must be used. In this case it is nec- 
essary to eStimate the area (and its variance) of the seam 
removed by the mining operationt in a given amount of time. 
One year is used as an example. 


The mine plan (described on page 38 above) has some 
benches with dates on them showing the limits of excavation. 
These are rarely at one year intervals or any other constant 
unit of time. In order to be able to estimate the area 
uncovered in one year, and to be able to use all the availa- 
ble data, the amount excavated was regressed against the 
time to excavate for all the dated benches. 


This method assumes that the dates on the map accurately 
show the extent of the excavation. There is no reason to 
Suppose that the mine surveyors have been inaccurate. The 
date could mark the time of a Survey rather than the time of 
the advance, but this seems unlikely. This method also 
assumes that the data shows the limits of coal removal and 
not just the limits of overburden removal. If the dates are 
for the latter, then the coal removal will usually lag by a 
fairly constant distance behind, and therefore, the length 
of coal removed along the bench will be about the same as 
the length of the overburden removed in the same time 
period. This possibility will increase the variability of 
the data. 


The amount excavated can be measured using three 
methods, in order of increasing complexity these are: 

1. Measure the length of the bench, then assume a constant 
width to convert to area. 

2. Measure the area excavated directly. 

3. Measure the volume of material excavated (mostly 
overburden), then assume a constant overburden depth to 
convert to area. 

The method used should be simple, but more importantly it 

should produce as low a variance for the area as possible. 


10.1 DIRECTLY MEASURING THE AREA OPENED UP 
The area excavated on the dated benches could be 
measured directly. However, the extra effort of planimetry 


was not justified for the following reasons. Firstly, 


+t A description of the mining methods can be found on 
page 34 above. 
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planimetering is not very accurate, particularly on areas 
with a large perimeter/area ratio, such as the benches. 
Secondly, as shown below, regressing the length of the bench 
was accurate enough. 


10.2 VOLUME AND LENGTH OF THE AREA OPENED UP 


The most reasonable method is to regress the volume 
excavated against the time required, since the equipment 
moves a volume, not an area. The volume of the dated 
benches was determined from the length along the bench, the 
average width determined below and the overburden thickness, 
all of which were recorded on the mine plan. 


The thickness of the seam has a low variance and is 
small (3 m) compared to the overburden thickness (40 to 
TOim)? consequently ,~it will not affect the variability of 
the total volume removed. 


The overburden thickness is a regionalized variable; 
therefore, its average thickness under each dated bench 
could be determined by kriging. However, the thickness 
varies only slightly under each bench, so the extra effort 
did not seem worthwhile. The overburden measured by all 
holes in and immediately around the bench was arithmetically 
averaged for the overburden thickness in the bench. The 
volume of material removed daily was calculated, then the 
mean, variance and coefficient of variation of the volume of 
the benches from each pit were estimated. Similarly, the 
coefficient of variation of the bench length was calculated 
and compared in Table 10.1. The variability of the bench 
volume is as large as the variability of the bench length 
for each pit; thus, no increase in accuracy 1S gained by 
measuring the volume of the benches rather than the length 
of the benches. 


10.3 REGRESSION OF THE LENGTH OF THE AREA OPENED UP 


The area excavated each year was estimated by 
regressing the curvilinear length along each bench (L) 
against the time (t) required to dig it (Figure 10.1). The 
average length per day of each dated bench is unimodal; 
therefore, a single machine or several with the same 
capabilities must have been used. The dated benches that 
Spanned more than one cut (4 out of a total of 13) were not 
Significantlyt below the regression line (two were above, 
two below); thus, the time to turn around the equipment and 
reverse the direction of excavation was not important and 
was not corrected for. 
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Table /10 <1 Coefficients of variation of the bench 


length and bench volume 


Pit Number of Coefficient of varia- 
e20n 

dated 

benches Length Volume 
Pit 1 5 59% 60% 
PitezAt = 2% el 
Bat. 3S 5 49% 58% 
All pits 3 49% 58% 


+ The thickness was constant over the all of pit 2A so the 
variability will not be different between the length and the 
volume. It was included in order to compare the variabilty 
between pits to within pits. 
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Figure 10.1 Regression of the bench length against the 


time to dig the bench 
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Two models were considered for the regression line, a 
linear model: 


Epa pte + sa Pale Oaeeal il 


and a power model: 


Eats WiCl2) 
Where: 
L = the curvilinear length along a bench. 
t = the time required to remove a bench. 
a, B, @ = regression parameters. 


The best fit regression parameters are listed in 
Table 10.2. The confidence intervals on r? have a large 
overlap; therefore, the models fit equally well. The best 
fit equations are: 


Es 
ie 


So bet Se KO 1307.13; 


and: 


tt 
R 


= aze°. 12 {10.4} 


The linear model has a positive value of a (5.5x107), which 
is both significantlyt non-zero and physically impossible. - 
Even though the linear model has the advantage of 
Simplicity, the power model was chosen for the following 
reasons. 

1. The best fit straight line 1S incorrect because of its 
impossible intercept (a). 

2. The data may be lognormally distributed since they are 
squeezed towards zero, which is an absolute bound on the 
data. Lognormally distributed data is linearly 
regressed by means of the power model. The power model 
is linearized by a log transformation to: 


intel =ee ini le sea HeOr ot 


3. The power curve has a more reasonable shape. Firstly, 
it passes through the origin (practically speaking). 
The true curve probably intersects the time axis at the 
time required to set up or turn around the operation at 
the end of a bench. Secondly, the curve is concave 
down, which implies that the longer the equipment is 
operating on a bench, the more likely a major breakdown 
becomes and the lower the average and instantaneous rate 
of advance becomes. 


The estimated advance in one year is 3.28x10* m, using 
equation {10.4}, p.226. The problem with using a power 


ta = 0.05 
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Table 10.2 


Parameter 


Odd Noha Oe 


mr OG 


Z2y), 


Regression parameters of the bench length 


with the time to dig the bench 


Linear model 


Sox 2 om 


6.2 m/day 


Power model 


47 m/day 
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model is that the variance on an estimated advance is 
asymmetric. The standard deviation is +1.04x10°? m and 
-0.79x10% m. Unfortunately, equation {9.7}, p.219 requires 
a symmetric variance. Therefore, the absolute values of the 
plus and minus standard deviations were averaged and squared 
for an eStimate (8.37x10° m*) of the symmetric variance on 
the length. This is probably a conservative estimate Since 
the variance uSing the linear model is an order of magnitude 
less: 7.44x10* m*. The average bench length and its error 
are Bisted in Table 1073. 


A few hundred bench widths were measured at evenly dis- 
tributed points on the mine plan. The statistics for the 
bench width were estimated using classical statistics. The 
average bench width and its error are listed in Table 10.3, 
Dezou. 


The average area excavated in one year was estimated 
using the preceding values for the length and width along 
with equation {9.73}, p.219, with the appropriate variable 
changes. The average and its error are listed in 
Table 10.3, p.229. The formula assumes that the length and 
the width are independent. In fact, they are not strictly 
independent, since in one year more of a narrower bench 
could be cut. This covariance should be small since the 
bench width is so regular. The width contributes 1.6% to 
the relative variance of the area with the rest coming from 
the length. 
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Table 10.3 Statistics of the bench length, width and 
area 

Mean Variance Coefficient Relative 

of variancet 


Variation 


Bench length excavated in one year 


S26 x10 am Geeahy cok Oe gms 28% 7.8% 


Average bench width 


31.8 m 1632 %m4 1 WaGs 


Bench area excavated in one year 


1IO4SH16 3 IMF 130200" oma 31% 9.4% 


+ The relative variances are smaller than the coefficients 
of variation because both are actually decimal fractions 
less than 1, but written as a percent for convenience. 
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APPENDIX 11 
SAMPLE OF DATA 


11.1 SAMPLE OF DATA FROM PIT 2A 


ID = Sample number; the first digit is the bench number, 

Coordinates = Metres north and east from the SW corner of 
NE-33-1-8-W2M. 

F = A flag indicating the data source: 1l=pit, 3=drill hole. 

N = The number of splits. 

All thicknesses are in metres; ***** indicates missing data. 


ID Coordinates FN Thickness Coal Parting Coal 
X "6 gross net bottom top 
201 1750.4 S44 s744 4 KAKKKK  KKKKK * KK KX 
2029677 38918 S233Seu Z 540390 q3e932 Si Oiro wy) alter Om m0 oalre 
29321272598 S0sSeai 2 44542 93-962 3.048 0.610 0.914 
20457171338 9284454 2 4.420 3.993 3h028 909427 0.974 
20594120008 O5641m4 2 42298 #324932 Se0us nO0ews66 30.914 
20692168 74:2 984¢43i1a%7 2 4.206 ****x EXKEK 0,244 *kKEX 
Z20/ 493673855 1042360% 2 44359. 422415 Se2os 20. 246" 00853 
ZOS%E16608%4 TO0S9E0 74. 4 334962 932962 33962 
209 4e164 7524 ' 40662154 2 AeStt 148206 So Dome soo FU. 65s 
240793634220 209.138 R4 -2 423% 932993 SHi70 902244 0.823 
2 161682 bite s5en> 4 348 74 330 7H 3687 1 
2b21arvG05eo2e 11469481 3 3.840 3.840 3.840 
293 sh 7591982 L144 260k 1 38932) H32932 3.932 
20454%7586e8e 12063081 4 34840 OS2670 324871 
205997579628 14239 37m 4 38962 932962 3.962 
2627457209" 12653891 4 34962 “39962 36°62 
2471447562989. 12963205 1 34962 434962 Be go2 
298491556749 43245004 4 30780° 436780 Se Oo 
2097845519389 1354e¢e4 4 35/780" 6332560 33780 
2204S 1546789 43639381 4 36566 436566 S566 
221761539292 14423784 4 365S7 E37 59% 345397 
2224171534764 1462368 1.4 34620 (233627, 36627 
2237875287982. 14723881 4 38629 8354627, 32627 
2244917527 7635200 139944. 4 33.237 (434234 S.23u 
225904520753 45302964 4 38284 "G3823i 34234 
301 1547606) 25533664 KKKKK KKEKKK KKK KK 
SO0Z2491550¢389 152096874 4 38356) 832653 3.356 
303711558420 14963364 4 36505 932505 Sie dey Ne 
804421565603 14666737 4 32656 B5e006 34658 
805107574680 1438 e350n 2 $6036 7852536 3,530 
SHON OP img! lhe t spat momeeiey p's Os fear ks a 36 Ow Oe D006 32506 
S07. Sos0e 3) Moo. wie SiO Osmo. SOG SMe ois: 
BUGSaetOS Veo seo owed B09) = 9S.Do/ 32057 
S0Ssn to04 2. 13207,50 173) Se eS Oo, Sof 
Sis E607. Ge. 1209 6) 1s a) SOC OLrS C06 B7566 
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11.2 SAMPLE OF DRILL HOLE DATA 


ID = The sample number, in sequential order, book-fashion, 
from the northwest corner of the study area 
Coordinates = Metres north and east from the SW corner of 
NE-33> 1-6-W2). 
F AVilag indicating ithe data source: /1=pit, S=drill hole? 
N The number of splits. 
All thicknesses are in metres; **** indicates missing data. 


iD SCoordas EP NeIhICKnesGeCoat et wCoat “Pt "Coal Pt Coal 


Xx vi gross net bottom top 
Mezo2d MoUSG 135) 383553235 e285 
Cmcrice, SUCHE Scho oa", CUNO LTS) 2s 4 
302226 3067 3 | HEX £EXK) £45% 
Sec oeun SU00 Ss 4 4, UG tee. 0 Oa40. t210.54 028270279" 02514047 
Se CO tO es eee Oe 0a 0 lOO oa oe ONO 2 a) 5 
GeecUet 206 ome OOo es Oe SUA 0. 24 «300 
Te Zeoe CIOONS Sh oO s on, BO seo. to 0,64 OS 2.04 
Cee 2 OOS oe S59 13.00 0650.15 0.35.0, 4202 822 
SEAS Zo 293) e824. 32 0St fac Omron 0,42 0.45 02547 047970 630 (0.6 
POGOe IO) Ste ec oS me WW. 2ZOP Oe TOO ol) 0.2 180s 76508507058 
eee 22h OS seer 2 One eo One Oe Oe oa. OG 
TAer VOOM 2 OCG so 3 tao ees Use7! OL B2. 0.67 1 tent eog 
PSO OCOD 4 Some ie Se OlO coy 1 O00) 24435 
eee cs ZOOL StS SOM SOO Oa o= Os FOU le UOmOn bores 1s 
POO l ce CO OE Sano SO OMS Tomson 22 lO onOme mie Z6 
BOS o2 22604 S22 cole SOP Wes he Oe 0 ee es 
POE COCe TOC SMO as SON see Seed Seu Sica 
homer 24 a 7B O4 3S. SSS ese Oa 30 01 SO on 0L 6 1. 1. 96 
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COOL ee Ogms LoS sO So 50 tl oO meen oo lO 2 OUe t.oG 
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APPENDIX 12 
SAMPLE OF COMPUTER PROGRAM OUTPUT 
12.1 OUTPUT FROM THE VARIOGRAM PROGRAM 
Veber 120. G RAM 
SAMPLE DATA FROM PIT 2A, AVERAGE VARIOGRAM 
WLTH -A -& ELD -OF 


90. DEGREES IN EACH DIRECTION 


DATA FILE = SAMPLE.DATA 
FORMAT = (i XteesGo no poms ok, 110K) (1k Fr es)) 
LAG IN METERS # ASO 0.0 Orme OMe VEO atone verete stats a teators 
NUMBER OF SAMPLES # 30 SAMPLE 
LOWER LIMIT FOR 2 # B73 40 SRAM es 
UPPER LIMIT FOR 2 # Bele dabale, (MP cies wie aie as 
GENERAL MEAN OF Z # 3.7572 360.000 
GENERAL VARIANCE OF Z # 0.0624 ry i eae 
DIST IN METRES N.PR DRIFT VARIOGRAM AVG DIST 
0 ---- 40 33 -0.046 0.0090 SHG 
40 ---- 80 47 -0.085 OVO 7 58.1 
80 ---- 120 48 -0.150 0.0261 96.7 
120 ---- 160 48 -0.143 02-0293 hetoh ey) 
160 ---- 200 29 -0.245 0.0478 180.0 
200 ---- 240 31 -0.289 0.0643 216.6 
240 ---- 280 34 -0.258 0.0546 259.5 
280 ---- 320 22 -0.339 0.0767 301.0 
320 ---- 360 24 -0.353 0.0842 38600 
360 ---- 400 oF -0.348 0.0889 379.04 
400 ---- 440 20 -0.419 0.1097 E2380 
440 ---- 480 15 -0.457 0.1366 459.0 
480 ---- 520 19 -0.470 0.1361 497.8 
520 ---- 560 14 -0.435 0.1086 542.2 
560 ---- 600 8 -0.438 0.1045 581.7 
600 ---- 640 10 -0.527 0.1589 617.6 
640 ---- 680 4 -0.617 0.11974 656.6 
680 ---- 720 2 -0.640 0.2067 686.4 
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12.2 OUTPUT FROM THE KRIGING PROGRAM 
DATA USE) (D 
DATA FILE USED = SAMPLE. DATA 
DAT FORMAT ANSE.) =aielxe oboe ik 2 UF load, 1k) 23k 1 Ok), (1k, F763) ) 
THE NUMBER OF SAMPLES IS: 30 
VEACR CO GR AGM.” MOD Ber iS 


THERE ARE 2 SPHERICAL MODELS 


COEFFICIENT= 0.0200 RANGE= DU a0 
COEFFICIENT= OFC a6 RANGE= 180.0 
NUGGET EFFECT= Ore SILL= 0.0358 
MAXIMUM RANGE=_ 180.0 


DOMAIN 


POLYGON FILE USED = PIT.2A 

POLYGON FORMAT USED = (7X,2(F16.5,1X)) 

AREA OF THE DOMAIN IS: 202309. 
AREA OF THE ENCLOSING RECTANGLE IS: 374986. 
EXTREMES OF THE ENCLOSING RECTANGLE ARE: 

XMIN= 1483.38 

XMAX= 1945.32 

YMIN= 841.66 

YMAX= 1653.42 


Poe Me ay ON Cea) 


BS lI MATEDOCELE oLZe. 1s: DES e) 3) 

(Eh) le) mph yA) glad Oy ah Bp 20.00 

COLUMNS= 28 

ROWS= 45 

MAXIMUM NUMBER OF ELEMENTS= 1260 
PROBABLE NUMBER OF ELEMENTS IN DOMAIN= 679 
EXTREMES OF THE ENCLOSING RECTANGLE ARE: 
XMIN= 1460.00 

AMAX= 1960.00 

YMIN= 820.00 

YMAX= 1660.00 

ACTUAL NUMBER OF CELLS IN THE DOMAIN IS: 503 AND 604 


MINOR WARNING: IER = 34 
1 DIGITS TESTEDAFOR ACCURACY < 
THES ACCURACY. TRS) tFATLED. 
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Regions Ue eee ks 


AREA OF DOMAIN IS: 
ZUZS0G00 
PERIMETER OF DOMAIN IS: 
na le ele) 
TOTAL NUMBER OF SIDES IS: 
78 


AVERAGE LENGTH OF A SIDE IS: 


27°36 


KRIGE RESULTS 


KRIGE Z VALUE IN THE DOMAIN: 


3.789875 
KRIGE VARIANCE: 
03002387 
THE SUM OF WEIGHTS: 
ORrooiee oS 
COVARIANCE OF THE DOMAIN: 
OF COn2797 
AVG COV OF SAMPLE TO DOMAIN: 
0.001234 
LAGRANGE PARAMETER: 
0.002344 


EXTENSION RESULTS 


VAR. 

.000001 
.000000 
.000000 
.000000 
.000000 


.000000 


ARITHMETIC MEAN OF Z IN THE DATA FILE: 


So Ow 
EXTENSION VARIANCE: 
0.002584 
AVG COV SAMPLE TO SAMPLE: 
0.003803 
COVARIANCE OF THE DOMAIN: 
O5.00) 2407 
AVG COV SAMPLE TO DOMAIN: 
0.001248 


0.000000 


0.000000 
0.000000 


aoe 

-000820 
sUO0OZT 
.000001 
.000136 
.000061 
.000054 


~000014 


.000136 
~000061 
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